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ABSTRACT 


In this thesis we treat the “small cell” problem encountered when building an attri- 
tion rate generator for large-scale manpower flow models, specifically for the USMC 
Officer Corps. Such models have a large number of low-inventory (i.e. small) personnel 
cells. This presents a dilemma: on one hand we want to preserve as much fidelity as 
possible in our work by preserving a great deal of detail in each cell; on the other hand 
our Statistical estimation techniques require larger cell sample sizes than intrinsically 
occur cell-by-cell in actual sample data. Our approach to producing stable attrition rates 
for such cells involves two efforts: (1) the aggregation of cells into groups that exhibit 
homogeneity of attrition behavior, and (11) the development of “shrinkage” estimation 
techniques for use in the individual groups. A heuristic algorithm is developed and 
tested to treat the aggregation problem. Empirical Bayes methods are developed to serve 
the multi-cell estimation requirements needed to preserve the fidelity. Cross validation 
techniques are used to verify these methods. 

The present work builds upon the results of previous studies; we integrate what was 


learned into a coherent package that is ready for use. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 


additional verification ıs at the risk of the user. 
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I. INTRODUCTION 


A. GENERAL 

The Officer Planning and Utility System (OPUS), a comprehensive and fully inte- 
grated manpower management system, is currently being implemented by the U.S. Ma- 
rine Corps (Decision System Associates, 1986). This system contains a set of 
computer-based manpower planning models and is used by the Officer Plans Section 
(MPP-30), Headguarters, U.S. Marine Corps, to produce several manpower planning 
documents. The system must be able to accurately predict personnel attrition, i.e., offi- 
cers leaving the service for purposes such as resignation, retirement, discharge, disability, 
or release. The forecasting of attrition is accomplished by the Marine Corps Officer 
Rate Projector (MCORP), developed by the Navy Personnel Research and Development 
Center (NPRDC), San Diego, California (NPRDC, 1985). 

The attrition rate generator developed by NPRDC calculates empirical attrition 
rates using historical data with user-defined weights and threshold parameters (Siegel, 
1983). This subjective input makes the current generator susceptible to unintentional 
misuse. 

In support of MCORP, Professor Robert R. Read of the Naval Postgraduate School 
has been working on the “small cell” problem: applying multiparameter statistical esti- 
mation schemes to estimating attrition when there is low personnel inventory, or small 
cells, which generally exhibit unstable empirical rates. 

A comment on terminology is in order. By attrition rate generator we mean meth- 
odology for estimating attrition probabilities for the various cells. The expression “em- 
pirical rates” refers to the ratio of leavers to inventory for each cell, unmodified by any 
information contained in “neighboring” cells. In contrast to this, the expression “em- 
pirical Bayes” refers to Bayes estimators whose prior parameters are estimated from 
data. 

Accurate forecasting of losses 1s extremely important to the manpower planner. 
Overestimating losses causes excess accessions, promotion delays, underutilized person- 
nel and increased costs, while underestimation causes personnel shortages and decreased 
readiness. The problem is compounded in that all but a few accessions must start at the 
bottom, 1.e., Second Lieutenant, and work their way up to the higher ranks only after 


manv vears of service. For example, if a shortage of Lieutenant Colonels arises, it can 


only be remedied by promoting more Majors, which has a rippling effect down the rank 


structure. 


B. BACKGROUND 

There have been seven Master's theses over the past four vears which have studied 
various aspects of the attrition estimation problem. A concise summary of these works 
is given by Read (NPS Report NPS55-88-006, 1988, pp.16-23). These studies can be 
grouped into three general areas: shrinkage methods, cell aggregation and peripheral 
studies. 

The application of a shrinkage method begins by identifying a number of personnel 
inventory cells, followed by the development of the empirical rates for individual cells 
and a weighted grand mean of these empirical rates. The final estimate for a cell is a 
convex combination of its empirical rate and the grand mean. There are numerous 
methods for accomplishing this, several of which have been applied in previous studies. 

Tucker was the first to investigate the application of these methods to attrition es- 
timation. He compared traditional estimators to the James-Stein estimator and the 
minimax estimator for a few selected grades and occupational fields. His results gave 
strong support to the James-Stein estimator; minimax was discarded as being too con- 
servative for small cell use. However, there remained pockets of instability for which 
goodness-of-fit tests failed. (Tucker, 1985) 

Following Tucker was Robinson, who introduced the Efron-Morris limited trans- 
lation shrinkage alternative to augment the James-Stein estimator. These methods were 
evaluated with a broader set of test cases. Robinson was able to confirm Tucker's re- 
sults, but the limited translation option failed to provide any consistent relief in the un- 
stable areas. (Robinson, 1986) 

The final application of shrinkage methods to estimating officer attrition rates was 
undertaken by Dickinson. He applied the previously used methods and an empirical 
Bayes estimator to a new and refined data base. Improved results were obtained, but 
the instability remained. Dickinson also performed some exploratory side studies dealing 
with the Freman-Tukey transform and the use of empirical Bayes methods that allow 
non-uniform shrinkage, both of which provided the impetus for the present study. 
(Dickinson, 1988) 

These three studies used ad hoc methods to deal with the second general area of 
study--cell aggregation. Aggregation of cells with low personnel inventory into sets of 


cells, often of larger inventory, is required when applying these shrinkage methods. The 


desire is to use cells which exhibit similar attrition behavior. Two previous studies have 
investigated this area. 

Amin Elseramegy used the Classification and Regression Trees (CART) program, 
which at the time was newly acquired by the Naval Postgraduate School, in an attempt 
to form aggregates of cells that exhibited homogeneous attrition behavior. Several dif- 
ficulties in using this program were encountered, e.g., because of insufficient memory 
allocation he found it necessary to partition the data base into nine sets and applv 
CART to each. The resulting aggregations were generally unusable. (Amin Elseramegy, 
1985) 

Major breakthroughs in cell aggregation were made by Larsen. He applied a hi- 
erarchical clustering algorithm to the new data base. The resulting rules for building 
aggregates are well defined and especially viable from an intuitive point of view. 
Larsen's work provides the framework for the cell aggregation method developed in 
Chapter II of this thesis. (Larsen, 1987) 

The remaining two theses of the seven were peripheral studies which applied alter- 
nate methods to attrition estimation. Hogan attempted multi-vear forecasting using 
exponential smoothing; the smoothing constants were rather unusual and extreme and 
his results inconsistent (Hogan, 1986). Yacin applied logistic regression in the attempt 
to develop an attrition rate scheme; the onlv new results were the identification of some 
areas that exhibited similar attrition behavior (Yacin, 1987). 

This thesis is the first to integrate the two main areas of study. Whereas previous 
studies of shrinkage methods have used ad hoc aggregation schemes, we now combine 
the implementation of a defensible aggregation method with empirical Bayes estimators. 
Moreover, these are applied to a larger and more refined data base. The results have 
been quite promising in that we have achieved greater stability in attrition rate esti- 
mation; we have defined guidelines for a heuristically appealing aggregation scheme; and 
we have acquired an increased understanding of the data base and developed more effi- 


cient wavs to use it. 


C. ORGANIZATION 

The remainder of this introductorv chapter provides a more detailed description of 
the small cell problem and the data base. The aggregation problem is discussed and the 
proposed aggregation method is presented in Chapter II. 

The shrinkage estimation methods, generally classified as empirical Bayes type esti- 


mators, are described in Chapter III. Several variations are presented to allow 


comparison and to gain further insight into their performance. Testing of these methods 
is important but for practical purposes must be carried out using sampling methods. 
The rationale used to select test cases, the cross validation techniques, and the measures 
of effectiveness used to evaluate the results are discussed in Chapter IV. A discussion 
of the results of the cross validation 1s also included in this chapter. 

Finallv, conclusions and recommendations based on these results are contained in 
Chapter V. 


D. SMALL CELL PROBLEM 

Marine Corps officers can be classified and thus partitioned by several attributes. 
The major partitioning of officers is bv grade, vears commissioned service (YCS), and 
militarv occupational specialtv (MOS). Grade describes the position an officer holds in 
the service. The numbers of officers in each higher grade have a pvramid structure, 1.e., 
there are more officers in the lower grades than in the higher grades. YCS is the total 
number of vears served since becoming a commissioned officer. There is a strong cor- 
relation between grade and YCS since an officer generally moves up the grade structure 
as he gains in YCS. MOS is a four-digit code identifying the specific skill for which a 
Marine is trained. MOS need not remain constant over an officer’s career, although 
most changes in MOS occur in the early vears of commissioned service. An officer has 
a single primary MOS, however as he develops new job skills he may be assigned one 
or more additional MOSs. 

For many purposes, partitioning by grade, YCS and MOS 1s sufficient. However in 
Some applications additional refinement by service component, commissioning source, 
sex, race or education level may be necessary. Service component consists of three cat- 
egories: regular officers, reserve officers, and reserves who have augmented to become 
regulars. It is strongly correlated with commissioning source, i.e., an officer receives a 
regular or reserve commission depending upon the commissioning source. Both affect 
an officer's initial service obligation, which is generally three to five years (except 
aviators, whose obligation ıs dependent upon the amount of flight training). Officers 
who receive a reserve commission normally serve three to four years active duty (except 
aviators), by the end of which thev must have either augmented into the regular force 
or are then separated from active service. 

These cross-classifications may be viewed as breaking the officer population into a 
multidimensional array, with each specific intersection of the classifications called a cell. 


The total number of possible cells is quite large, on the order of 10%. Many of these cells 


are structurally infeasible in that no officer could possibly fit the cell characteristics, e.g., 
there are no Majors with two years commissioned service. The total officer inventory 
of approximately 20,000 officers is partitioned by the remaining feasible cells, some cells 
have inventory as large as 50, however most have less than five. An officer's character- 
istics are dynamic, 1.e., as an officer moves through the grade, YCS and MOS structure 
he moves from one cell to another. As a result, the inventory of the feasible cells is also 
dynamic and fluctuates between zero and low inventory (less than five) over time. 

These sparselv populated cells have very unstable empirical attrition rates. For ex- 
ample, a cell whose inventorv is two officers of which one leaves the service during a 
given time period yields a 50% empirical attrition rate, whereas a cell whose inventory 
is one officer who remains in service during the same time period yields a 0% attrition 
rate. Itis obvious that neither of these empirical estimations provides a usable attrition 
rate. Furthermore, these two rates could change dramatically during the next time pe- 
riod, tvpifving their instability. 

Even when more modern estimation techniques (e.g. shrinkage) are applied, these 
small cells can still create statistical instability, thereby producing intolerably variable 
attrition rate estimates. The problem then is how to deal with these low inventory cells, 


or “small cells” in order to achieve stabilitv. 


E. DATA BASE 

In this thesis we benefit from a refined data base compiled by NPRDC and made 
available to the Naval Postgraduate School in 1987. This data base, used bv Larsen in 
his aggregation work (Larsen, 1987), was not available for the previous estimation 
studies at NPS. 

The new data base provides more detailed information about the officer population. 
The grade structure now allows separation of Limited Duty Officers (LDO) as well as 
Warrant Officers (WO) from unrestricted officers. Officers who have failed selection to 
the next higher grade can also be identified. YCS is listed instead of length of service 
(LOS), which became ambiguous when dealing with officers who have prior enlisted 
service. MOS can now be broken out completely into 236 MOSs or summarized by the 
39 occupational fields. Service component and commissioning source are both new 
categories. Other new categories that are not considered here are education level, race, 
additional MOSs, and military schools completed. Larsen gives a complete description 


of the classifications (Larsen, 1987, pp.66-82). 


The data base also allows attrition to be broken out by retirement, release, dis- 
charge, resignation, etc., but for our immediate purpose we are only concerned with the 
total number of losses for anv reason. 

This refinement of the data presents a dichotomy: we can now break the data into 
more definitive cells to search for homogeneous attrition behavior and stability in esti- 
mation, but this leads to an even greater number of low inventory cells. 

The new data base contains ten years of inventory and attrition data from the period 
1977-1986, a significant improvement from the previous seven year data base covering 
the period 1977-1983. The inventory data is now obtained from quarterly vice yearly 
snap-shots of the officer population. The attrition data is annualized, 1.e., the attrition 
count for a cell reflects the number of personnel who leave the service at any time during 
the year. Attritions are credited to the cell which the officer occupies at the time he 
leaves. 

Two problems arise from this quarterly versus annualized data. First, it is possible 
for a cell to record zero inventory via the snap-shots, yet be credited with one or more 
attritions. To avoid this situation, the cell inventory used in all calculations is defined 
to be the larger of the inventory and the attrition count. This ensures that the inventory 
for a cell is as least as large as its recorded number of leavers. (This override occurs in- 
frequently; a more sophisticated treatment would require significant model enhance- 
ment.) Second, to use the inventory and attrition data together we must divide the 
inventory data by four. This poses a philosophical problem when invoking a binomial 
model: the sample size may not be integral. However, for our application the usual 


mean and variance formulas are usable and can still serve in the interpolative sense. 


o 


H. CELL AGGREGATION 
A. GENERAL 


The aggregation problem takes on new meaning with the use of shrinkage estima- 
tors. Originally, aggregation had only one concern: how to pool cells together into a 
single cell in order to meet a user-defined minimum inventory threshold. This single 
aggregated cell was then used to determine the attrition rate estimate for the original, 
unaggregated cell. In this way an estimated value for a cell is obtained by using the 
grand mean for many cells. 

The empirical Bayes multiparameter estimation techniques provide a way to com- 
promise, using both the stability of a grand mean and the specific information of an in- 
dividual cell. Now we pool cells together and obtain a number of cells that meet the 
user-defined minimum inventory threshold. It is important to note that we should be 
able to use a lower inventory threshold with empirical Bayes, thus retaining individual 
cell behavior to a greater extent. It is also important to use cells with homogeneous 


attrition behavior in the aggregation process. 


B. BACKGROUND 

The aggregation method currently used by MCORP 1s called the Small Cell Override 
Methodology (NPRDC, 1985, Appendix H). It is used to solve the original aggregation 
problem, 1.e., if a cell is below the user-defined threshold, then cells are adjoined to the 
original cell until the threshold is met. The process for selecting cells for adjunction 1s 
rather crude, and large-scale with only a few levels (prior to using the entire officer 
corps). The attrition rate estimate for the original cell is the empirical rate from this 
aggregated cell. 

To begin the process, the user defines a cell for which an attrition rate estimate 1s 
required bv grade, YCS and MOS. The user also defines the minimum cell inventory 
threshold (and other parameters which are not relevant here). If the cell he identifies 
meets the threshold, no aggregation is required and the empirical attrition rate 1s deter- 
mined. If the cell is below the threshold, additional cells must be added until the 
threshold is met. 

This search for additional cells occurs by expanding by YCS and MOS, with grade 
remaining fixed throughout. Expanding in this sense means changing the YCS or MOS 


parameter to identify the additional cells to be added to the original cell. Initially, the 


single cell is expanded by YCS. For example, if the original cell's grade/YCS/MOS was 
Capt/7/0802, the cells identified by Capt/6/0802, Capt/8/0802, etc., are added sequen- 
tially until the threshold is met. This YCS expansion has an upper bound at the 20 YCS 
point; an obvious boundary for attrition behavior due to retirement eligibility. If the 
original cell's YCS is above 20, then 20 would serve as the lower YCS bound. 

If the threshold is not met after the YCS expansion, the override method starts over 
with the original cell and expands by MOS. Each MOS belongs to one of nine MOS 
groups which are defined along traditional Marine Corps functional areas, e.g., all heli- 
copter pilot MOSs are grouped together as are all combat support MOSs. MOS ex- 
pansion adds those cells identified by the MOSs in the same MOS group as the MOS 
of the original cell for the original YCS and grade. If the threshold is not met, all the 
MOSs in the MOS group are expanded by YCS in the same manner as the YCS expan- 
sion discussed previously. 

If this MOS group and YCS expansion is unsuccessful, the override method starts 
over with the original cell and expands by all MOSs for the original grade and YCS. If 
necessary, all the MOSs are expanded by YCS as before. 

Cell aggregation using this expansion method can potentially include all MOSs and 
YCS bounded only at the 20 year point. The desire to aggregate using cells with ho- 
mogeneous attrition behavior is obviously compromised. Larsen provides a more com- 
prehensive description of the current method (Larsen, 1987, pp.16-22). 

Larsen examined attrition behavior in the MOS and YCS structure. He applied a 
hierarchical clustering algorithm in an attempt to find MOSs and YCSs that displayed 
homogeneous attrition behavior. He confirmed the belief that YCS is an important 
factor. The YCS expansion bounds he proposed reflect points at which officers reach 
the end of their initial service obligation as well as when they are eligible for retirement, 
which makes them especiallv viable from an intuitive point of view. Larsen also found 
that some MOSs did not cluster strictly by functional areas. This was especiallv signif- 
icant in the aviation community. Whereas the previous data base allowed aviators to 
be considered only as one occupational field, the refined MOS information was able to 
identify six distinct homogeneous groups of aviators. 

Larsen uses these results to define more refined MOS groups and YCS boundaries. 
To avoid the giant expansion leap from MOS group to all MOSs, he proposed a hier- 
archy of small MOS groups, large MOS groups and major MOS groups developed by 
observing which MOSs tend to exhibit similar attrition behavior. Homogeneity is 


greatest within the MOS group, and becomes successively worse as we move to the large 


MOS group and then the major MOS group. Each MOS is assigned to a small MOS 
group. Small MOS groups combine to make a large MOS group, and large MOS groups 
combine to make a major MOS group. 

Each small MOS group is assigned a set of YCS expansion bounds. Due to the 
different attrition behavior of the small MOS groups with respect to YCS, three different 
sets of YCS expansion bounds are proposed. 

Initial expansion is by YCS within the specified boundaries, with grade and MOS 
held constant. If more expansion 1s required, we retain this aggregated cell and expand 
by small MOS group for the original grade and YCS. If the aggregated cell is still below 
the threshold, the MOSs in the small MOS group are expanded by YCS. Subsequent 
expansion to large MOS group and YCS, and major MOS group and YCS is accom- 
plished until the threshold cell inventory 1s met. 

Unlike the current expansion method, expansion using Larsen’s proposed method 
will not cross defined MOS groups or YCS bounds to ever include all MOSs and YCSs 
bounded only at the 20 vear point. Larsen provides a more detailed description of his 


recommended expansion rules (Larsen, 1987, pp.45-61). 


C. EXPANSION METHOD 

We now address the methods used to obtain the cells required for use with empirical 
Bayes estimation techniques. Expansion continues to mean finding more cells to be 
used, however we no longer simply add these cells to the original cell to form a single 
aggregated cell. The cells identified bv the expansion process are now aggregated to- 
gether to produce a number of cells. After the discussion of the expansion process in 
this section, an actual aggregation scheme is introduced in the next section. 

To begin the estimation process, the manpower planner defines a specific cell by 
grade, YCS and MOS. The attributes service component and commissioning source are 
also included as possible cell descriptors for the purposes of this study. All other de- 
scriptors listed in the data base--sex, education level, additional MOSS, race and military 
schools--are ignored. Loss types are considered as a combined total, i.e., in this study 
we do not discriminate among the various types of losses. The first three user-defined 
descriptors--grade, YCS and MOS--are single-value inputs. The last two descriptors, 
service component and commissioning source, can be single values, or either one of them 
can be treated as a vector of values for each single cell. This vector is collapsed (total 
the components) during the aggregation process, i.e., all records which meet any of the 


vector's values are included in the same cell. As in the previously described expansion 


methods, only YCS and MOS change during expansion, the remaining cell descriptors 
remain constant. 

To use shrinkage techniques, the amount of expansion required not “only depends 
upon the minimum cell inventory threshold but also upon a new input parameter: the 
threshold number of cells. These two parameters are denoted: 


l. T, - cell inventory threshold. The minimum average inventory for a cell obtained 
by averaging the cell inventory over the ten years of data. 


2. K, - threshold number of cells. The minimum number of aggregated cells whose 
inventory exceeds 7). These aggregated cells are the input cells for the empirical 
Baves technigues. 

For example, if 7, = 5.0 and K, = 10, the expansion algorithm continues until at 
least ten aggregated cells, each with average inventory 5.0 or larger are obtained. Since 
we are concerned primarily with the small cell problem the values of 7, and K, used are 
selected to range from five to 30. It is also presumed that T, is less than or equal to 
K,. These threshold values can certainly exceed 30 for other applications, however the 
resulting cells are not considered small and their attrition behavior most likely would not 
be as unstable, therefore not requiring special attention. 

Prior to explaining the expansion process, we first define the MOS groups and YCS 
bounds. We have adopted much of Larsen’s work in this area; many of the changes are 
minor but are necessary for implementation purposes. 

The general idea of a hierarchy of MOS groups is repeated, as shown in Table 1. 
Each MOS belongs to a small MOS group, a large MOS group and a major MOS group. 
Listed are 14 small MOS groups, which combine to make six large MOS groups, which 
combine to make four major MOS groups. For example, small MOS groups one and 
two form large MOS group one, and small MOS groups three through six form large 
MOS group two. Large MOS groups one and two, which collectively contain small 
MOS groups one through six, make up major MOS group one. Major MOS group one 
contains only ground MOSs, and major MOS group two contains only aviation MOSs. 
Major MOS groups three and four are special cases as discussed below. 

A subjective decision was made to keep the ground MOSs in groups defined along 
the more traditional functional areas. This is reflected in small MOS groups one through 
six. For estimation purposes it is advantageous if the cell inventories are not too vari- 
able in size (Carter and Rolph, 1974, p.882). It is also desirable to avoid having too 
many MOSs in each small MOS group. This allows the expansion to occur more grad- 


ually, and is especially important for small values of 7, and K,. Asa result, MOS 0302 
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Table 1. MOS GROUPS 


Group Name MOSs MOS MOS | MOS 
Group | Group | Group 

L le 
[Combat Support [0802 1302 1802 1803 | 2 
| Combat Servicel | 0180 0202 2502 2602 | 3 | 
Combat Service 2 | 3415 4002 4302 5803 | 4 | 
Combat Logistics | 0402 3002 3060 3502 6002 | 5 | 
[Air Control | 7204 7208 7210 7320 | 6 | 


7501 7511 7522 7542 7543 
Fixed Wing Pilots 7545 7576 i 3 


F-18 Pilots 192198 (923 


7556 7557 7562 7564 7565 2 
Rotarv Wine Pilots + 7566 7587 SR 


7508 7509 7563 7581 7583 
“Naval Flight Officers + 7584 7585 7586 7588 > 


0101 0201 0301 0401 0801 
| 1301 1801 2501 2601 300] 
Basic Ground 3401 3501 4001 4301 4401 
5801 6001 7201 7301 9901 

7580 7597 7598 7599 
o 7500 7510 7520 7540 7550 


1560 7575 


1302 E: 
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(infantry) is placed alone in a small MOS group. This MOS contains approximately 
15% of the total officer population, and therefore its respective cells normally contain 
large inventory. The MOSs in small MOS group two also contain fairly large inventory, 
therefore are grouped together and their first expansion is with MOS 0302. The 
remaining ground MOSs in small MOS groups three through six have similarly small 
Inventory. 

The aviation small MOS groups (seven through ten) remain relatively unchanged 
from Larsen's recommendations. MOS 7564 (CH-53 pilot), was removed from a ground 
MOS group and added to small MOS group nine, which reflects its functional area. 
MOSs 7551 (C-9 pilot), 7552 (TC-4C pilot), 7555 (UC-12B pilot) and 7559 (CT-39 pilot) 


were deleted since they are not primary MOSs. MOS 7530 (basic pilot VMFA (F-4)) 
was deleted since it is not a current MOS. (MCO P1200.7G, 1988) 

Officers who have not acquired sufficient schooling or field experience to qualify for 
a primary MOS listed in small MOS groups one through ten are gathered together as 
basic officers or students in small MOS groups 11-13. These officers are generally sec- 
ond lieutenants or junior first lieutenants with three or fewer YCS. They are disregarded 
for the remainder of the study because their attrition rates are extremely low; probably 
because none of the officers in these groups have reached the end of their initial 
obligations. 

MOS 4402 (lawyers) is considered a special case and is not addressed in this study. 

All MOSs listed in Table 1 are primary MOSs for unrestricted officers as listed in 
the current Military Occupational Specialties Manual (MCO P1200.7G, 1988). It would 
be a logical and relatively simple extension of this table to create additional groups 
containing LDO and WO MOSs. These grades are not considered in this study and 
therefore their respective MOSs are excluded from the table. 

Several of these seemingly ad hoc decisions to alter Larsen’s recommended MOS 
groups are due to the YCS expansion bounds shown in Table 2. Every effort was made 
to group MOSs with similar YCS expansion bounds to allow for feasible implementation 
of the expansion algorithm. This is especially applicable when expanding to large and 


major MOS groups. 


Table 2. YCS EXPANSION BOUNDS 


Small 
MOS Group MOS Bounded YCS Groups 
Groups 
"ing 
Fixed Wing Pilots, F-15 Pilots, 7, 8, 14 (1-6, 8-19) (7) (20-25) (26) 


Lawv CTS 


Rotarv Wing Pilots, Naval Flight 
EES D Se (1-5, 8-19) (6,7) (20-25) (26) 


All Others 1-6, 11-13 | (1-3, 6-19) (4,5) (20-25) (26) 
















The YCS expansion bounds reflect the maximum expansion allowed from the initial 
YCS defined by the user. For example, if the original cell's grade/YCS/MOS is 
Capt/9/7501, we see from Table 1 that this MOS belongs to small MOS group seven. 


Thus its YCS expansion bounds are listed on the first line of Table 2. The value of nine 


1⁄2 


for YCS falls in the first YCS range, thus we could expand using all YCSs from one 
through 19, excluding seven. If the YCS for this original cell had been seven, no YCS 
expansion would be allowed. 


These YCS expansion bounds are used with the MOS groups to define the addi- 


tional cells which can be used with the original cell to obtain the required number of 


cells, K, each with minimum average inventory, Tọ The expansion stages are: 


1. 


Stage 1 - Locate the small MOS group which contains the user-defined MOS. The 
initial cells are those specified by the MOSS in this group for the user-defined YCS, 
grade, service component and commissioning source (grade, service component and 
commissioning source remain fixed throughout the expansion process and thus are 
not repeated). These cells are aggregated to obtain cells with average inventorv 
greater than or equal to To. After aggregation, if the number of cells is greater than 
Ky, stop, otherwise go to Stage 2. 


Biases - Expand Oy incrementine YCS (YCS-1, YES+ 1, YCS-2, YES 2, etc.) 
within the YCS bounds listed in Table 2 for all MOSs in the small MOS group. 
After each YCS increment, aggregate the cells to obtain cells with average inven- 
tory greater than or equal to 7,. After aggregation, ifthe number of cells is greater 
than A,, stop, otherwise continue to increment by YCS. If the YCS bounds are 
reached before obtaining enough aggregated cells, retain the cells identified in 
Stages 1 and 2 and go to Stage 3. 


Stage 3 - Expand to the large MOS group for the single user-defined YCS. Ag- 
gregate the cells to obtain cells with average inventory greater than or equal to 7). 
After aggregation, if the number of cells is greater than K,, stop, otherwise go to 
Stage 4. 


Stage 4 - Expand bv incrementing YCS for the large MOS group. After each YCS 
increment, aggregate the cells to obtain cells with average inventory greater than 
or equal to 7). After aggregation, if the number of cells is greater than Kj, stop, 
otherwise continue to increment by YCS. If the YCS bounds are reached before 
obtaining enough aggregated cells, retain the cells identified in Stages 1 through 4 
and go to Stage 5. Á 


. Stage 5 - Expand to the major MOS group for the single user-defined YCS. Ag- 


gregate the cells to obtain cells with average inventory greater than or equal to 7.. 
After aggregation, if the number of cells is greater than K,, stop, otherwise go to 
Stage 6. 


. Stage 6 - Expand by incrementing YCS for the major MOS group. After each Y CS 


increment, aggregate the cells to obtain cells with average inventory greater than 
or equal to 7). After aggregation, if the number of cells is greater than Aj, stop, 
otherwise continue to increment by YCS. If the YCS bounds are reached before 
obtaining enough aggregated cells, stop. No more expansion is allowed. Inform 
the user that the thresholds are unattainable. Do not cross any major MOS group 
or YCS bounds. 


Two points about the expansion process are emphasized. First, we retain the cells 


identified by all previous stages as we progress to the next stage. As stated before, the 


degree of homogeneity decreases as we move from small to large to major MOS groups. 


Thus we want to locate as manv cells from the small MOS group as possible before we 
expand to the large MOS group, and then locate as many cells from the large MOS 
group as possible before expanding to the major MOS group. The YCS expansion for 
each group may be different, e.g., the small MOS group may be expanded by all YCSs 
within the given YCS range, but the large MOS group may only be expanded by a few 
YCSs before the thresholds are attained. 

The second point is that, When aggregating cells, any aggregation that was per- 
formed previously is discarded and all cells currently identified are pooled and made 
available for aggregation. This affords the aggregation algorithm greater flexibility and 
could create more aggregated cells than if the aggregated cells from previous stages were 
left intact, thereby keeping the amount of expansion to a minimum. 

The aviation small MOS groups (seven through ten) make up the only major MOS 
group (two) that contains different YCS bounds, i.e., small MOS groups seven and eight 
have different YCS expansion with regard to year six than do small MOS groups nine 
and ten. To implement the expansion algorithm in a computer program, this difference 
is overcome by using the YCS bounds for the original user-defined MOS. For example, 
suppose MOS 7501 from small MOS group seven is the original MOS. If MOS expan- 
sion continues into the major MOS group, the MOSs in small MOS groups nine and ten 
would follow small MOS group seven's YCS expansion bounds. 

In summary, this method of grouping MOSs should provide greater homogeneity 
among cells which are used in estimating attrition rates. Unlike the current method, 
ground and aviation MOSs are never used together. The YCS bounds provide a logical 
and effective way to treat periods of different attrition behavior. However, the greater 
the expansion the less homogeneous the cells become, which should be kept foremost in 


mind when setting the threshold parameters. 


D. AGGREGATION METHOD 

While the expansion steps are being undertaken in order to achieve the threshold 
levels specified by the user, those cells with inventory less than T, must be gathered up 
into larger, aggregated cells whose combined inventory exceeds T}. In order to limit the 
expansion to as few additional MOSs and YCSs as possible, we desire to maximize the 
number of aggregated cells obtained at any stage of the expansion. 

The term maximization suggests the possible use of linear programming (LP). While 
an LP would ensure maximization, this would not be a trivial problem to solve, i.e., the 


LP relaxation would almost certainly fractionate cells, using their inventory in more than 


N 


one aggregated cell. This is not allowed since a cell may be assigned intact to only one 
aggregated cell. Thus an integer LP would be reguired which would typically contain 
500 or more integer variables. This method would not be expedient in terms of computer 
usage, especially considering the potential number of integer LPs that may have to be 
solved for a single estimation cycle. 

While we are trying to maximize the number of aggregated cells, it would be satis- 
factory to obtain close to the maximum if we could preclude the expense in computer 
time required by an integer LP. For this reason, a heuristic “greedy” algorithm was de- 
veloped. Complete descriptions of the heuristic algorithm and the LP formulation are 
contained in Appendix A. The performance of this heuristic 1s discussed along with the 
results of the empirical Bayes methods in Chapter IV. 


III. ESTIMATION METHODS 


A. GENERAL 
Once the cell aggregation phase is completed, we begin the attrition rate estimation 


process. The following notation is used to define the cell data 


K 
T 


number of cells 1) 
number of years of data. ( 


Then fori = l,- K andi = UIM! 


` 


Wé 
Y; (1) 


inventory of cell / in year t 
number of attritions in cell i in year 1. 


(2) 


The cell data is assumed to be independent binomial, Le, Filz Bin(N; (t), p;) . A 
success is defined to be an attrition, 1.e., an officer from that cell leaves the service during 
the vear. The empirical attrition rate for cell í is given by the Maximum Likelihood 
Estimator (MLE) 


EO x 
Pi = D (3) 


I 


This estimate of p works well for cells with large inventory, but not those with small in- 
ventory, which is most often the case in our application. 

The MLE has been shown to be dominated by shrinkage methods for K > 3 (Carter 
and Rolph, 1974; Efron and Morris, 1975; Casella, 1985). These methods find a grand 
mean or central attrition rate for the group of cells and a shrinkage factor for each cell. 
Previous theses have primarilv used a common shrinkage factor for all cells; we now al- 
low this shrinkage factor to vary from cell to cell. Each cell’s MLE is shrunk towards 
the central rate by its shrinkage factor. In this way, attrition information from one cell 
“spills over” into other cells. 

The shrinkage methods are developed under the theoretical assumption that the data 
is normally distributed. Most of the previous studies using empirical Bayes methods 


have used independent normal data with constant variance (Efron and Morris, 1972, 


1973, 1975; Dickinson, 1988). Some applications have used binomial data, using a 
transformation to make it behave more like normal data (Carter and Rolph, 1974; Efron 
and Morris, 1975). In Carter and Rolph's estimation of fire alarm probabilities, trans- 
formation of the binomial data did not have a large effect on the results (Carter and 
Rolph, 1974). Our application allows us to investigate the impact of the transformation 
when applied with more extreme values of p. 

Six variations of the empirical Bayes method are applied to the attrition rate esti- 
mation problem. The first four are similar in that they use the same iterative procedure 
to compute the amount of shrinkage for each cell. Of these four, two are on the trans- 
formed scale and two on the original scale. Each scale includes two methods of com- 
puting the cell variance: one method where the variance is time dependent and the other 
where it is time independent. The two variance calculations, if they produce like results, 
provide supporting evidence for the assumption that the data is independent and iden- 
tically distributed over time. This assumption is certainly questionable, since an officer 
who remains in a given MOS will move through the YCS and grade cell structure in a 
predictable manner. As a result, variance that is constant in time (time independent) 
may not perform as well as one that allows for time variation. The fifth method uses a 
different iterative procedure to determine the amount of shrinkage and is addressed 
separately in paragraph III.C.. The final method breaks the cell data into its vector 
components (service component or commissioning source) before shrinkage technigues 


are applied and is addressed in paragraph IIILD.. 


B. EMPIRICAL BAYES 
1. Transformed Scale 
We begin our application of empirical Baves methods on the transformed scale 
in an effort to overcome some of the weaknesses in our assumptions. The transforma- 
tion we use is the Freeman-Tukev transform, a modification of the basic arcsin trans- 
formation for binomial data. Its purpose is to stabilize the variance at one and make 


the data behave more like normal random variables. The form used is 


EE `` 2 4 l 
Al) = > V M (t) + 5 EIE -1) + arcsin( 2 KOT =) , (4) 
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Now, let 


X. 
XT;(t) = O) for t = 1,...,T; (5) 
NORS 


except when M, (r) = 0 (no inventory in year t), in which case XT, (t) does not exist and 


we reduce T; by one. The time average of the transformed values for cell í becomes 


XTB = TF LATO (6) 
t 


We now need to compute the variance of these time averages. Two methods 
are used: the first calculation is time dependent, 1.e., the variance changes over time, the 
second is time independent. 

The transform stabilizes the variance at one for large values of n and non- 
extreme values of p. These requirements on n and p are often violated in our application, 
therefore we have many combinations of n and p for which the variance is less than one. 
Dickinson was able to discover an interpolative formula which provides a good approx- 
imation for the variance of the transformed values, X, (£), for small values of n and p, and 


K2 3 (Dickinson, 1988, pp.8-11). This variance is given by 
Var(X,()) = muintl, MA; (1) (7) 


where V(X, (t)) is found by solving 


l 


HK) = a(X (AW) +0)? H+ C-1)? (8) 


with E 
Ga /N(D5(3) (9) 
and 
a = 1.6835 b = —.8934 be 9881, (10) 


Equation (8) obviously breaks down if X;(r)+ C< 1. When this occurs, we set 


X, (1) + C = 1.001 and continue. The effect is to use a small but positive variance. The 


value of one in Equation (7) dominates for about X, (1) + C > 2.2. The variance of the 
time average is then 


sd LD Var(X, (p) 
Var(XTB) = ` 2 Var(XT (1) = E O NOS (11) 
The second method of computing the variance is the more familiar one. Con- 


tinuing from Equation (6), the variance of the transformed values is given by 


Var(XT) = FFT, (XT, (À) — XTB} . (12) 


The variance of their average is therefore 


Var(XTB) = + Var(XT). (13) 
i 
Regardless of which variance calculation we use, the same iterative algorithm is 
used to determine the empirical Baves estimate for each cell. This estimate, XEB,, is 
found by solving 
A Var(XTB;) 


er (14) 


ren) EE TAĦT B) 


where XEB,, XTB, and Var( XTB) are cell specific, XBB is the (weighted) grand mean or 
central attrition rate, and A is the variance of the prior distribution of the cell means. 
These latter two values must be estimated simultaneously using the following iterative 
algorithm. 


We initialize the algorithm with 4 = 0 and store the previous value of A by 


Now compute the (weighted) grand mean, XBB. Let 


l 
Ses NN. 16 
" A + Var( XTB) (16) 


and 





7 K (17) 
Ð E 
ja 
Then 
K 
XBB = ) yXTB,. (18) 
i=] 
The updated value of A becomes 
K 
K= 1— ) a(XTB;— XBB) 
A- A4- EL. (19) 


as (XTB,— XBB) 


K 
=] 


* 


t 


If 4 <0, set 4=0 and exit. This represents the case when there is 100% shrinkage to- 
ward the grand mean. If 4> 0, then check | 4 — 4,| <e (e.g., e = .0001). If false, return 
to Equation (15) for another iteration. If true, the iterations have converged. Exit with 
the current values of 4 and XBB for use in Equation (14) to solve for the ZER, 

Close study of Equation (14) shows that the amount of shrinkage changes from 
cell to cell since the variance terms are generally not equal. Specifically, cells with higher 
variance are shrunk more than those with lower variance. In addition, if A is small the 
shrinkage is greater towards XBB. As A = oo, the shrinkage is minimal and the indi- 
vidual cell means dominate. 

Once the XEB, are determined, these values must be transformed back to the 


original scale. We use 


A 


Pi = (1 +sin(XEB)). (20) 
2. Original Scale 


We return to the assumption of binomial data for original scale calculations. 


As in the transformed scale, two methods to calculate the variance are used. We begin 
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with 
A Y; 
AI) = pt) = ko. (21) 


As before, if N, (1) = 0 (no inventory in year 1), XT, (1) does not exist and we reduce T, 
by one. This leads to the time average for cell í as 


= de = rn 
XTB = FLAT = F ÊO: (22) 

The variance calculation which is time dependent, i.e., changes over time, is 
given by 


ME” P; (Ð(1 = P; (1) (23) 


Var(XTB) = 2 Var(XT; (1)) = CNO 


t 


We return to Equation (15) with these variance values to perform the iterative 
algorithm for finding the empirical Baves estimate, XEB,, given by Equation (14). Since 
we are already in the original scale, the transformation given in Equation (20) 1s ignored, 
Re p, = XEB. 

A problem arises while performing the iterations if a cell has Y, (+) =0 Y z (zero 
attrition for every year). In this case, the variance given by Equation (23) equals zero. 
When this value is used in Equation (16), the formula for «, becomes undefined. We 
resolve this problem using the Laplace Law of Succession. Assume that 
Y,(t) ~ Bin(N, (1), p) and let 


o O+! < NO- YO 


AGE dë 24 
Pi Ee C Nj() +1 (24) 


be the estimates as prescribed by this law, i.e., Bayes estimator using uniform prior. 
Then 


€ rl N; (1) — a) 
i A N l N; l 
al | _ Pë _ Ron mori $ (25) 


2 


If Y, (t) = 0, then 


i id 
( Y, (1) ) N; (r) +] N (t) + 1 l 
GA se 


NO) EE EE o 


N; (1) On a 
This value is used as the summand in Equation (23) whenever Y,(t) = O (zero attrition 
in any year). 


For comparison purposes we again compute an alternate variance which 1s time 


independent. Continuing from Eguation (22), let 


> Y,() 


at 


i= =. (27) 
» N; (1) 
t 
The alternate variance is given by 
l p (1 — p) l 
Kat = Ü RR. — . 28 
(XTB) Zz (XT; (Ð) = > rh (28) 


The problem with cells that have Y, (1) =0 V ¢ (zero attrition for every vear) also 
occurs here, since the variance given by Equation (28) would equal zero. Using the same 


concept as before, we obtain the formula 


GENEE 
(1+ Ý no) EEN 


However in this case, this variance formula is necessary only if all years have zero 


Var(XTB) = (29) 


attrition. 
As before, we return to Eguation (15) to perform the iterative algorithm for 


finding the empirical Baves estimate, XEB,, given bv Equation (14). 
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C. EFRON-MORRIS METHOD 

This method is a modification of the iterative algorithm used to estimate A and 
XBB given by Efron and Morris (Efron and Morris, 1973, pp.127-129). It differs from 
the method given by Equations (14) through (19) in that it allows the variance of the 
prior, A, to change from cell to cell. It also gives greater weight to the cells with low 
variance, and reduces to the James-Stein estimator when the cell variances are constant. 

Onlv one scenario for this method is considered, corresponding to the initial trans- 
formed scale, time dependent variance case. Thus, Equations (4) through (11) are re- 
peated, and we begin from the point where we are entering the iterative algorithm. To 
simplify the following equations, let D, = Var(XTB,) as given by Equation (11). 

We initialize the algorithm with 4,=0 and SP,=0 (previous values of S) for 
r= l1,..., K. Let 





Dir = HH00—— 30 
| A; + D; E 
and 
St 
Vi = K Š (31) 
» e 
= 
Then 
K 
X = ) näTB, (32) 
j=] 
and 
Sar (33) 
Now set i= l and let 
8; — D. f 
SN, = ES WEE V (34) 


2 
ja (A; + Dj) 
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and 


sD, = > ——. (35) 
ji (A; + Dj) 


We then use the Newton-Raphson iteration method to solve 


ET EED SN 
= ee = dA). (36) 
3+(4,+ D) SD, 


First set AP, = 4, for i=1,..., K (previous values of 4). The updated value for A, 


becomes 


AALEN, 


A 


(37) 
If 4,<0, set 4,=0, let i=i+1, and return to Equation (34). If 4,> 0, then test 
|A,— AP,|> e. If true, return to Equation (36). If false, let ¿= i+ 1 and return to 
Eguation (34). 

In either case after incrementing i, if /= K + 1, we exit and test | S, — SP <e Vi. 
If false, we update SP, S, and return to Equation (30) with updated values of A. If 


true, the iterations have converged and we must finalize the estimators, XEM, Let 


l 





d = 3+(4A; +D) )  — (38) 
i i i 2 AN 
and 
= D; 
B = (1-— . 39 
| d Jaz - 
li Bb > 1, setið =E 1;.0r iF B < 0 set 2 SOM MEN 
YER E E (40) 


This equation 1s comparable to Equation (14), which was used to determine the trans- 


formed scale estimates, YEB,, using the previous iterative algorithm. The quantity B, is 
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the amount of shrinkage toward the grand mean, X. The corresponding quantity in 
arl. 


od Va(XTB) | 
To obtain p,, the XEM, must be transformed back to the original scale using XEM, 


in place of XEB, in Equation (20). 


t 


Equation (14) is 


D. VECTOR METHOD 

This method is similar to the previous methods in the sense that it uses the aggre- 
gated cells produced to meet the defined threshold levels. However, prior to the esti- 
mation process, we now partition each aggregated cell bv either service component or 
commissioning source, thus obtaining cells whose elements are vectors. The procedure 
given bv Efron and Morris, modified to compensate for the assumed variance of the time 
averages, provides the framework for this method (Efron and Morris, 1972, pp.341-344). 

The separation by service component or commissioning source reguires us to define 


a third index: the components of the vector. Let 


P = number of service components/commissioning sources. (41) 


K= k= and =a E 


Ni (2) 
Vi) 


inventory of cell í and vector component / in year r 


(42) 


number of attritions in cell í and vector component j in year L. 


Before we had K cells with scalar information, but now we need a K x P matrix where 


p 


P 
> Nyl) = ND) and ) Yy) = KO). (43) 
V=] 


=) 


A requirement for this method is that K> P+ 2, for reasons that will soon become 
obvious. 

We begin by defining X,, (1) as the transformed value for N,,(1) and Y,, (1) as given 
by Equation (4). Continuing in similar manner as the previous transformed scale 
methods, let | 

XT;;(t) = r. ER son d; 


44 
JN; (1) + 5 ` Së 


except when A,, (r) = 0, in which case AT,; (f) does not exist and we reduce 7,, by one. 
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The time averages of the transformed values are then 


XTB; = FLT (e). (45) 
t 


Here we obtain a vector of grand mean values, with each of the P grand means defined 
by 


K 
l 
XBB, = K LATE D (46) 


i=! 


The transformed scale estimate, 6,,, is then found by solving 


jis 


ô = XBB;+ (Ip — (K— P - 2)S”)(XTB;;— XBB) (47) 


where 1, is the identity matrix of order P, and S” is found as discussed below. Reversal 
of the ij index in this and subsequent equations simply means the transpose of the 
K x P matyix to aP < R matin, 


To solve for S7', we begin bv defining 


c T 
where 


The V,, matrix 1s the modification required by our application. (The multiplication in 
Equation (49) is element-wise as opposed to normal matrix multiplication.) The Efron 
and Morris method was developed under the assumption that XTB,- N(6,,,1) , 


whereas We are using 


a, AA _ 
XTB,, cl, = E K, is) (50) 


A 


provided that AT,, (1) has variance of one. Therefore 


l l 
EE mar e 51 


We use the requirement that the P x P matrix resulting from the operations within 
the brackets in Equation (47) must be nonnegative definite to solve for S without 
having to actually compute its inverse. We proceed by doing an eigenanalysis of S, 
which is seen be Equation (48) to be a real symmetric matrix. We form the diagonal 
matrix E, which has the eigenvalues, e, as its diagonal elements, and the matrix I’, which 
has the corresponding eigenvectors as its columns. For any e, <(K-— P — 2), we replace 


it with the value (K — P— 2). The eigenanalysis provides us with the solution to 
i (52) 


Post-multiplying by FF, we obtain 
Secs ET (53) 
since T is ortho-normal and therefore IT T” = I,. We then have 
a= Tem) + PE Tr (54) 


which is easily solved since E"! is found by replacing the diagonal elements of E by their 
reciprocal. This solution for S is then used in Equation (47) to solve for the trans- 
formed scale estimates, ó,,. To obtain the attrition rate estimate for a cell, Di , WE USE 


the inversion formula given by Equation (20) with ö,, in place of ÄEB. 
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IV. CROSS VALIDATION = 


A. GENERAL 

The six estimation methods discussed in Chapter III are evaluated using cross valı- 
dation of the data base. This consists of successively holding out one year's data while 
the other nine years are used to estimate that year’s attrition rates. Three measures of 
effectiveness (MOEs) are used to evaluate the validity of our assumptions and the per- 
formance of the estimation methods. Two of these are original scale MOEs--mean ab- 
solute deviation (MAD) and chi square statistic. The third is a transformed scale 
MOE--mean squared error (MSE). Test cases are chosen as input. The results of the 


cross validation are then discussed. 


B. MEASURES OF EFFECTIVENESS E 
1. Mean Absolute Deviation 

The MAD is probably the most useful MOE to the manpower planner. Our 
version is augmented to display overestimation and underestimation information. Along 
with the MAD we observe the magnitude of our errors in both directions, which is es- 
peciallv important since the cost of overestimating is generallv not the same as the cost 
of underestimating. While it does not provide a specific value or standard to gauge the 
performance of our estimation methods, it does provide very useful insight into tenden- 
cles to consistently underestimate or overestimate. 

For comparison of the estimation methods, we desire a MAD measure that does 
not depend upon cell inventories, vet still displays the overage/underage information. 
For these reasons, we use the attrition rate estimates, p,, as opposed to the estimated 
number of attritions, (p, e N,(1)) (where t = validation year), in our MAD calculations. 
For those estimates obtained in the transformed scale, the XEB, are inverted back to the 
original scale using Eguation (20) prior to calculating this MOE. 


We define the empirical attrition rate for cell i in validation year ¢ as 





Pr N, (0) (59) 


except when A; (r) = 0 (no inventory in cell í for the validation year). In this case we do 


not compute the cell's deviation from the estimated attrition rate since it would 
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artificially create an overage situation. Therefore, we reduce X by one and continue with 
the remaining cells (the reduced value of Kis then used in the following formulas). 


The MAD measures generated for each validation year are 


e = fraction of cells with underage (56) 


where X, is the number of cells which have underage, 


D (p; = py > 
EL oe = fraction of MAD due to underage (57) 
» | pi — Dil 
i 
and 
l A 
MAD = FL Le — Pil - (58) 


We also calculate the average MAD over the validation years. Here we use a 
weighted average, since the number of cells may have been different in some validation 
years, 1.e., a reduced value of K was used in these years. The weighted average takes the 


form 


> K, MAD, 
t 
> K, 
I 


where K, is the (possibly reduced) number of cells used in validation year z. 
2. Chi Square 


The chi square test is used as an indicator of how well the binomial model serves 


Avg MAD = (59) 


as a description of the attrition process. The test statistic is 


(Y, () — p, N, ())° 


60 
iN; (2) D (1 — p) u 


X () = 
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where r is the validation year. As with the MAD calculations, if N. (t) = O we reduce K 
by one and continue. Additionally, if p,= 0 or 1, the denominator equals zero and the 
summand is undefined. The same course of action is used if this occurs--reduce K by 
one and continue. Those estimates obtained in the transformed scale are inverted back 
to the original scale prior to using Equation (60). 

This MOE can be used as a gauge. The chi square statistic given bv Equation 
(60) has expected value K and variance 2K. We are looking for a X' value that is less 


than two standard deviations to the right of the mean, or 


Xi < K+2/2K . (61) 


A weighted average chi square is computed in the same manner as the weighted 
average MAD in Equation (59). However, if the number of cells and thus the degrees 
of freedom, K, are different over the validation years a problem arises in determining the 
degrees of freedom for the weighted average. We solve this dilemma by assuming that 
the weighted average chi square has the original value of K degrees of freedom. 

3. Mean Squared Error 2 

The MSE is used to check the validitv of our theoretical basis. It is the average 
sguared deviation of the estimated rate from the actual rate, both rates on the trans- 
formed scale. The actual rate used is the transformed validation year data. The MSE 


Is defined as 


L(5,1) = = KL D (62) 
where 
d = XEB; 
H = XT;(t) (t= validation year). (63) 


Again, if N, (£) = 0, we reduce K by one and continue. A weighted average MSE is also 
computed similar to Equation (59). 

The MSE also has a standard to gauge our model. Using Equations (5) and (12) 
we can compute a baseline variance for any given validation cell. The MSE for that cell, 
when compared to the baseline value, provides a figure to gauge the value of using 


shrinkage estimators instead of the cell averages, XTB, . There is considerable variability 
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in these ratios, ranging from 20% to 100%, but 80% appears to be a fair median figure. 
For example, for cell variances computed from Eguation (12) running about 0.15, the 
MSE hovers about 0.12. 
4. Vector Method MOEs 

The MOEs discussed above require slight modification before being applied to 
the vector method described in paragraph III.D.. Recall that this method uses K cells 
with service component or commissioning source broken out into a vector of length P. 
An attrition rate estimate, 0,, is obtained for each of the K x P matrix components. 
Thus now we have KP estimates which are compared to the corresponding empirical 
rates for the validation vear. Eguations (55) through (63) are modified by replacing all 
i subscripts with jj, replacing all summations over i by double summations over i and j, 


and replacing all instances of K by the product KP. 


C. TEST CASES 

The selection of test cases takes on great importance since they provide the foun- 
dation for comparison of these methods. It would be impossible to test every permuta- 
tion of input parameters; therefore we seek a representative fraction of these which 
would give an accurate account of the performance of our aggregation and estimation 
methods. Because we are using a different data base from previous theses on estimation 
methods, no attempt to duplicate their test cases was made. 

An approach based upon Latin square experimental design principles was used to 
select 30 test cases for the first five estimation methods. The test cases for the vector 
estimation method are addressed later. In determining the test cases, we randomized 
when possible and intervened to force pairings only when necessary. To begin, we se- 
lected values for the input parameters--7,, K,, grade, YCS and MOS. Service compo- 
nent and commissioning source are ignored for these test cases, i.e., all classifications 
of both are accepted. 

To ensure proper representation from small MOS groups one through ten, one MOS 
from each group was randomly selected: 0302, 1802, 2502, 4002, 3060, 7204, 7545, 7523, 
7557 and 7563. Since YCS and grade are strongly correlated, these parameters were se- 
lected jointly. To ensure each YCS range within the bounded YCS groups was repres- 
ented along with a fair representation of grades, four grade/YCS pairs were selected: 
beget Gs, Capty YES, LtCol/20 YCS and LtCol(failed select)/26 YCS. The two 
threshold parameters were also selected jointly, resulting in ten pairs (74/K,): 30.0/30, 
20.0/30, 20.0/20, 10.0/30, 10.0/20, 10.0/10, 5.0/30, 5.0/20, 5.0/10 and 5.0/5. 
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With these choices in place, it was necessary to combine them to define the actual 
test cases. It was decided to limit the grade/YCS pairs for these cases to 1Lt/4 YCS and 
Capt/7 YCS due to the large values for the first five threshold pairs. With ten MOSs 
specified, we sought ten test cases. Thus, each of the first five threshold pairs was listed 
twice. Each of the four aviation MOSs was randomly assigned to one of the five 
threshold pairs; the six ground MOSs were then randomly assigned to the remaining 
pairs. The two grade/YCS pairs were then randomly assigned within a set of common 
threshold pairs, e.g., for the two cases with 7, /K, of 30.0/30, one was randomly assigned 
1Lt/4 YCS, the other was then assigned Capt/7 YCS. 

All four grade/YCS pairs would be used with the five remaining threshold pairs. 
Thus 20 more test cases were generated, with each of the five threshold pairs listed four 
times. Each MOS was randomly assigned to two distinct threshold pairs, ensuring that 
the large and major MOS groups were evenly spread throughout the pairs. The four 
grade/YCS pairs were assigned in random order to each set of common threshold pairs, 
ensuring that they were evenly spread across large and major MOS groups. The 30 test 
cases are summarized in Table 3. 

The input parameters for six vector test cases were selected from the 30 test cases: 
Nos. 2, 3, 6, 10, H and 20. A small number of vector test cases was initially chosen to 
investigate the possible advantages of the vector method. If this method appeared to 
be favorable, then further testing would be conducted. 

The six test cases contain a cross section of the input parameters. Thev include 
three ground and three aviation MOSs, and use each of the first three grade/YCS pairs 
twice. The grade/YCS pair of LtCol(FS)/26 YCS was not used because of its extremely 
low inventory numbers, which when broken out into a vector would have been of little 
exploratory use. These cases also include six different 7, /K, pairs. 

Each of the vector test cases is used twice: first with service component and then 
with commissioning source as the vector component. All three service components-- 
regular, augmented regular and reserve--were used as vector components. Rather than 
use all 15 commissioning sources (these 15 are listed below Table 4) as vector compo- 
nents (many of them would contain little or no inventory) five commissioning sources 
for the ground test cases and five for the aviation test cases were chosen. These five were 
determined to be the sources which contain the largest percentage of inventory for the 
respective ground or aviation MOSs. The specific commissioning sources selected along 


with the other vector test case input parameters are summarized in Table 4. 
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Table 3. TEST CASES FOR METHODS 1-5 


BEER os ia [ou 
2 | 300 | 30 [7523 | 8:3:2 | 7 | Cap Í | 
A. Lan | 30 | 3060 | 5:2:ı | 7 | Capt _ 
4 | 200 | 30 | 7563 [10:4:2 | a | mt — 
Mk u ste — 1 Ce _ 
6 | 200 | 20 | 7557. l mld | lit |, 
os Pepe a 
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(S: L: M = Small MOS Group : Large MOS Group : Major MOS Group) 
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Table 4. TEST CASES FOR VECTOR METHOD 


Service Comm 


= (5 [5 [55/57 r | mm 
> 300 | so [3 | 7 [cap [ai [12351 
Ta | 2 [z [| á | m [125 | o 
s |o |» (| a | m Im [ass 
Do [100 [20 [2002 e 
iwo | o [ 2 [42 | a |m e iron 
un [00 [io [2502 [a [ca [123 | am > 
oo o pas [20 [uco [um [13710 













Service Component: Commissioning Sources used: 

1 - regular 1 - U.S. Naval Academy 

2 - augmented regular 2 - Platoon Leader Class-Aviation 
3 - reserve 3 - Platoon Leader Class-Ground 


5 - Aviation Officer Candidate 

7 - Officer Candidate Course-Ground 
10 - Enlisted Commissioning Program 
11 - NROTC-Scholarship 
Commissioning Sources not used: 

4 - Platoon Leader Class-Law 

6 - Marine Aviation Cadet 

8 - Officer Candidate Course-Law 

9 - Officer Candidate Course- Women 
12 - NROTC-Ground College 

13 - NROTC-Aviation College 

14 - NESEP 

15 - AIl Other 
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D. RESULTS 

For each test case, we apply the aggregation method to meet the threshold levels and 
then execute the estimation methods. While the ultimate use of these methods is to 
obtain an attrition rate estimate for the original cell, inspection of these estimates would 
be of little value in evaluating and comparing the estimation methods. Thus the output 
from the program takes the form ofthe MOES. 

The inclusion of the entire output from every test case would not only be cumber- 
some but would provide an inadequate means of comparing the methods. Therefore the 
output is summarized in Tables 5 through 7. They contain the output from the first five 
estimation methods only; the output from the vector test cases is presented later. Sam- 
ple output for the six methods is contained in Appendix C. 

The results summaries list the test case number, the cell inventory threshold, Tp, and 
the actual number of cells used, K. The level of expansion required to achieve these 
parameter levels is also listed. For example, test case one had a cell inventory threshold 
of 30.0, and 24 aggregated cells were obtained by expanding the small, large and major 
MOS group by YCSs four and five. The value for Kis often different from the threshold 
number of cells, K,, listed in Table 3. When K is less than K,, maximum expansion oc- 
curred and the threshold was unattainable. When K is greater than K,, the expansion 
was the least amount possible to remain above the thresholds. From these test cases 
we can see that it is difficult to meet the threshold number of cells exactly. 

The results summaries then list the weighted average MOES for each of the five es- 
timation methods. The first row within each test case contains the MAD values, the 
second row the chi square values and the third row the MSE values. The maximum 
desired chi sguare value as given by Eguation (61) is listed in parentheses, e.g., for test 
case one this value is 37.9. This affords easier comparison of the chi sguare values for 
the different methods. The values of MSE for both original scale methods are blank 
because MSE is a transformed scale MOE only. 

Before discussing the results in general, some additional comments about specific 
test cases are necessarv. Test case four could not be executed by the Efron-Morris 
method. This is because one of the aggregated cells contains zero attrition for all ten 
vears of data. Asa result, the iterative algorithm does not converge. 

Test case 12 was not possible because not even one aggregated cell meeting the cell 
inventory threshold was obtained with maximum expansion in major MOS group two. 
This extremely low inventory problem was generally true for all test cases involving the 
LtCol(FS)/26 YCS pair. Test cases 17, 21 and 24 obtained only three aggregated cells 
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with maximum expansion in major MOS group one. Test case 30 was also from major 
MOS group two, and obtained only one aggregated cell meeting the cell inventory 
threshold. As a result, test case 30 was changed to Major(FS)/18 YCS so that results 
for these low thresholds could be obtained. 

As the thresholds became low (test cases 23-30) cells with inventory much larger 
than 7, were being obtained prior to any aggregation. This was especially true for the 
ILt/4 YCS and Capt/7 YCS pairs. To avoid masking the results of low inventory 
thresholds by actually using large inventory cells, the service component/commissioning 
source parameters for test cases 26, 27 and 29 were changed. Rather than accepting all 
classifications of these parameters, only one classification for each was accepted. Thus 
test case 26 was executed with regular/USNA, test case 27 was executed with augmented 
regular/PLC-ground, and test case 29 was executed with regular/NROTC-scholarship. 

We now focus our attention on the results of these test cases with respect to the 
MOEs. The weighted average MAD figures vary little within each test case over the five 
methods. This suggests that the total deviation from the validation year empirical 
attrition rate was the same for all methods. However, this figure does not identify 
whether the deviations were overestimations or underestimations. The fraction of MAD 
from underage (not listed in the results summaries) was studied to gain more insight into 
this important consideration. For each of the 29 test cases (no results for test case 12), 
a weighted average fraction of MAD from underage was computed for each of the five 
estimation methods (weighted by the number of cells just as the weighted averages for 
the MOEs). A weighted average of these 29 values was then computed. This overall 
weighted average indicates the tendency of the method to underestimate or 
overestimate--averages above 0.5 indicate a tendency to underestimate; averages below 
0.5 indicate a tendency to overestimate. The author is unaware of any information 
comparing the relative costs of underage and overage. Hence, as a default, we look for 
values of 0.5, which is a balance between overestimation and underestimation. The av- 
erages calculated were: TS1=0.426; TS2=0.436; OS1=0.481; OS2=0.512; and 
EM = 0.452. Although the MAD figures were generally the same for all methods, these 
averages indicate that the tendencies to overestimate or underestimate may not be the 
same. The original scale methods seem to have achieved more balance than the trans- 
formed scale methods. 

The chi square results were not entirely consistent across test cases nor across 
methods within a test case. These results are discussed first by comparisons between test 


cases; then by comparisons within test cases. 
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Of the first ten test cases, onlv three (Nos. 3, 5 and 8) had weighted average chi 
square values within the acceptable range. These three test cases expanded only into the 
large MOS group, whereas of the seven test cases which were unacceptable, all but one 
(No. 10) expanded into the major MOS group. Of the last 20 test cases, onlv four had 
chi square values outside the acceptable range (Nos. 13, 14, 15, and 18). Of these four, 
two expanded into the large MOS group, and two into the major MOS group. All the 
test cases with unacceptable values had either 1Lt/4 YCS or Capt/7 YCS pairs. They 
were fairly well spread across MOS groups. None of the test cases with lower thresholds 
(Nos. 19-30) had unacceptable chi square values. This suggests that lower thresholds, 
which result in less expansion, achieve more acceptable results with respect to this MOE. 

To investigate this claim further, different combinations of threshold levels for test 
cases seven and nine were used. The results are contained in Table 8. These results re- 
inforce the claim that lower thresholds are in fact better, since in both cases the chi 
square results improved as the thresholds, and therefore the level of expansion, were re- 
duced. When comparing these extra test cases, keep in mind that the chi square values 
for each set of threshold values should be compared to the acceptable range for that 
specific number of cells; comparisons across test cases with different values of K are not 
valid. An important aspect of the argument for lower thresholds is that the thresholds 
must be considered jointly. For example, in test case seven, with a 7,/K pair of 10.0/6, 
the chi square values were nearly acceptable, whereas with 5.0/19 they were clearly un- 
acceptable. Thus we should be most aware of the value 7, x Ko. 

We now turn our attention to comparing the chi sguare values within a test case. 
Several test cases had chi sguare values that varied significantly over the estimation 
methods (Nos. 4, 6 and 22). These three all had ILt/4 X'CS pairs and were large MOS 
group four. In test cases four and six the chi square values for the transformed scale 
methods were not too much larger than the desired maximum; the values for the original 
scale methods were significantly larger than the desired maximum. In test case 22 the 
chi sguare values for the transformed scale methods were acceptable, however the values 
for the original scale methods were again significantly larger than the desired maximum. 

Several other test cases had varving chi sguare values, but to a lesser degree. In test 
case eight, only the transformed scale methods had acceptable chi sguare values, the 
original scale methods exceeded the desired maximum, although not by a significant 
margin. In test case 19, only the OS2 method exceeded the desired maximum. All 
methods for test cases 16, 25, and 27 were within the desired maximum, however the chi 


square values varied to a large degree over the five methods. Thus, using the chi square 
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MOE, it appears that the transformed scale methods were generally the same, and as a 
group outperformed the original scale methods. 

The MSE was used only with the transformed scale methods, and thus no compar- 
ison with original scale methods can be made. The values for this MOE were generally 
egual between methods within a test case, and acceptable overall. 

The results for the vector method test cases are summarized in Table 9. The table 
lists the value K x P (instead of K) because this is the number of estimates obtained and 
compared to the validation year empirical attrition rates with this method. 

Test case two with commissioning source as the vector component had to be modi- 
fied because only seven aggregated cells were obtained. Asa result, K= P+2, and the 
vector method could not be conducted. Therefore, commissioning source three (Platoon 
Leader Class-ground) was deleted as a vector component and the test case run with only 
four commissioning sources (1, 2, 5 and 11). Test case 20 with service component as the 
vector component was infeasible. By starting with a low cell inventory threshold (5.0), 
when the cells were broken out into the vector components their inventorv became ex- 
tremely low. As a result, when validating year five, two of the cells had zero inventory 
for all of the remaining nine vears for service component three (reserve). Thus the value 
for XTB, , given by Equation (45), becomes undefined and the method cannot be 
completed. 

The results of the vector method with respect to the MOEs is similar to the results 
of the previous five estimation methods. This method produced acceptable MAD and 
MSE values, but its chi square values were fairly inconsistent. Test cases two and ten 
had unacceptable chi square values for both vector components; test case six had an 
unacceptable chi square value with service component as the vector compunent. Recall 
that these test cases also had unacceptable values with the vector component collapsed. 

No fair comparison with the first five estimation methods can be made using the 
summarized results. Obviously the MAD and MSE quantities will be larger since we are 
comparing three to five times as many estimates to empirical rates with the vector 
method. Thus a different evaluation technique must be used. 

The vector method 1s designed to take advantage of any correlation between the cells 
when broken out into vector components. To see if this is occurring we must look at 
the matrix S as given by Equation (54). In all of the vector test cases, this matrix was 
essentially diagonal, indicating little correlation between the vector components. In ad- 
dition, all of the eigenvalues, which become the elements of the diagonal matrix E, were 


less than (A — P — 2). Therefore, the eigenvalues were replaced by this quantity and the 
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Table 5. SUMMARY OF RESULTS (CASES 1-10) 
, Expansion 
MAD 0.102 | 0.101 | 0.102 | 0.102 | 0.103 
i X(37.9) | 98.79 | 99.28 | 100.99 | 100.96 | 101.23 
MSE 0.079 | 0.078 0.080 
S:(7) MAD 0.091 | 0.090 | 0.092 | 0.091 | 0.091 
2 30.0 L:(7) X(16.0) | 28.16 | 27.68 | 28.71 | 28.68 | 27.73 
M:(7) MSE 0.077 | 0.076 0.076 
S:(1-3,6-19) | MAD 0.056 | 0.056 | 0.055 | 0.055 | 0.055 
3 20.0 | 31 | L:(6-8) X(46.7) | 37.49 | 37.80 | 39.58 | 42.85 | 38.38 
MSE 0.062 | 0.062 0.062 









30.0 | 24 





MAD 0.029 | 0.029 | 0.026 | 0.024 
X(36.6) | 39.34 | 47.36 | 95.41 | 139.14 





-5.8-19) | MSE 0.035 | 0.037 


20.0 
S:(1-3,6-19) | MAD 0.048 | 0.049 | 0.047 | 0.047 | 0.048 
5 ZU 02 20 P?) 19210) 1 20.75 | 20:83 | 24.34 | 24.20 | 21.31 
MSE 0.048 | 0.049 0.049 
IE ) | MAD 0.029 | 0029 | 0.027 | 0.025: | 0.029 
20.0 | 22 :(1-5. X(35.3) | 40.42 | 50.70 | 90.93 | 127.86 | 65.20 
: MSE 0.037 | 0.039 0.040 


MAD op ME 0 | 0150 1 015011 0165 
X(48.0) | 111.67 | 113.23 | 112.64 | 113.64 | 114.28 
MSE 0.137 0157 0.140 


S:( MAD 0.039 | 0.038 | 0.039 | 0.038 | 0.059 

28 | L:(1-3, X(43.0) | 34.37 | 34.70 | 45.16 | 49.60 | 35.54 
MSE 0.039 | 0.040 0.040 

S:(7) MAD OA ANS A 01179011 0.115 

IS ML) X(24.6) | 38.85 | 37.50 | 39.59 | 39.62 | 40.26 
MII MSE SZ o; 0.152 
S:(4.5) MAD 0.128 | 0122 $ 05129 | 0125 | 0.127 

mL (2) X(41.7) | 64.33 | 65.92 | 65.62 | 66.61 | 62.64 
MSE 0.157 | 0137 0.135 


TS1 - Transformed scale, time dependent variance 
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TS2 - Transformed scale, time independent variance 
OS1 - Original scale, time dependent variance 

OS2 - Original scale, time independent variance 
EM - Efron-Morris method 
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matrix E always had (K — P — 2) as its diagonal elements. Because these results indi- 
cated no worthwhile improvements over the first five methods, no further testing of the 
vector method was conducted. 

Finally, the performance of the heuristic aggregation algorithm listed in Appendix 
A was also evaluated. For each test case, the total inventory of cells below T, was 
summed and this value divided bv 7,. The integer part of this number provides an upper 
bound on the number of aggregated cells that can be obtained by any aggregation 
technique. This upper bound was compared to the actual number of aggregated cells 
produced by the algorithm. The algorithm achieved the maximum in 71.4% (20 of 28) 
of the test cases. It achieved one less than the maximum in 21.4% (6 of 28) of the test 
cases, and two less than the maximum in 7.2% (2 of 28) of the test cases (only 28 test 
cases required aggregation: No. 12 was infeasible; No. 14 all cells were above 7, ). This 


performance is acceptable for our application. 
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Table 6. SUMMARY OF RESULTS (CASES 11-20) 
Expansion 
zs ler Ver Te (e [ost Le [em 
S:(20-25) MAD 0.109 | 0.106 | 0.108 | 0.107 | 0.107 
11 11 | L:(20-25) NO RMS | 19052 29965120557 | 1895 
MSE 0.141 | 0.142 0.139 
S:(26) MAD 
12.1 10:0 14:20) X(0.0) 
M:(26) MSE 
S:(7) MAD VAIM | OTT BEE | 0.117 1.0.TI5 
13 14 | L:(7) X(24.6) | 38.85 | 37.50 | 39.59 | 39.62 | 40.26 
M:(7) MSE 0.130 | 0.127 0.132 
S:(4.5) MAD DEE CLOSE 011571 0.115 
14 | 10.0 L:(4,5) N(19.9) |#57.80 E On aie fl} 57.80 TRES 
MSE 0.102 | 0.103 0.104 
S:(4,5) MAD 0.137 | 0.140 | 0.139 | 0.139 | 0.139 
KA 5:0 32 | Lid, 5 X(48.0) | 69.96 70. 22 | O ET E 
MSE Ie Os 0.155 
S:(20-25) MAD 032 25 OA 101 11102121 | 0.12) 
L:(20-25) X(48.0) A 41 | 40.38 | 44.35 | 47.84 | 42.27 
M:(20-22) MSE 0.146 | 0.146 0.148 
S:(26) MAD 0.166 | 0.166 | 0.168 0.168 
5.0 E:(26) X(7.9) 3.66 3157 3.88 4.02 
M:(26) MSE 0.0.189} 0.188 0202 
al (6. 7) MAD 0.141 | 0.141 | 0.142 | 0.141 | 0.141 
RS 3.0 29 :(6,7) X(44.2) | 79.04 | 78.23 | 78.61 | 81.29 | 78.84 
Y (7) MSE 0.176 | 0.176 0.177 
S:(2,3,6-10) | MAD 0.085 | 0.082 | 0.080 | 0.084 | 0.080 
5.0 24 X(37.9) | 28.95 | 29.28 | 30.91 | 43.36 | 35.44 
MSE 0122 7202123 0.132 
É (20- 25) MAD OMi 0110.1133 í 0:106 | 0.110 
BOS 5.0 AE 2025) de ES) IG 22 1729 [65.92 2537.1 1828 
00) MSE 0.113 | 0.114 0.114 


TS1 - Transformed scale, time dependent variance 
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TS2 - Transformed scale, time independent variance 
OS] - Original scale, time dependent variance 

OS2 - Original scale, time independent variance 

EM - Efron-Morris method 
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SI 7. SUMMARY OF RESULTS (CASES 21-30) 
Expansion 
EHE be [To Te Te 
S:(26) MAD 0.167 | 0.166 | 0.168 | 0.168 | 0.168 
5.0 L:(26) X(7.9) 3.66 3.97 3.88 3.76 4.02 
M:(26) MSE 0.189 | 0.188 0.202 
S:(1-5,8-19) | MAD 0.032 | 0.042 | 0.024 | 0.022 | 0.032 
SS REH 22 | L:(3,4) A(35.3) | 23.56 | 15.03 | 63-10 | 92550] 2289 
MSE 0.059 | 0.041 0.056 
S:(20-25) MAD > Gë E Ü 0.121 x fe 
23 | 50 L:(20) u 5) 9.28 
o 141 0. G 0 137 
MAD 0.166 | 0.166 | 0.168 | 0.168 | 0.168 
24 | 5.0 L Gei X (7.9) 3.66 3:97 3.88 3.78 4.02 
M:(26) MSE 0.189 | 0.188 0.202 
S:(1-6,8-19) | MAD 0.050 | 0.051 | 0.020 | 0.020 | 0.048 
2 A, L:(1-6) X(18.9) | 5.09 4.78 15.35 | 16.54 | 4.85 
MSE 0.044 | 0.037 0.039 
S:(1-3,6-19) | MAD 0.138 | 0.138 | 0.138 | 0,135 (“089 
26 | 5.0 L:(6,7) X(20.4) | 13.61 | 13.65 | 14.09 | 16.10 | 14.28 
MSE 0.184 | 0.186 0.189 
S:(6-8) MAD 0.070 | 0.069 í 0.055 | 0.050 | 0.068 
2 510 X(12.9) | 4.64 4.55 6.88 9.41 4.60 
MSE 0.075 | 0.074 0.074 
S:(20,21) MAD 0.114 | 0.116 | 0.118 } 0.113 | 0.116 
282150 7 X(14.5) | 7.59 7.67 7.99 8.74 8.37 
MSE 0.140 | 0.143 0.147 
S:(4,5) MAD 0.170 | 0.169 | 0.168 | 0.167 | 0.169 
29 | 50 5 X(11.3) | 6.36 6.23 7.64 8.25 ! 
MSE 0.211 | 0.210 
S:(1-6,8-19) | MAD 0.112 | 0.110 | 0.110 
30 | 5.0 6 | L:(17-19) AZ ONIS 6.40 7.31 
MSE 01156 | 0,136 


TS1 - Transformed scale, time dependent variance 













TS2 - Transformed scale, time independent variance 
OS1 - Original scale, time dependent variance 

OS2 - Original scale, time independent variance 

EM - Efron-Morris method 
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Table 8. SUMMARY OF RESULTS (CASES 7 AND 9 EXPANDED) 
z xpansion 

rf [| ge [o [o [o [om [os 
S:(4,5) 

7-1 32 | L:(4,5) X(48.0) | 111.67 | 113.23 | 112.64 | 113.64 | 114.28 
M:(4) 
S:(4,5) 

7-2 17 LA) X(28.7) | 45.31 | 45.99 | 45.90 | 46.58 | 43.68 
S:(4,5) o 

7-3 X(12.9) | 13.08 | 12.98 | 13.41 | 13.53 | 13.03 
S:(4,5) 

7-4 | 5.0 L:(4) X(31.3) | 49.28 | 49.53 | 49.23 | 50.84 
S:(4,5) 

ES | 5.0 7 X(18.9) | 14.42 | 14.27 | 14.82 | 14.93 | 14.18 
S:(4) 

7-6 | 5.0 X(9.7) 9.20 8.88 922 9:20 
LH?) X(24.6) | 38.85 | 37.50 | 39.59 | 39.62 | 40.26 
M:(7) 
S:(7) 

9-2 | 10.0 À(9.7) 1055 TOTS 710.58 2210.57 OM 
S:(7) 

9-3 | 5.0 18:17) X(28.7) | 44.95 | 42.87 | 44.40 | 46.03 | 46.15 
B 


TS1 - Transformed scale, time dependent variance 

















TS2 - Transformed scale, time independent variance 
OS1 - Original scale, time dependent variance 

OS2 - Original scale, time independent variance 
EM - Efron-Morris method 
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Table 9. SUMMARY OF RESULTS (VECTOR METHOD) 
Vector Vector 
Component Method 
MAD 0.149 
2 0.0 SC X(37.9) 61.03 
MSE 0.229 
MAD 0,222 
2 30.0 28 CS X(43.0) 115.34 
MSE 0.398 
MAD 0.158 
3 20.0 95 SC X(120.3) 110.96 
0.219 
02157 
3 20.0 155 CS X(190.2) 132.537 
0.201 


Ä 0.062 
X(89.0) 109.19 
0.074 


0.083 
20.0 X(128.3) 58.58 
0.081 


0.212 

SC X(106.5) 251.63 

0.368 

MAD 0.321 

10 10.0 120 CS X(151.0) 707.85 
0.625 


0.193 
X(60.3) 44.86 
0.214 


TID 
10.0 CS X(81.9) 53.55 
0.225 

















































SC - service component 


CS - commissioning source 


V. CONCLUSIONS AND RECOMMENDATIONS 
A. CONCLUSIONS 


The results indicate that the desired stability in estimating attrition rates for low in- 
ventory cells has been achieved with the aggregation and estimation methods presented 
in this study. The use of “shrinkage” methods applied to well selected groups of cells 
allows for the achievement of quality estimates of attrition in the face of low inventory 
numbers for the individual cells. 

None of the six estimation methods stood out as a clear favorite. The vector 
method did not provide any additional benefits using service component or commis- 
sioning source as vector components. Since it is a more complicated method and has 
the potential to become unsolvable with zero inventory vector components, it appears 
to be the least favorite. Perhaps more success would be obtained with alternative clas- 
sifications for the vector component. 

The Efron-Morris method also involves more computational effort than the first 
four empirical Bayes methods. Its performance was very much similar to the trans- 
formed scale, time dependent variance method since the only difference between them 
is the iterative algorithm used to determine the amount of shrinkage. The Efron-Morris 
method has the potential to become unsolvable when a cell has zero attrition for every 
year--a distinct possibility when dealing with low inventory cells. This suggests that it 
is the least favorite of the first five methods. i 

Of the remaining four methods, there seems to be only small difference between the 
time dependent variance and the time independent variance methods on the same scale. 
In test cases where the chi sguare values were marginal or unacceptable, the time de- 
pendent variance methods were usually better. In these same test cases, the transformed 
scale methods performed better than the original scale methods. Therefore, if one 
method was to be singled out as best, it would be the first method: transformed scale, 
time dependent variance. 

The tendency to overestimate or underestimate as shown by the weighted average 
fraction of MAD from underage may also be a consideration when selecting a method. 
An analysis of this type must weigh the costs of overestimating versus the costs of 


underestimating, which generally are not the same. This type of an analysis 1s beyond 
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the scope of this study. In addition, further testing of the methods would be required 
to gain a more accurate estimate of this tendency. 

The threshold levels also seem to strongly influence the performance of the esti- 
mation methods. It appears that expansion past the large MOS group begins to detract 
from homogeneous attrition behavior. While further study would be required to identify 
optimal threshold levels, it 1s apparent that both thresholds should not exceed 20.0, and 
the value of 7, x K, should not exceed 100. 

A method for dealing with cells whose inventory 1s much greater than 7, must be 
developed. In some test cases, cells with inventory three or more times as large as 7, 
were obtained and used in the estimation process. This did not seem to affect the results, 
as they were present in almost all test cases. These cells could be disaggregated into 
multiple cells with inventory closer to the threshold, although the effect of this has not 


been determined. 


B. RECOMMENDATIONS 

The proposed aggregation method should be implemented as a method of identifying 
additional cells to be used in the attrition rate estimation process. This method provides 
greater homogeneity of attrition behavior among cells over the current method. 

The empirical Baves estimation methods developed in this study are recommended 
for use in estimating the attrition rates for low personnel inventorv cells. 

Overall, the empirical Bayes estimation methods when combined with the proposed 
aggregation method have achieved the stability in attrition rate estimation that is re- 


quired to provide a foundation for manpower planning. 
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APPENDIX A. AGGREGATION ALGORITHMS 


A. HEURISTIC ALGORITHM 


l. 


The heuristic algorithm for aggregating cells is as follows: 


Given a set of cells, S, and the (time average) inventory of each cell, INV,, partition 
S into two subsets as follows: 

S Gs ces: INK > Ty) 

Sy CES IW < Th) 


Put the cells in S, into the set of aggregated cells, K. 


Order the cells in $, according to size of their inventorv: 
INV = INV, = + + + < INV n= |S,| 


. Start with c,, the cell in S, with the largest inventory. Find the smallest cell in S,, 


c', that when united with c, the resulting total inventory will meet or exceed 7, . 
Combine its data with c,, put c, into K, and remove c' from S, (the modified set S, 
will now be referred to as S; ). Repeat the procedure with c,_,, and so forth. 


If no cell in 5; when combined with the current largest cell, c,_, , exceeds 7,, use 
the next largest cell, c, .,, and remove c, , , from 5, . This will create an aggregated 
cell that is still below threshold. Return to the procedure in Step 4 of trving to find 
Geer ro such cell is contained in S: , use C-z, and so forth. 


Continue this procedure until the sum of all the cells remaining in $; is less than 
T,. These cells are sequentially added to the aggregated cells in K in Step 7. 


Add the largest cell in S; to the smallest cell in K, and update its average inventory. 
Add the next largest cell in S; to the current smallest cell in K, and update the in- 
ventory. Continue until all cells in S; have been used. 


We now have | K| aggregated cells which exceed the threshold, 7,, to use in the 


attrition rate estimation procedure. 


B. 


INTEGER LINEAR PROGRAM 


The formulation as an integer linear program 1s as follows: 


Index Use 
6 cell (before aggregation) 
a aggregated cell 
Given Data 
INP average inventory of cell c 
T threshold cell inventory 
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Binarv Variables 


X l indicates use cell c in aggregated cell a 

2; l indicates use cell a for aggregation 
Formulation 

MAX `> Z, 


subject to | 
E 1 Ve (each cell used at most once) 
> INV, + X,, 2 h.e Z, Va (aggregated cell must have size > 7, ) 
Nees Z, € (0,1) 
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APPENDIX B. COMPUTER PROGRAMS di 
A. GENERAL 


A computer program written in FORTRAN is used to conduct the cross validation 
using the methods developed in this thesis. Although the program consists of 33 sub- 
routines, 6 function subroutines, and over 2000 lines of code, it can be easily summarized 
by breaking it into the two areas of the thesis: cell aggregation and estimation methods. 

The main program and aggregation subroutines (listed in paragraph B) read the in- 
put parameters and execute the expansion and aggregation methods discussed in Chap- 
ter II. An existing program written bv Luis Uribe, an independent contractor under the 
direction of Professor Read, underwent extensive modification to fulfill these require- 
ments. The input parameters--7,, K, , MOS, YCS, grade, service component(s), and 
commissioning source(s)--are read by Subroutine GETPAR either in the interactive 
mode via the terminal or by using MC87 EXEC (listed in paragraph E). Uribe uses an 
innovative method to estimate the amount of expansion required to meet the threshold 
parameters. This approach precludes the requirement to read the data base after each 
step in the expansion process which would be extremely computer time intensive. In- 
ventory information is extracted from the data base and stored in a separate data file for 
each pay grade (a sample data file and the program used to create it are listed in para- 
graph F). The data file is accessed via the user’s A-disk which is significantly faster than 
accessing the data base through MVS. Subroutine READET reads the appropriate data 
file for the specified grade and constructs a table of cells for those records that are in the 
same major MOS group as the user defined MOS, and meet the service component and 
commissioning source parameters. All YCSs are accepted, since the extent of expansion 
is not yet determined. Function NCEVAL screens this table using the current level of 
expansion and estimates the number of aggregated cells with average inventory greater 
than or equal to 7, that will be obtained. If this number is less than A,, Subroutine 
EXPAND begins the expansion stages as described in paragraph H.C.. After each in- 
crement of expansion, NCEVAL screens the table and estimates the number of aggre- 
gated cells that will be obtained. This loop through EXPAND and NCEVAL continues 
until the estimated number of aggregated cells meets the threshold, K, The estimated 


number of cells and the level of expansion are then displaved on the terminal screen. 
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The user may elect to go forth and read the data base to determine the actual number 
of aggregated cells obtained, or may elect to change the level of expansion. 

The level of expansion 1s changed through the variable AGGPCT. This variable 
estimates the effectiveness of the heuristic aggregation method listed in Appendix A. 
To estimate the number of aggregated cells that will be obtained, NCEVAL compares 
the cells which meet the expansion criteria to the minimum inventory threshold, 7. 
Those that are greater than T, will obviously produce one aggregated cell. The inventory 
of those that are less than 7, is summed. The estimated number of cells is then the total 
of the number of cells greater than 7, and AGGPCT times the sum of the cell inventory 
below T, divided by T, The variable AGGPCT is initially set at 0.9, but can be inter- 
actively changed via the terminal. By increasing the value of AGGPCT we can decrease 
the level of expansion; by decreasing the value of AGGPCT we can increase the level 
of expansion. 

Once we decide to go forth and read the actual data base, Subroutine READER 
extracts records meeting the expansion criteria developed using EXPAND and 
NCEVAL and pools them into cells. Subroutine AGGREG aggregates these cells to 
meet the average inventory threshold, 7, . The actual number of aggregated cells ob- 
tained is then compared to the threshold number of cells, K,. Again, the user has the 
option of changing the level of expansion to obtain more or fewer cells, or continuing 
on to the estimation process. 

The first five estimation methods are called by SUBROUTINE MC87BZ (listed in 
paragraph C). The estimation methods are contained in separate subroutines: EBTSI, 
EBTS2, EBOSI, EBOS2 and EMIS (EB-empirical Bayes; EM-Efron-Morris; 
TS-transformed scale; OS-original scale; 1-time dependent variance; 2-time independent 
variance). The iterations required by the first four methods are conducted in Subroutine 
EBITER; the Efron-Morris iterations are conducted in Subroutine EMITER. The 
MOEs are then computed by Subroutines MSE and OSMOE. 

If the vector method is to be used, Subroutine BK DOWN then breaks the cells out 
by their vector components (a vector of length three for service component; a vector of 
length five for commissioning source). The vector estimation method is contained in 
Subroutine MC87V (listed in paragraph D). Since all of its computations are unique, 
this subroutine is self-contained with the exception of the transformation formula, which 


is contained in Function FTTV. 
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W 


C --- PROGRAM TO CONDUCT AGGREGATION AND ESTIMATION METHODS 
€ 
C --- PARAMETER MXY MUST BE UPDATED TO REFLECT EXACT NO. YEARS OF DATA 
C --- PARAMETER MXP IS MAX LENGTH OF 3RD DIMENSION P-VECTOR 
C --- PARAMETER MXK IS MAX NUMBER OF AGGREGATED CELLS (MAX NO) 
PARAMETER (MXX=600, MXY=10, MXP=6, MXK=50) 
PARAMETER (NMS=81, NG=14, NLG=6, NMG=4, NYB=4, NYE=18, NYEG=4) 
E 
INTEGER ST1,ST2,LYR 
INTEGER SYCS(31), NYCS 
INTEGER SYCSG(31),SYCSL(31),SYCSM(31),NYCSG,NYCSL,NYCSM 
INTEGER SMOS(30), NMOS 
INTEGER SVCMP(5), NSC 
INTEGER SCSRC( 16),NCSR 
INTEGER SGRD 
INTEGER*2 MOSGR(2,NMS), YCSB(NYE,NYB,NYEG), VYC(NYE) 
INTEGER*2 LGRP(NG), MGRP(NLG) 
REAL INV(MXX,MXY) , Y(MXX ,MXY), SINV(MXX,MXY) , SY(MXX,MXY) 
INTEGER DATA(MXY) 
C --- ARRAYS FOR MC87BZ 
REAL XTB(MXX) , VXTB(MXX) , XEB(MXX) , A(MXX) 
C --- ARRAYS FOR MC87V 
REAL XTBJI(MXP,MXK), DELTA(MXP,MXK), X(MXP,MXK) 
REAL XVYR(MXP,MXK), VYRINV(MXP,MXK), VYRY(MXP,MXK) e 
REAL BSTAR(MXP,MXP), S(MXP,MXP), GAMMA(MXP,MXP) 
REAL XBBJ(MXP), EVAL(MXP) 
E 
INTEGER*2 PTRTBL(MXX, 2), INDX(MXX), MKG(MXX), RETTBL(MXX,3) 
INTEGER*2 PTBL(MXX, 3), BKTBL(MXX,3) 
REAL AVINV(MXX), RETINV(MXX) 
DATA MKG/MXX*O/ 
C --- ASSIGN MOS TO SMALL, LARGE AND MAJOR MOS GROUP 
BR MOSGR/013,1, 020,2, 027 2, 038 2, 039.2, 
* 005,3, 007,3, 049,3, 052,3, 
* 074,4, 079,4, 085,4, 101,4, 
* 016,5, 060,5, 064,5, 076,5, 111,5, 116,5, 
* 132,6, 134,6, 135,6, 139,6, 
“AS 7. 147 7, 150,7, 153,7, 154,7, 155,7, 12057, 
* 149,8, 151,8, 
* 160,9, 161,9, 164,9, 166,9, 167,9, 168,9, 178,9, 
“ET io ie 10, 175,10, 176,10, 177,10, 179,10, 128310, 
* 145,10, 165,10, 
“EN ii 006,11, 012,11, 015,11, 019,11, 026,11, 037777. 
“Es  05111/0059/:11,0070,11,8 075,11, 078,11, 084,11, 
O ROO 138,11, 217,11, 
“ET 2, 187,12, 188,12, 189,12, 
* 142,13, 146,13, 148,13, 152,13, 156,13, 163,13, 169,13, 
* 088,14 
DATA LGRP/1,1,4*2,3,3,4,4,5,5,5,6/ 
DATA MGRP/1,1,2,2,3,4/ 
C --- CREATE YCS EXPANSION BOUNDS 


Sp 
ale 
dh 


7 17*0, 


1. ,2 os A SE 10, no 12, ioe Ss 056, 17,18, 19, 150, 


20, 21, 22, 23, 24, 25, 12: 0, 
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MAIN PROGRAM AND AGGREGATION SUBROUTINES 


DARS ol. 2 o 0. 819105 12 13 MERE 16,17,18,19, 


26, 1750, 


MC800010 
MC800020 
MC8 00030 
MC800040 
MC8 00050 
MC8 00060 
MC800070 
MC800080 
MC8 00090 
MC800100 
MC800110 
MC800120 
MC800130 
MC8 00140 
MC800150 
MC800160 
MC800170 
MC800180 
MC800190 
MC800200 
MC800210 
MC800220 
MC800230 
MC800240 
MC800250 
MC800260 
MC800270 
MC800280 
MC800290 
MC800300 
MC800310 
MC800320 
MC800330 
MC800340 
MC800350 
MC800360 
MC800370 
MC800380 
MC800390 
MC800400 
MC800410 
MC800420 
MC800430 
MC800440 
MC8 00450 
MC8 00460 
MC800470 
MC8 00480 
MC800490 
MC800500 
MC800510 
MC800520 
MC800530 
MC800540 


* 6,7,16*0, 00 20,21,22,23, 24,25, 1240 DEE ON MC800 
* 1,2,3,4,56, B,9,10, LI 02 19, 4,15, DO RAIN MC800 
* 7.1770, 20,21 22.23, 25 25 12%0, 26 HON MC800 
* 1.2.3, 6.7.B,9 10,11, 12,18 14,15. 16 EE on MC800 
* 4,5,16%0) 00 21,22 28 2425 120) 2 6 N H MC800 
C --- INITIALIZE INVENTORY AND ATTRITION ARRAYS MC800 
DO 1 I=1,MXX MC800 
DO ET MC800 
SINV(1,J)=0 MC800 
SY(1,J)=0 MC800 
INV(1,J)=0 MC800 
Y(1,J)=0 MC800 
2 CONTINUE MC800 
1 CONTINUE MC800i 
C --- DEFINE FILE FOR OUTPUT MC8001 
CALL EXCMS('FILEDEF 11 DISK MC87 OUTPUT A') MC800) 
C --- FIRST/LAST YEAR OF DATA ON TAPE. UPDATE WHEN NECESSARY MC800) 
STi=77 MC800 
NYR=MXY MC800 
LYR=ST1+NYR-1 MC800: 
C --- INITIAL VALUE FOR AGGREGATION ESTIMATION PERCENTAGE MC800; 
AGGPCT=0. 9 MC800; 
je ep MC800; 
C --- GET INPUT PARAMETERS MC800; 
CALL GETPAR(AIMIN,NO ,NMOS , SMOS ,NYCS,SYCS,SGRD, MC8007 
* NSC,SVCMP, NCSR,SCSRC, IGR,MOSGR,NMS, ISFLAG) MC800£ 
C --- MAJOR GROUP IS MG, LARGE GROUP LG, GROUP IGR, YCS BLOCK IY MC800£ 
LG=LGRP( IGR) MC800£ 
MG=MGRP( LG) MC800£ 
WRITE(6,*) ' ' MC800£ 
WRITE(6,*) '---- GR,LG,MG=' ,IGR,LG,MG MC800€ 
WRITE(6,*) ' ' MC800€ 
C --- READ EVALUATION TABLE. SELECT ONLY RECS PASSING SELECTION CRITERIAMC8008 


Saas READET(RETTBL,RETINV,MXX,NRET, SGRD, NSC,SVCMP, NCSR,SCSRC, 


MG,LGRP,MGRP, 'MOSGR „NMS ) 

5 K 0 

IGX=IGR 

LGX=0 

MGX=0 

NYCSG=1 

SYCSG(1)=SYCS(1) 

NYCSL=1 

SYCSL(1)=SYCS(1) 

NYCSM=1 

SYCSM(1)=SYCS(1) 

NCTOT=0 


NC=NCEVAL(AIMIN, IGX, LGX, MGX, NYCSG, SYCSG, RETTBL, RETINV,NRET,MXX, 


* LGRP,MGRP,NMS,AGGPCT, IGR,LG) 


C --- DO WHILE NCTOT<NO & RC=0 (EXPAND AS LONG AS NO NOT MET) 
10 IF(NC .GE. NO) THEN 
WRITE(6,*) '$GG EVAL NC,SYCSG=' ,NC, (SYCSG(II),II=1,NYCSG) 
GO TO 60 
ENDIF 


IF(NYCSG. EQ. 1) THEN 
CALL GETVYC( SYCSC1) , LG, YCSB,NYE,NYB,NYEG,VYC) 
WRITE(6,*) '=== VYC=' ,(VYC(T), I=1 ‚NYE) 
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MC8008 


MC8009? 
MC8009! 
MC8009) 


MC8010% 
MC80101 
MG8011( 


ENDIF 

CALL EXPAND(NYCSG,SYCSG,VYC,NYE,IGX,LGX,MGX,LG,MG,RC) 

IF(IGX .EQ. 0) GO TO 20 

NC=NCEVAL(AIMIN, IGX, LGX,MGX, NYCSG, SYCSG, RETTBL, RETINV, NRET, MXX, 
* LGRP,MGRP,NMS,AGGPCT, IGR,LG) 

COT TOO 

20 NCTOT=NC 
WRITE(6,*) '$$G EVAL NC,SYCSG=' ,NCTOT,(SYCSG(II),II=1,NYCSG) 


--- EXPAND TO LARGE MOS GROUP 
WRITE(6,*) ' ' 
WRITE(6,*) '=== EXPANDING BY LARGE GROUP: ' ,LGX 
NC=NCEVAL(AIMIN,IGX, LGX,MGX ,NYCSL, SYCSL,RETTBL,RETINV,NRET,MXX, 
* LGRP,MGRP,NMS ,AGGPCT,IGR,LG) 
30 IF((NCTOT+NC) .GE. NO) THEN 


MC801110 
MC801120 
MC801130 
MC801140 
MC801150 
MC801160 
MC801170 
MC801180 
MC801190 
MC801200 
MC801210 
MC801220 
MC801230 
MC801240 
MC801250 


WRITE(6,*)'$LL EVAL NC, SYCSL=",(NCTOT+NC), (SYCSL(II), II=1,NYCSL)MC801260 


GO TO 60 
ENDIF 
IF(NYCSL. EQ. 1) CALL GETVYC(SYCS(1),LG,YCSB,NYE,NYB,NYEG,VYC) 
CALL EXPAND(NYCSL, SYCSL, VYC,NYE, IGX,LGX,MGX,LG,MG,RC) 
M Gr PO o) GO TO 40 
NC=NCEVAL( AIMIN, IGX,LGX,MGX,NYCSL,SYCSL,RETTBL,RETINV,NRET,MXX, 
* LGRP,MGRP,NMS,AGGPCT, IGR,LG) 
CORTO 30. 
40 NCTOT=NCTOT+NC 
WRITE(6,*) "SSL EVAL NC,SYCSL=' ‚NCTOT,(SYCSL(II),II=1,NYCSL) 


--- EXPAND TO MAJOR MOS GROUP 
WRITE(6,*) ' ' 
WRITE(6,*) '=== EXPANDING BY MAJOR GROUP: ' ,MGX 
NC=NCEVAL(AIMIN, IGX, LGX, MGX, NYCSM, SYCSM, RETTBL, RETINV, NRET,MXX, 
* LGRP,MGRP,NMS,AGGPCT, IGR,LG) 
50 IF((NCTOT+NC) .GE. NO .OR. RC .NE. 0) THEN 


MC801270 
MC801280 
MC801290 
MC801300 
MC801310 
MC801320 
MC801330 
MC801340 
MC801350 
MC801360 
MC801370 
MC801380 
MC801390 
MC801400 
MC801410 
MC801420 
MC801430 


WRITE(6,**)'SMM EVAL NC,SYCSM=" ,(NC+NCTOT) , (SYCSM( II), II=1,NYCSM)MC801440 


GO TO 60 
ENDIF 
IF(CNYCSM. EQ’ 1) CALL GETVYC(SYCS(1),LG,YCSB,NYE,NYB,NYEG, VYC) 
CALL EXPAND(NYCSM, SYCSM, VYC ,NYE ,IGX,LGX ,MGX,LG,MG,RC) 
NC=NCEVAL( AIMIN, IGX, LGX,MGX,NYCSM ,SYCSM, RETTBL, RETINV,NRET ,MXX, 
* __ LGRP ,MGRP ,NMS , AGGPCT , IGR , LG) 
GO TO 50 


--- EXPANSION FINISHED 
60 ABRES. NES 0) THEN 


MC801450 
MC801460 
MC801470 
MC801480 
MC801490 
MC801500 
MC801510 
MC801520 
MC801530 
MC801540 


WRITE(5,*)'**** REOUIRED NO MAY NOT BE MET: NO,NC=',NO,(NC+NCTOT)MC801550 


ENDIF 
--- ALLOW USER TO CHANGE EXPANSION LEVEL 
WRITE(5,%) 'ESTIMATED NUMBER OF CELLS =' ,NC+NCTOT 
70 WRITE(5,*) 
WRITE(5,*) 'ENTER 1 TO CALL READER, 0 TO CHANGE EXPANSION' 
READ(5,*) NPICK1 
IF(NPICK1 .EQ. 1) THEN 
GO TO 80 
ELSE 
WRITE(5,*) 'AGGPCT IS CURRENTLY =', AGGPCT 
WRITE(5,*) 'ENTER NEW VALUE FOR AGGPCT' 
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MC801560 
MC801570 
MC801580 
MC801590 
MC801600 
MC801610 
MC801620 
MC801630 
MC801640 
MC801650 
MC801660 


READ(5,*) AGGPCT MC801: 








GO TO 5 MC801: 

ENDIF MC801: 

80 WRITE(5,%) 'CALLING READER' MC801! 

Ë MC801' 
C --- USER ELECTS TO READ THE DATA BASE - DETERMINE MOS EXPANSION LEVEL MC801. 
CALL GETMOS(SMOS , NMOS ,MGX, LGX,MG,LG, IGR,MOSGR,LGRP,MGRP, MC801 

* NMS,NG,NLG) MC801 

C --- READ THE DATA BASE AND CREATE THE CELLS MC801.| 


CALL READER( DATA, INV,Y,MXX,NMOS ,NYCSG,NYCSL,NYCSM,NSC ,NCSR,NYR, MC801 | 
* SMOS,SYCSG, SYCSL, SYCSM, SGRD, SVCMP, SCSRC,NRC, PTRTBL, LGX, MGX, IGR, MC80T32 
* LG,MGRP,LGRP,MOSGR,NMS,NG,NLG, ICYCLE,NPT,PTBL, ISFLAG,SINV,SY) MC801. 


C --- PERFORM CELL AGGREGATION TO MEET INVENTORY THRESHOLD MC801; 
CALL AGGREG( INV, Y ,MXX,NYR,SMOS,SYCSG, MC801£ 

* NRC, NRCOLD,PTRTBL, INDX, AVINV, AIMIN,MKG) MC801£ 

C --- ALLOW USER TO CHANGE EXPANSION LEVEL MC801€ 
WRITE(5,*) "NUMBER OF CELLS =' ,NRC MC801€ 

90 WRITE(5,*) MC801E 
WRITE(5,*) 'ENTER 1 TO CONTINUE, 0 TO CHANGE EXPANSION' MC801£ 
READ(5,*) NPICK2 MC801£ 
IF(NPICK2 .EQ. 1) THEN MC801£ 

GO TO 100 MC801£ 

ELSE MC801£ 
WRITE(5,*) 'AGGPCT IS CURRENTLY =", AGGPCT MC801S' 
WRITE(5,*) "ENTER NEW VALUE FOR AGGPCT' MC801$ 
READ(5,*) AGGPCT MC801$ 
ICYCLE=ICYCLE+1 MC8019 

GO TO 5 MC8019) 

ENDIF MC8019' 

C MC8019 
C --- USER ELECTS TO CONDUCT ESTIMATION MC8019: 
100 CONTINUE MC8019' 
WRITE(11,201)'EXPANSION INFORMATION: ' MC8019- 
WRITE(11,203)'ACTUAL NO. OF CELLS USED= ' ,NRC MC8020. 
WRITE(11,202)'MOS GROUP f4' ,IGR,' YCS''S USED=', MC8020. 

* (SYCSG(I),I=1,NYCSG) MC8020: 
IF(LGX .GT. 0) THEN MC8020: 
WRITE(11,204)' LARGE MOS GROUP #' ,LG,' YCS''S USED=', MC8020: 

* (SYCSL(1),1=1,NYCSL) MC8020 

ELSE IF(MGX .GT. 0) THEN MC8020) 
WRITE(11,204)'LARGE MOS GROUP #',LG,' YCS''S USED=', MC8020) 

* (SYCSL(I),I=1,NYCSL) MC8020) 
WRITE(11,204)'MAJOR MOS GROUP #',MG,' YCS''S USED=', MC8020) 

* (SYCSM(I),I=1,NYCSM) MC80211 
ENDIF MC8021 | 
C --- PERFORM ALL BUT VECTOR ESTIMATION METHODS IN MC87BZ _ MC8021:| 
CALL MC87BZ(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) MC8021:) 
C --- VECTOR METHOD--BREAK CELLS INTO VECTOR, CONDUCT ESTIMATION MC802 1º 
IF(ISFLAG .GT. 0) THEN MC8021:; 
CALL BKDOWN(PTBL,NPT,PTRTBL,NRCOLD,INDX,MKG,MXX,MXY, MC802 14 
* SINV,SY,INV,Y,BKTBL,NBK) MC8021; 
CALL MC87V(INV,Y,MXX,NYR,NRC,XTBJI,DELTA,X,XVYR,VYRINV,VYRY, MC8021€ 
* BSTAR,S,GAMMA, XBBJ ,EVAL,MXP,MXK, BKTBL,NBK,NSC ,NCSR,ISFLAG) MC8021¢ 
ENDIF MC8022( 

C MC80221 
201 FORMAT(/1X,A) MC8022: 
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202 FORMAT(1X,A,12,A/1X,18(13)) MC802230 
203 FORMAT(1X,A,12) MC802240 
204 FORMAT(1X,A,11,A/1X, 18(13)) MC802250 

END MC802260 

C MC802270 

WHEW A de eos rd dee He Mee He Wee se He He He Fede Were Me Fe se He je Tele ge fe te de He He HICH HH e fe fe fe e e fe eck de de decke de decke dee dee eck MT SD 2 280 

C MC802290 

SUBROUTINE EXPAND(NYCSX,SYCSX,VYC,NYE, IGX,LGX,MGX,LG,MG,RC) MC802300 

C --- EXPAND YCS IF FEASIBLE, ELSE EXPAND MOS TO LG/MG MC802310 

INTEGER SYCSX(31), NYCSX MC802320 
INTEGER*2 VYC(NYE) MC802330 

C --- FIND POSITION OF ORIGINALLY REQUESTED SYCS(1) MC802340 
1Y=0 MC802350 

DO 10 I=1,NYE MC802360 
IF(SYCSX(1) .EQ. VXC(I)) IY=I MC802370 

10 CONTINUE MC802380 
IF(IY.EQ.0) GO TO 30 MC802390 

C --- FIND NEAREST NON-ZERO YCS TO USE FOR EXPANSION MC802400 
DO 20 I=1,NYE MC802410 
J=IY-I MC802420 

IF(J. GE. 1) THEN MC802430 
IF(VYC( J). GT. 0) GO TO 50 MC802440 

ENDIF MC802450 

J=IY+I MC802460 

IF( J. LE. NYE) THEN MC802470 
ME) GT. 07460910 50 MC802480 

ENDIF MC802490 

20 CONTINUE MC802500 
30 CONTINUE MC802510 

C --- NO MORE YCS EXPANSION POSSIBLE. SEE IF MOS EXP. FEASIBLE MC802520 
IF( IGX. GT. 0) THEN MC802530 

GNS EXPAND FROM GROUPS TO LARGE GROUP LGX MC802540 
IGX=0 MC802550 

LGX=LG MC802560 

ELSE IF(LGX.GT.0) THEN MC802570 

C --- EXPAND FROM LARGE GROUP LGX TO MAJOR GROUP MGX MC802580 
LGX=0 MC802590 

MGX=MG MC802600 

ELSE MC802610 

RC=1 MC802620 

ENDIF MC802630 
RETURN MC802640 

C MC802650 
C --- EXPAND WITH YCS IN POSITION J & CLEAR VYC(J) MC802660 
50 CONTINUE MC802670 
NYCSX=NYCSX+1 MC802680 
SYCSX(NYCSX)=VYC( J) MC802690 
VYC(J)=0 MC802700 

END MC802710 

C MC802720 

ee ET Here a EK IERE MORO27 30 

C MC802740 

FUNCTION NCEVAL(AIMIN, IGX, LGX, MGX, NYCSX, SYCSX, RETTBL, RETINV, MC802750 
NRET,MXX, LGRP,MGRP ,NMS , AGGPCT,IGR,LG) MC802760 

C --- COMPUTE ESTIMATED NO. CELLS TO BE OBTAINED WITH CURRENT EXPANSION MC802770 
INTEGER SYCSX(31),NYCSX MC802780 
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INTEGER*2 LGRP(14),MGRP(6) MC802. 





INTEGER*2 RETTBL(MXX, 3) MC802% 

REAL RETINV(MXX) MC802% 
NCEVAL=0 + MC802% 

IFC IGX. EQ.0 . AND. LGX. EQ. O . AND. MGX.E0.0) RETURN MC802€ 
TAINV=0. 0 MC802€ 

DO 100 I=1,NRET MC802€ 

C --* SCREEN ON YCS MC802€ 
DO 10 J=1,NYCSX MC802€ 
IF(RETTBL(I,2) TEOM 16: AT) GO Ta 15 MC802E 

10 CONTINUE MC802€ 

GO TO 100 MC802$ 

C === SCREEN ON MOS BY GROUP, L.GRP OR MG DEPENDING ON IGX,LGX,MGX MC8029 
15 CONTINUE MC8029 
MOS=RETTBL(I,1) MC8029 
IGP=RETTBL(I,3) MC8029 
LGP=LGRP(IGP) MC802$ 

IF(MGX .GT. 0) THEN MC8029 
IF(MGRP(LGP) .EQ. MGX) THEN MC8025 

IF(LGP .NE. LG) GO TO 80 MC8029 

ENDIF MC8029 

ELSE IF(LGX .GT. 0) THEN MC8030. 

IF(LGP .EQ. LGA) THEN MC8030 

IF( IGP .NE. IGR) GO TO 80 MC8030 

ENDIF MC8030. 

ELSE MC8030. 

IF(IGP .EQ. IGX) GO TO 80 MC8030 

ENDIF MC8030. 

GO TO 100 MC8030. 

80 CONTINUE MC8030. 

C === ACCEPTED MC8030: 
IF(RETINV(I) .GE. AIMIN) THEN MC8031 
NCEVAL=NCEVAL+1 MC8031 

ELSE MC8031. 
TAINV=TAINV+RETINV( 1) MC8031. 

ENDIF MC8031. 

100 CONTINUE MC8031i 

C --- FINAL ESTIMATE 15 NGEVAL MC8031! 
IF(AIMIN. GT. 0) NCEVAL=NCEVAL + AGGPCT*TAINV/AIMIN MC80311 

END MC80311 

C MC8031: 
ee dee dee de de dee e ele de de ee ee dee ee ede de De deed ee ede de ede ede te dee de de ede de ee de dee dee ede deed [01006 A 
C MC8032 
SUBROUTINE GETVYC(SYCS, LG, YCSB,NYE,NYB,NYEG, VYC) MC8032: 
INTEGER*2 YCSB(NYE,NYB,NYEG), VYC(NYE), LGEX(6) MC8032: 
INTEGER SYCS MC8032: 

DATA LGEX/4,4,1,2,4,3/ MC8032: 

C --- L INDICATES LAST DIMENSION IN YCS EXPANSION TABLE MC8032€ 
L=LGEX( LG) MC80328 

C --- FIND TO WHICH YCS BLOCK SYCS BELONGS AND MAKE COPY IN VYC MC8032t 
DO 10 J=1,NYB MC8032¢ 

DO 20 I=1,NYE MC8033( 
IF(SYCS . EQ. YCSBOR EN MC8033: 

DO 30 K=1,NYE MC80334: 
VYC(K)=YCSB(K,J,L) MC80333 


30 CONTINUE MC8033¢ 
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20 
10 


RETURN 
ENDIF 
CONTINUE 
CONTINUE 


WRITE(6,*) "kee YCS NOT FOUND IN YCSB TABLE YCS=' ,SYCS 


END 


MC803350 
MC803360 
MC803370 
MC803380 
MC803390 
MC803400 
MC803410 


RAAK IK dek eck dee kee RR EIE Re RR RR RR KIKKER RR RR RR RR KERK ze ie ze ese s MC 8 O 3 2 0 


M = ECH TABLE WITH ALL EXISTING COMBINATIONS FOR SELECTION CRITERIA 
C --- ACCEPT RECS WITH MATCHING PG,MG,CS,SVC. 


C 


20 


21 


30 


10 
999 


100 


a READET(RETTBL,RETINV,MXX,NRET, SGRD, NSC,SVCMP, 


NCSR,SCSRC, MG, LGRP, MGRP, MOSGR ,NMS) 


ACCEPT ALL YCS 
INTEGER SVCMP(5), NSC, SVC 
INTEGER SCSRC(16),NCSR, CS 


INTEGER SGRD, PG = 


INTEGER MOS,YCS 

INTEGER*2 MOSGR(2,NMS), MGRP(*) ,LGRP(*) 
INTEGER*2 RETTBL(MXX, 3) 
REAL RETINV(MXX), AI 
NRET=0 
DO 10 1=1,999999 

READ( 10+SGRD, 100, END=999) PG,MOS,YCS,SVC,CS, NRECS,AI 
IF(PG .NE. SGRD) GO TO 10 
IGR=IGFIND(MOS, MOSGR,NMS) 

LG=LGRP( IGR) 

IF(MGRP(LG) .NE. MG) GO TO 10 

DO 20 J=1,NSC 

IF(SVC . EQ. SVCMP(J)) GO TO 21 

CONTINUE 

Cor Tor 10 

CONTINUE 

DO 30 J=1,NCSR 

IF(CS .EQ. SCSRC(J)) THEN 
CALL ACCEPT(MOS,YCS,IGR,RETTBL,MXX ,NRET,RETINV,AI) 
Co TO 10 
ENDIF 
CONTINUE 


CONTINUE 

CONTINUE 

IF(NRET .GT. MXX) THEN 
WRITE(6,7) '**:** ERROR - TOO MANY RECORDS IN RETTBL' 
STOP 

ENDIF 

FORMAT(12,14,13,12,13,14,F7. 2) 

END 


MC803430 
MC803440 
MC803450 
MC803460 
MC803470 
MC803480 
MC803490 
MC803500 
MC803510 
MC803520 
MC803530 
MC 803540 
MC803550 
MC803560 
MC803570 
MC803580 
MC803590 
MC803600 
MC803610 
MC803620 
MC803630 
MC803640 
MC803650 
MC803660 
MC803670 
MC803680 
MC803690 
MC803700 
MC803710 
MC803720 
MC803730 
MC803740 
MC803750 
MC803760 
MC803770 
MC803780 
MC803790 
MC803800 
MC803810 
MC803820 


dede dede de de dede de dede de de dee e de de dede de de de e IE MK e de de e ED KM EE Ee EE EA TEM EE EE IE EE III MC LO03830 


C 


C= == AQEBPI ENTRY. 


SUBROUTINE ACCEPT(MOS ,YCS,IGR, RETTBL,MXX ,NRET,RETINV,AI) 
ACCUMULATE IF ALREADY SAME COMBINATION IS PRESENT 


INTEGER MOS,YCS 
INTEGER 2SRETTBECHRX, 3) 
REAL RETINV(MXX), AI 

DO 10 I=1,NRET 


S7 


MC803840 
MC803850 
MC803860 
MC803870 
MC803880 
MC803890 
MC803900 


IF(MOS. EQ. RETTBL(I,1) . AND. YCS. EQ. RETTBL( 1,2) ) THEN 


RETINV(I)=RETINV(I) + AI 
RETURN 
ENDIF 
10 CONTINUE 
C --- NEW COMBINATION 

NRET=NRET+1 

RETTBL(NRET, 1)=MOS 

RETTBL(NRET, 2)=YCS 

RETTBL(NRET, 3)=IGR 

RETINV(NRET)=AI 


ISFLAG) 


,SMOS( 1) 


END 
C 
keck kee e eck cke ck EE este EE EIER EE IE II IE EIE IE RR RE EK EIE MO 80 40 
ë 
, SUBROUTINE GETPAR(AIMIN,NO ,NMOS , SMOS,NYCS,SYCS,SGRD, 
NSC , SVCMP, NCSR, SCSRC, IGR, MOSGR, NMS, 
C --- e. SELECTION CRITERIA FROM USER AND VALIDATE 
INTEGER SYCS(31), NYCS 
INTEGER SMOS(20), NMOS 
INTEGER SVCMP(5), NSC 
INTEGER SCSRC(16),NCSR 
INTEGER SGRD 
INTEGER*2 MOSGR(2,NMS) 
WRITE(5,*) ' ENTER THRESHOLD MIN. FOR AVERAGE INVENTORY’ 
READ(5,*) AIMIN 
WRITE(5,*) ' ENTER THRESHOLD MIN. FOR NUMBER OF CELLS' 
READ(5,*) NO 
WRITE(5,*) ' THRESHOLDS TO USE AIMIN, NO=' ,AIMIN,NO 
C 
WRITE(5,%) ' ENTER MOS (ONLY 1 ACCEPTED)' 
NMOS=1 
READ(5,*) SMOS(1) 
WRITE(6,*) ' MOS SELECTED: ', SMOS(1) 
IGR=IGFIND( SMOS(1), MOSGR,NMS) 
WRITE(6,*) ' GROUP TO USE:', IGR 
IF(IGR. EQ. 0) THEN 
WRITE(S,*) '*#***%* ERROR - INVALID MOS SELECTED: ' 
STOP 
ENDIF 
E 
WRITE(5,*) ' ENTER YCS (ONLY 1 ACCEPTED) ' 
NYCS=1 
READ(5,*) SYCS(1) 
WRITE(6,*) ' YCS SELECTED: ', SYCS(1) 
C 
WRITE(5,*) ' ENTER GRADE' 
READ(5,*) SGRD 
WRITE(6,*) ' GRADE SELECTED', SGRD 
E 


WRITE(5,*) ' ENTER NO. OF SVC. COMPS & ARRAY (1-3, 4=1+2, S=ALL)' 


READ(5,%) NSC, (SVCMP(I), I=1,NSC) 
C --- EXPAND 4 TO 1,2 AND 5 TO 1,2,3 
DO 10 I=1,NSC 
IF(SVCMP(I).EQ.4 .OR. SVCMP(I).EQ.5) THEN 
NSC=SVCMP(I)-2 
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. MC8041! 






MC8 040. 


MC8040. 
MC8040: 
MC8040. 
MC8040 
MC8040 
MC8041! 
MC8041 
MC8041. 


MC804 1: 
MC8041. 
MC804 1 
MC8041. 
MC80411 
MC804 1! 
MC8042( 
MC8042: 


MC80439; 
MC80440. 
MC8044 1' 
MC80442 
MC80443 
MC80444 
MC80445 
MC80446 


15 
10 


dil 


CO soe 


22 


20 
26 


C --- 


C ess 


C 
101 
102 
103 
104 


C 


DO 15 J=1,NSC 
SVCMP( J)=J 
CONTINUE 
GEE 11 
ENDIF 
CONTINUE 
CONTINUE 
WRITE(6,*) ' SERVICE COMPONENTS SELECTED', (SVCMP(I), I=1,NSC) 


WRITE(5,*) ' ENTER NO. COMM. SOURCES AND ARRAY (1-15, 16=ALL)' 
READ(5,*) NCSR, (SCSRC(1), I=1,NCSR) 
IF 16 IS SELECTED THEN EXPAND ARRAY TO COVER ALL 1-15 
DO 20 I=1,NCSR 
IF(SCSRC(I) .EQ. 16) THEN 
NCSR=15 
DO 25 J=1,NCSR 
SCSRC( J)=J 
CONTINUE 
GO TO 26 
ENDIF 
CONTINUE 
CONTINUE 
WRITE(5,*) ' COMM. SOURCES SELECTED: ', (SCSRC(I), I=1,NCSR) 


MC804470 
MC804480 
MC804490 
MC8 04500 
MC804510 
MC804520 
MC804530 
MC804540 
MC804550 
MC804560 
MC804570 
MC804580 
MC804590 
MC804600 
MC804610 
MC804620 
MC804630 
MC804640 
MC804650 
MC804660 
MC804670 
MC804680 
MC804690 
MC804700 


FLAG TO DETERMINE WHICH OF SVC OR CS WILL BE USED AS 3RD DIMENSIONMC804710 
WRITE(5,*) 'SELECT 3RD DIM. TO USE: O=NONE, l=SVC, 2=COMM. SOURCE 'MC804720 


READ(5,*) ISFLAG 
WRITE INPUT PARAMETER INFO TO OUTPUT FILE 
WRITE(11,101) 'TEST CASE INPUT PARAMETERS: ' 
WRITE(11,102) 'INVENTORY THRESHOLD= ",AIMIN, 

'THRESHOLD NO. OF CELLS= ",NO 
WRITE(11,103) 'MOS= ',SMOS(1),'YCS= ',SYCS(1),'GRADE= ',SGRD 
WRITE(11,104) 'SERVICE COMPONENTS= ',(SVCMP(I),I=1,NSC) 
WRITE(11,104) 'COMM SOURCES= ',(SCSRC(I),I=1,NCSR) 
WRITE(6,*) | '3RD DIMENSION= ' , ISFLAG 


FORMAT( 1X, A) 

FORMAT(1X,A,F4.1,7X,A,12) P 
FORMAT(1X,A,I3,2(5X,A,I2)) 

FORMAT(1X,A,15(13)) 

END 


MC804730 
MC804 740 
MC804750 
MC804760 
MC804770 
MC804780 
MC804790 
MC804800 
MC804810 
MC804820 
MC804830 
MC804840 
MC804850 
MC804860 
MC804870 
MC804880 


dede dede dede dede dede de dede dedo de dede dede dede dede dede dede de de dede dese dede de dede de dede de dede dee dee dese dede EIE IERE KI IESER III MO BOLROYO 


C 


, 
SÉ 


C =-~ 


C -=~~ 


SUBROUTINE GETMOS(SMOS ,NMOS ,MGX,LGX,MG,LG,IGR,MOSGR,LGRP,MGRP, 
NMS , NG , NLG) 
BUILD SMOS ARRAY BASED UPON EXPANSION 
INTEGER SMOS(30) 
INTEGER*2 MOSGR(2,NMS), LGRP(NG), MGRP(NLG) 
NMOS=0 
IF(MGX .GT. 0) THEN 
HAVE EXPANDED TO MAJOR MOS GROUP 
DO 10 I=1,NMS 
IGP=MOSGR( 2,1) 
LGP=LGRP( IGP) 
IF(MGRP(LGP) .EQ. MG) THEN 
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MC804900 
MC804910 
MC8 04920 
MC804930 
MC804940 
MC804950 
MC804960 
MC804970 
MC804980 
MC804990 
MC805000 
MC805010 
MC805020 


NMOS=NMOS+1 MC805( 











SMOS(NMOS )=MOSGR( 1,1) MC805( 

ENDIF MC805( 

10 CONTINUE MC805( 
RETURN MC805( 

ELSE IF(LGX .GT. 0) THEN MC805¢ 

C --- HAVE EXPANDED TO LARGE MOS GROUP MC805C 
DO 20 I[=1,NMS MC805] 
IGP=MOSGR( 2,1) MC8051 

IF(LGRP( IGP) .EQ. LG) THEN MC8051 
NMOS=NMOS+1 MC8051 

SMOS( NMOS )=MOSGR( 1,1) MC8051 

ENDIF MC8051 

20 CONTINUE MC8051 
RETURN MC8051 

ELSE MC8051 

C --- HAVE EXPANDED TO SMALL MOS GROUP MC8051! 


DO 30 I[=1,NMS MC8052. 
IF(MOSGR(2,1) . EQ. IGR) THEN MC8052. 


NMOS=NMOS+1 MC8 052: 

SMOS(NMOS )=MOSGR(1,1) MC8052. 

ENDIF MC8052. 

30 CONTINUE MC8052: 
RETURN MC8052! 
ENDIF MC8052' 


END MC8052: 


MC805 2! 
dever de de dede dede dede dede dede dee de dede dede dede dede dede de dede dede dede dede dede dede dede dede dede de dede dede dede de de dede de dede dede dede dede de dede Ú ( 01210 ho BL 
C MC8 053) 

FUNCTION IGFIND(MOS, MOSGR,NMS) MC8053. 

C --- FIND LOCATION OF MATCHING MOS IN GROUP TABLE. RETURN GROUP NO MC805 34 
INTEGER*2 MOSGR(2,NMS) MC80534 

DO 10 I=1,NMS MC8053: 
IFCMOSGRCL 1). EO: MH 050 EIN MC8053€ 
IGFIND=MOSGR( 2,1) MC805 34 

RETURN MC8053€ 

ENDIF MC80535 

10 CONTINUE MC8054C 
IGFIND=0 MC80541 

END MC80542 

C MC80543 
ERREICHT Kaal (s IT ATA 


C MC80545: 
SUBROUTINE READER( DATA, INV ,Y ,MXX ,NMOS ,NYCSG,NYCSL,NYCSM,NSC,NCSR, MC80546 

* NYR,SMOS,SYCSG,SYCSL,SYCSM,SGRD,SVCMP,SCSRC,NRC,PTRTBL,LGX,MGX, MC80547 

* IGR,LG,MGRP,LGRP,MOSGR,NMS,NG,NLG, ICYCLE,NPT,PTBL, ISFLAG,SINV,SY)MC80548: 


REAL INV(MXX,NYR),Y(MXX,NYR), SINV(MXX,NYR),SY(MXX,NYR) MC80549; 
INTEGER*2 PTRTBL(MXX, 2), PTBL(MXX,3) MC80550: 
INTEGER SYCSG(*), SYCSL(*), SYCSM(*) MC805511 
INTEGER SMOS(*), NMOS MC80552: 
INTEGER SVCMP(**), NSC MC80553 
INTEGER SCSRC(*) ,NCSR MC80554 
INTEGER SGRD MC80555 
INTEGER TYPE, YCS,PG,MOS,SEX,CS,EDLV,SVC,MOS1,MOS2, RACE“ MC8055 6! 
INTEGER DATA(NYR) MC80557! 
CHARACTER*7 CITLS MC805581 
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INTEGER*2 MOSGR(2,NMS),LGRP(NG),MGRP(NLG) 


C 
C --- REWIND DATA FILE AND RESET INV,Y IF CYCLING THRU READER 
IF(ICYCLE .GT. 1) THEN 
REWIND 1 
DO 6 I=1,MXX 
DO 5 J=1,NYR 
INV(1,J)=0. 0 
Y(1,J)=0.0 
SINV(1,J)=0. 0 
SY(1,J)=0. 0 
5 CONTINUE 
6 CONTINUE 
ENDIF 
C --- READ RECORD AND STORE IN MATRIX 
ICR=0 
NRC=0 
NPT=0 
ICNT=0 
IYNO=0 
IYR=0 
C 
1 READ(1,101,END=999) TYPE,YCS,PG,MOS,SEX,CS,EDLV,SVC,MOS1,M0S2, 
d RACE ‚CITLS ‚DATA 
ICR=ICR+1 
C --- CHECK IF RECORD MEETS SELECTION CRITERIA. OTHERWISE REJECT. 
C --- COLLECT TYPES O=INVENTORY, AND 1-5 ALL LOSSES 
IF(TYPE. GT. 5) GO TO 999 
C 
C --- SCREEN FOR GRADE 
IF(PG .NE. SGRD) GO TO 1 
C 
C --- SCREEN FOR MOS 


IGP=IGFIND(MOS,MOSGR,NMS) 
IECIGP EO O) CO TO l 
LGP=LGRP( IGP) 
IF(MGX .GT. 0) THEN 
C --- HAVE EXPANDED TO MAJOR MOS GROUP 
IF(LGP .EQ. LG) THEN 
DO 10 I=1,NYCSL 
IF(YCS .EQ. SYCSL(1)) THEN 
it 
GO TO 60 
ENDIF 
10 CONTINUE 
CO TO 1 
ELSE IF(MGRP(LGP) .EQ. MGX) THEN 
DO 20 I=1,NYCSM 
IF(YCS . EQ. SYCSM(1)) THEN 
II 
GO TO 60 
ENDIF 
20 CONTINUE 
Go Toi 
ELSE 
GO TO 1 
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MC805590 
MC805600 
MC805610 
MC805620 
MC805630 
MC805640 
MC805650 
MC805660 
MC805670 
MC805680 
MC805690 
MC805700 
MC805710 
MC805720 
MC805730 
MC805740 
MC805750 
MC805760 
MC805770 
MC805780 
MC805790 
MC805800 
MC805810 
MC805820 
MC805830 
MC805840 
MC805850 
MC805860 
MC805870 
MC805880 
MC805890 
MC805900 
MC805910 
MC805920 
MC805930 
MC805940 
MC805950 
MC805960 
MC805970 
MC805980 
MC805990 
MC806000 
MC806010 
MC806020 
MC806030 
MC806040 
MC806050 
MC806060 
MC806070 
MC806080 
MC806090 
MC806100 
MC806110 
MC806120 
MC806130 
MC806140 


ENDIF 
ELSE IF(LGX .GT. 0) THEN 
C --- HAVE EXPANDED TO LARGE MOS GROUP 
IF(IGP .EQ. IGR) THEN 
DO 30 I=1,NYCSG 
IF(XCS .EQ. SYCSG(I)) THEN 
IY=I 
GO TO 60 
ENDIF 
30 CONTINUE 
CO TO Í 
ELSE IF(LGP .EQ. LGX) THEN 
DO 40 I=1,NYCSL 
IF(YCS . EQ. SYCSL(I)) THEN 
IY=I 
GO TO 60 
ENDIF 
40 CONTINUE 
CO TO 
ELSE 
EO rom 
ENDIF 
ELSE 
C --- HAVE EXPANDED TO SMALL MOS GROUP 
IF(IGP .EQ. IGR) THEN 
DO 50 I=1,NYCSG 
IF(XCS . EQ. SYCSG(I)) THEN 
IY=I 
GO TO 60 
ENDIF 
50 CONTINUE 


ENDIF 
60 CONTINUE 


DO 70 I=1,NMOS 
IF(MOS .EQ. SMOS(1)) THEN 
IM=I 
GO TO 80 
ENDIF 
70 CONTINUE 
WRITE(6,*) '*** ERROR IN MOS SCREENING ***' MOS 
WRITE(6,*) ‘NMOS,SMOS=' ,NMOS,(SMOS(I),I=1,NMOS) 
GO TO 1 


C --- SCREEN FOR SERVICE COMPONENT 
80 CONTINUE 
DO 90 I=1,NSC 
Tr SVGa EO) SVCME(IT)) THEN 
IS=I 
GO TO 100 
END IF 
90 CONTINUE 
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MC8061 





MC8063. 
MC8063 
MC8063. 
MC8063: 
MC8063 


fae = 
100 


110 


120 


C --- 


999 


C --- 


GO TO 1 


SCREEN FOR COMMISSIONING SOURCE 
CONTINUE e 
DO 110 I=1,NCSR 
IF(CS .EQ. SCSRC(I)) THEN 
IR=I 


GO TO 120 
END IF 
CONTINUE 
GO TO 1 


CONTINUE 
RECORD ACCEPTED - INSTALL IT IN INV,Y,SINV,SY, PTRTBL AND PTBL 


ICNT=ICNT+1 
IF(ISFLAG. EQ. 1) THEN 


IW=IS 

ELSE IFC ISFLAG. EQ. 2) THEN 
IW=IR 

ELSE 
IW=-99 

ENDIF 


MINV=GINV(PTRTBL, MXX,NRC, IM,IY,-99) 
MV=GINV(PTBL, MXX,NPT, IM, IY,IW) 
IF(TYPE. EO. 0) THEN 
CALL INSINV(PTRTBL,MXX,NYR,NRC,MINV, IM, IY, -99, INV, DATA) 
CALL INSINV(PTBL, MXX,NYR,NPT,MV, IM,IY, IW,SINV,DATA) 
ELSE 
CALL INSY(MXX,NYR,MINV,Y,DATA) 


CALL INSY(MXX,NYR,MV, SY. DATA) 
IYR=IYR+1 
IF(MINV. EQ. 0) THEN 
WRITE(6,*) '*** ERROR IN DATA BASE. LOSS W/O INV. REC. ' 


WRITE(6,122) 'Y**:', MOS,YCS,PG,EDLV,SVC,RACE, 
(DATA(IT), IT=1,NYR) 
IYNO=IYNO+1 


ENDIF 
ENDIF 
GO TO 1 
CONTINUE 
WRITE(6,%) ' ' 
WRITE(6,*) 'TOTAL RECORDS READ ER 
WRITE(6,*) 'TOTAL INV. MOS/YCS COMBINATIONS=' ,NRC 
WRITE(6,*) 'TOTAL INV. MOS/YCS/IW COMBINATIONS=' ,NPT 
WRITE(6,*) ‘TOTAL RECORDS ACCEPTED ENT 
WRITE(6,*) 'TOTAL LOSS RECORDS ACCEPTED IR 
WRITE(6,7) 'TOTAL LOSS RECORDS NOT MATCHED =) NO 


TERMINATE IF NO DATA COLLECTED 

IF(NRC .EQ. 0) THEN 
WRITE(6,*) '**** NO DATA MEETS SELECTION REQS' 
STOP 

ENDIF 
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MC806710 
MC806720 
MC806730 
MC806740 
MC806750 
MC806760 
MC806770 
MC806780 
MC806790 
MC806800 
MC806810 
MC806820 
MC806830 
MC806840 
MC806850 
MC806860 
MC806870 
MC806880 
MC806890 
MC806900 
MC806910 
MC806920 
MC806930 
MC806940 
MC806950 
MC806960 
MC806970 
MC806980 
MC806990 
MC807000 
MC807010 
MC807020 
MC807030 
MC807040 
MC807050 
MC807060 
MC807070 
MC807080 
MC807090 
MC807100 
MC807110 
MC807120 
MC807130 
MC807140 
MC807150 
MC807160 
MC807170 
MC807180 
MC807190 
MC807200 
MC807210 
MC807220 
MC807230 
MC807240 
MC807250 
MC807260 


101 ° EORMAT( 312,15. 1414. Pr OSS E) MC807. 








121 FORMAT(A8,1316) MC807: 
122 FORMATCAS,7I6, 5%, Io) MC8O7; 
131 FORMAT(I4, 216) MC807: 
132: FORMATCI4, 31661050) MC807. 
END MC807: 
MC807: 


eke ee ee Je deve deve deve sete We He Fete Fete qe Je Je ge ge qege gee ge ge fe ge ge ge gee ge gee gee qe fe ese gee e Je ge ge fe ge gede gee ge ke je ge g7e 7e ge3e 3esesçM C 8 O 7 ° 


C MC807: 
FUNCTION GINV(PTBL, MXX,NPT, IM,IY,IW) MC807: 
C --- FIND LOCATION OF INVENTORY ENTRY FOR MOS,YCS,SVC/CS COMBINATIONS MC807: 
C --- 3RD DIMENSION CHECKED ONLY IN CASE IW>0 MC807: 
INTEGER*2 PTBL(MXX, *) MC807: 


DO 10 I=1,NPT MC80 74 


IF(PTBL(I, 1) .EQ. IM . AND. ' MC8074 

* PTBL(I, 2) .EQ. IY ) THEN = MC8074 
IFCIW.LT.O .OR. (IW.GT.O .AND. PTBL(I, 3).EG. IW)) THEN ` MC8074 

GINV=I MC8074 

RETURN MC8074 

ENDIF MC8074 

ENDIF MC8074 

10 CONTINUE MC8074 


GINV=0 MC8074 
END MC8075. 


C MC8075 
kkikikikikikikikikiki kiżiżikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikikiki żidu (01010 W ASS 
C MC 8075: 
SUBROUTINE INSINV(PT,MXX,NYR,N,K,IM,IY,IW,INV,DATA) MC8075: 

C --- ACCUMM INTO KTH ENTRY. INSTALL IN POINTER TABLE IF NOT PRESENT MC8075 
REAL INV(MXX, NYR) MC8075' 
INTEGER*2 PT(MXX, *) MC8075: 
INTEGER DATA(NYR) MC8075: 

IF(K .EQ. 0) THEN MC8075 

C sos ADD NEW ENTRY MC8076: 
=N+1 MC8076: 

IF(N .GT. MXX) THEN MC8076: 

WRITE(6,*) '**>* ERROR - TOO MANY INV. COMBINATIONS' ,N MC8076: 

STOP MC8076: 

ENDIF MC8076: 

=N MC80761 

PT(K, 1)=IM MC8076: 

PT(K, 2)=IY MC8076: 

IF(IW.GT.O) PT(K, 3)=IW MC807641 

ENDIF MC8077( 

DO 130 IT=1,NYR MC8077: 
INV(K,IT)=INV(K,IT) + . 25*FLOAT(DATA(IT)) MC8077: 

130 CONTINUE MC8077: 
END MC807 7: 


G MC8077: 
Tetede dede dede de dede dede dede dede dede dede dede dedede dede dedo dede dede dede dede dede de dede dede de dede dede dede de dede dede de dede dede dede dede de detector MCB8O0 776 
C MC8077, 
SUBROUTINE INSY(MXX,NYR,K,Y,DATA) MC80771 
C --- ACCUMM INTO KTH ENTRY FOR LOSS MC8077€ 
REAL Y(MXX, NYR) MC8078( 
INTEGER DATA(NYR) MC8078) 


IFFR) RETURN MC80784. 
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DO 10 IT=1,NYR MC807830 

SR UT) = + DATA(IT) MC807840 

10 CONTINUE MC807850 
END MC807860 

C MC807870 
kkkkkkkkkkliklkiðklklkevklkklkkkkilikkklkilkiaoklkiklikkklilikkklklkiklkikilkkiklilklkUlklik kk ngMC 807880 
C MC807890 
_ DEERE AGGREG( INV, Y,MXX, NYR, SMOS, SYCSG, MC807900 

NRC, NRCOLD, PTRIBL, INDX, AVINV, AIMIN,MKG) MC807910 

> See Som AVERAGE INV. & SORT MC807920 
REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) MC807930 
INTEGER*2 PTRTBL(MXX, 2), INDX(MXX) ,MKG(MXX) MC807940 
INTEGER SYCSG(*), SMOS(*) MC807950 
REAL*8 TINV,TY MC807960 

C MC807970 
C --- RESET MKG (NECESSARY WHEN CYCLING THRU AGGPCT VALUES) MC807980 
DO 10 [=1,MXX MC807990 
MKG(I)70 MC808000 

10 CONTINUE MC808010 
TINV=0 MC808020 

TY=0 MC808030 

DO 100 I=1,NRC MC808040 

ba-—— FIX INV. ENTRIES LOWER THAN CORRESP. LOSSES & COMP. AVG INV. MC 808050 
AI=0 MC808060 

DO 201 J=1,NYR MC808070 
TINV=TINV+INV(1,J) a MC808080 

Tye TYE YET i MC808090 
IF(INV(I,J).LT.Y(I,J)) INV(I,J)=Y(I,J) MC808100 
AI-AIHINV(I,J) MC808110 

201 CONTINUE MC808120 
AVINV(I)=AI/NYR MC808130 

INDX( 1)=1 MC808140 

100 CONTINUE MC808150 
WRITE(6,*) '===== TOTAL INV,Y=',TINV,TY MC808160 

C MC808170 
C --- SORT ASCENDING BY AVG INVENTORY MC808180 
CALL SORT2(AVINV, INDX,NRC) MC808190 

C MC808200 
NS1=0 MC808210 

C --- DISPLAY TABLE IN SORT SEQUENCE MC808220 
CALL DSPTBL(INV,Y,AVINV,PTRTBL,INDX,AIMIN,NRC,MKG,MXX,NYR, MC808230 
SYCSG, SMOS ) MC808240 

C MC808250 
DO 200 K=NRC,1,-1 MC808260 
IF(AVINV(K) . GE. AIMIN) THEN MC808270 

Doo MARK AS MEMBER OF SET 50 MC808280 
MKG(K)=32767 MC808290 

ELSE MC808300 

GM = INITIAL COUNT OF MEMBERS OF SET Si MC808310 
NS 1=K MC808320 

GO TO 202 MC808330 

ENDIF MC808340 

200 CONTINUE MC808350 
202 CONTINUE MC808360 

C --- DO AGGREGATIONS WITHIN SET Sl UNTIL NO MORE POSSIBLE (KF GE 0) MC808370 
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KF=-1 MC808: 


--- DO WHILE KF<0 MC808: 
300 IF(KF.GE.0) GO TO 310 MC8084 
CALL AGGI(AVINV,INDX,MKG,NS1, INV,Y,MXX,NYR, AIMIN,KF) MC808/ 

GO TO 300 MCBO84 

310 CONTINUE MC8084 
--- DISPLAY TABLE AFTER 1ST AGGREGATION MC808: 
CALL DSPTBL(INV,Y,AVINV,PTRTBL, INDX,AIMIN,NRC,MKG,MXX,NYR, MC8084 

* SYCSG,SMOS ) MC8084 
IF(NS1. EQ. NRC) THEN MC8084 
WRITE(6,%) '***%*%* SET SO EMPTY. NO CELLS ABOVE THRESHOLD' MC8084 

STOP MC8084 

ENDIF MC808£ 

--- DO AGGREGATIONS FROM SET 51 INTO SET 50 UNTIL NO MORE POSSIBLE MC808£ 
KF=1 MC8085 

--- DO WHILE KE20 MC808£ 
320 IF(KE. LE AGO MEONTO 0 MC8085 
CALL AGG2(AVINV,INDX,MKG,NS1,NRC,INV,Y,MXX,NYR, KF) MC8085 

GO TO 320 MC8085 

330 CONTINUE MC8085 
--- DISPLAY TABLE AFTER 2ND AGGREGATION MC8085 
CALL DSPTBLCINV,Y,AVINV,PTRTBL, INDX,AIMIN,NRC ,MKG ,MXX,NYR, MC8085 

* SYCSG,SMOS ) MC808€ 
--- MOVE VALUES GE AIMIN TO BEGINNING OF ARRAYS MC8086 
CALL CMPRS( INV,Y ,MXX,NYR,NRC ,NRCOLD , AIMIN, AVINV) MC8086 

--- DISPLAY TABLE AFTER MOVING VALUES. MC8086 
DO 400 K=1,NRC MC8086 
WRITE(6,122)K,AVINV(K), (INV(K,J),J=1,NYR) MC8086 
WRITE(6, 129) elt che MC8086 

400 CONTINUE MC8086 
122 FORMAT(/15,14X,F8.3, 6X, 10672) MC8086 
123 FORMAT( Do 110 7.2) MC8086 
END MC8087 

e gl ele dee oe de dee gie ee ce MC8087 


SUBROUTINE AGGI(AVINV, INDX,MKG,NS1,INV,Y,MXX,NYR,AIMIN,KF) MC8087 






--- DO ONE PASS OF AGGREGATION MC8087 
REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) MC8087; 
INTEGER*2 INDX(MXX) , MKG(MXX) a MC8087: 
KF=0 MC80871 
CI=0 MC8087: 
DO 10 I=NS1,1,-1 MC8087) 

IF(MKG(1). EO. 0) THEN MC8087;| 


IF(KF.EQ.0) THEN MC8088: 

--- THIS WILL BE THE COLLECTING CELL MC8088: 
KF=I MC8088: 
CI-AVINV(I) MC8088: 

ELSE MC8088) 
IF(CI+AVINV( I). LT. AIMIN) THEN MC8088! 

--- ACCUM. WITH CELL KF TEMPORARILY. SET TEMP. POINTER -KF MC8088: 
CI=CI+AVINV( I) MC8088: 

MKG( I )=-KF MC8088. 

ELSE MC8088' 

= FIND SMALLEST CELL TO ADD MC8089! 
CALL AGGIA(AVINV,MKG,I,CI,AIMIN,KF,MXX) MC8089 

ENDIF MC8089: 
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IF(CI. GE. AIMIN) THEN 
MAKE THIS AGGREGATION PERMANENT AND EXIT 
AVINV(KF)=CI 
CALL AGGIB(INDX,MKG,KF , INV,Y,NYR,MXX) 
NS1=NS1-1 
MKG(KF)=32767 
KF=-1 
RETURN 
ENDIF 
ENDIF 
ENDIF 
10 CONTINUE 
C 
IF(KF.EQ.O) RETURN 
C --- CLEAR TEMPORARY POINTERS LEFT. 
DO 20 1=1,NS1 
IF(MKG( 1). LT.0) MKG(1)=0 
20 CONTINUE 
END 
C Tee Dee dee ee Dee ee eer 
SUBROUTINE AGGIA(AVINV,MKG,ILAST,CI,AIMIN,KF,MXX) 
C --- FIND SMALLEST CELL TO ADD AND SET TEMPORARY POINTER 
REAL AVINV(MXX) 
INTEGER*2 MKG(MXX) 
DO 10 I=1,ILAST 
IF(MKG(I).EQ.O) THEN 
IF(CIHAVINV(I) .GE. AIMIN) THEN 
CI=CI+AVINV(I) 
MKG(I)=-KF 
RETURN 
ENDIF 
ENDIF 
10 CONTINUE 
WRITE(6,*) 
STOP 
END 
C EE EE e de kk dee decke 
SUBROUTINE AGGI1B( INDX,MKG,KF,INV,Y,NYR,MXX) 
C --- MAKE AGGREGATION PERMANENT 
REAL INV(MXX, NYR), Y(MXX, NYR) 
INTEGER*2 INDX(MXX) ,MKG(MXX) 
K-INDX(KF) 
DO 10 I=1,KF-1 
IF(MKG(I) .LT. 0) THEN 
IF(MKG(1).NE. -KF) STOP 777 
MKG( I )=KF 
L=INDX( I) 
DO 20 J=1,NYR 
INV(K,J)=INV(K,J)+INV(L,J) 
Y(K,J)= Y(K,J)+ Y(L,J) 
20 CONTINUE 
ENDIF 
10 CONTINUE Së 
END 


e wl ml wl sl ele sl mis tl. Le se el, 
C Tee de ze dee dee e ee dee 


"Wi ERROR IN AGGIA. NO VALUE FOUND rt! 
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THIS WAS AN UNSUCCESSFUL AGGREG. 


MC808930 
MC808940 
MC808950 
MC808960 
MC808970 
MC808980 
MC808990 
MC809000 
MC809010 
MC809020 
MC809030 
MC809040 
MC809050 
MC809060 
MC809070 
MC809080 
MC809090 
MC809100 
MC809110 
MC809120 
MC809130 
MC809140 
MC809150 
MC809160 
MC809170 
MC809180 
MC809190 
MC809200 
MC809210 
MC809220 
MC809230 
MC809240 
MC809250 
MC809260 
MC809270 
MC809280 
MC809290 
MC809300 
MC809310 
MC809320 
MC809330 
MC809340 
MC809350 
MC809360 
MC809370 
MC809 380 
MC809390 
MC809400 
MC809410 
MC809420 
MC809430 
MC809440 
MC809450 
MC809460 
MC809470 


10 
12 


20 


30 


SUBROUTINE AGG2(AVINV, INDX,MKG,NS1,NRC, INV, 


Y,MXX,NYR, KF) 


DO ONE PASS OF AGGREGATION FROM SET 51 TO SET SO 


ON EACH PASS ONE ELEMENT OF S1 IS TAKEN & ADDED TO SMALLEST OF SO 


REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) 
INTEGER*2 INDX(MXX) , MKG(MXX) 
KF20 
FIND ELEMENT OF S1 
DO 10 I=1,NS1 
IF(MKG(1). EO. 0) THEN 
KF=I 
0) Wo) 
ENDIF 
CONTINUE 
CONTINUE 


(ONLY THOSE WITH POINTER MKG( 1)=0) 


IF KF STILL 0 THEN NO MORE ELEMENTS IN 51 LEFT 


IF( KF. EQ. 0) RETURN 


FIND SMALLEST ELEMENT OF SO AND ADD TO IT. 
ISM=NRC 
SMALL=AVINV( ISM) 
DO 20 I=1, NRC 
IF(MKG(1). EQ. 32767) THEN 
IF(AVINV(I). LT. SMALL) THEN 
ISM=I 
SMALL=AVINV( I ) 
ENDIF 
ENDIF 
CONTINUE 
JOIN ELEMENT KF TO ELEMENT ISM 
AVINV( ISM)=AVINV( ISM) + AVINV(KF) 
MKG(KF)=ISM 
L=INDX(KF) 
K=INDX( ISM) 
DO 30 J=1,NYR 
INVCK,J)=INVCK,J)+INV(L,J) 
Y(K,J)= Y(K,J)+ Y(L,J) 
CONTINUE 
END 


Sevevesesevevededevevedesedeve 


SUBROUTINE CMPRSC INV,Y,MXX,NYR,NRC,NRCOLD ,AIMIN, AVINV) 


REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) 
COMPRESS INV,Y IN PLACE. 
NRCOLD=NRC 

NRC=0 


DO 10 I=1,NRCOLD 


AI=CAINV( INV, I, MXX,NYR) 
IF(AI .GE. AIMIN) THEN 
TRANSFER ACTIVE CELL I ---> NRC 
NRC=NRC+1 
AVINV(NRC)=AI 
DO 20 J=1,NYR 
INV(NRC ,J)=INV(I,J) 


Y(NRC,J)= Y(I,J) 
20 CONTINUE 
ENDIF 
10 CONTINUE 
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ONLY WITH MKG(1)=32767 


MOVE ALL ROWS GE AIMIN TO TOP 





MC8096: 
MC8096: 
MC8096: 
MC8096: 


MC8097; 
MC8097: 
MC809 74 
MC8097: 


MC80996; 
MC80997' 
MC809981 
MC80999ı 
MC810001 
MC81001( 
MC810021 
MC810031 


END 
KATTA ek 
FUNCTION CAINV(INV,I,MXX,NYR) 
REAL INV(MXX, NYR) 
--- COMPUTE AVERAGE INVENTORY FOR ROW I 
CAINV=0 
DO 10 J=1,NYR 
CAINV=CAINV+INVC I,J) 
10 CONTINUE 
CAINV=CAINV/NYR 
END 
úk Te Te RS 
SUBROUTINE DSPTBL(INV,Y,AVINV,PTRTBL, INDX,AIMIN,NRC, MKG, MXX,NYR, 
* SYCSG,SMOS ) 
--- DISPLAY TABLE IN SORT SEQUENCE 
REAL INV(MXX, NYR), Y(MXX, NYR), AVINV(MXX) 
INTEGER*2 PTRTBL(MXX, 2), INDX(MXX), MKG(MXX>) 
INTEGER SYCSG(*) 
INTEGER SMOS(*) 
INTEGER IATT(2) 
CHARACTER*1 STI 
WRITE(6,121) 
WRITE(6,*) 'INV. THRESHOLD MIN. VALUE=" ,AIMIN 


WRITE(6,*) ' I INDX AVG MKG INVENTORY/LOSSES' 
DO 200 K=1,NRC 

STIs' ' 

I-INDXCK) 

AI=AVINV(K) 

IF(AI .LT. AIMIN) STIs'$' 

IATT(1)=SMOS(PTRTBL(I,1)) 

IATT(2)=SYCSG(PTRTBL(I,2)) 


MC810040 
MC810050 
MC810060 
MC810070 
MC810080 
MC810090 
MC810100 
MC810110 
MC810120 
MC810130 
MC810140 
MC810150 
MC810160 
MC810170 
MC810180 
MC810190 
MC810200 
MC810210 
MC810220 
MC810230 
MC810240 
MC810250 
MC810260 
MC810270 
MC810280 
MC810290 
MC810300 
MC810310 
MC810320 
MC810330 
MC810340 
MC810350 


WRITE(6,122)K,I,AI,MKG(K), STI, (INV(1,J),J=1,NYR) , (IATT(J),J=1,2) ,MC810360 


de PTRTBL(1,1),PTRTBL(1,2) 
WRITE(6,123) ( Y(I,J),J=1,NYR) 
200 CONTINUE 


121 FORMAT(///) 
122 FORMAT(/215,F8.3,19,1X,A2, 10F7.2, 5X, 615) 
123 FORMAT( 30X, 10F7.2) 
END 
ese ese ee ee Ese 
SUBROUTINE SORT2(Y,INDX, N) 
EOS SORT USING SHELL ALGORITHM “eens: 
--- SORTS ON Y AND DOES SAME REORDERING ON INDEXES INDX 
REAL Y(N) ,TEMP 
INTEGER GAP 
INTEGER*2 INDX(N), ITEMP 
LOGICAL EXCH 


GAP=(N/2) 
5 BENT (GAP NE.0)) GO TO 500 
10 CONTINUE 
EXCH=. TRUE. 
K=N-GAP 


DO 200 I=1,K 
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MC810370 
MC810380 
MC810390 
MC810400 
MC810410 
MC810420 
MC810430 
MC810440 
MC810450 
MC810460 
MC810470 
MC810480 
MC810490 
MC810500 
MC810510 
MC810520 
MC810530 
MC810540 
MC810550 
MC810560 
MC810570 
MC810580 
MC810590 


100 
200 


500 


C 


Fede dedede deve dese He seve se Weve qe eee qeq de ee de de eee dee de de dee dede de dede dede dee dee de de FE e ee e e ee kee He Fe Te TENE HE Ke KK IKMC 81 O7 | 


C 


KK=I+GAP 
IF(. NOT. (Y(1). GT. Y(KR))) GO TO 100 
TEMP=Y( 1) 
Y(1)=Y(KK) 2 
Y(KK)=TEMP 
ITEMP=INDX(I) 
INDX(I)=INDX(KK) 
INDX(KK)=ITEMP 
EXCH=. FALSE. 
CONTINUE 
CONTINUE 
IF (. NOT. (EXCH)) GO TO 10 
GAP=(GAP/2) 
GO TO 5 
CONTINUE 
RETURN 
END 


SUBROUTINE BKDOWN(PTBL,NPT,PTRTBL,NRC, INDX,MKG,MXX,MXY, 
*  SINV,SY,INV,Y,BKTBL,NBK 


C --- BREAKDOWN AGGREGATED VALUES BY THE 3RD DIMENSION SVC/CS 


C --- 


10 
Ca 


C --- 


ion 
102 
103 
104 


REAL INV(MXX,MXY) ,Y(MXX,MXY) , SINV(MXX,MXY),SY(MXX,MXY) 
INTEGER*2 PTRTBL(MXX, 2), INDX(MXX), MKG(MXX) 

INTEGER*2 PTBL(MXX, 3), BKTBL(MXX,3) 

REAL*8 TINV,TY 

NBK=0 

TRAVERSE MKG ARRAY AND BUILD BKTBL 

DO 10 I=1,NRC 

IF(MKG(I). NE. 32767) THEN 

ICELL=MKG( 1) 
ELSE 
ORI 

ENDIF 

IX-INDX(I) 

IM-PTRTBL(IX,1) 

IY=PTRTBL(IX,2) 

CALL BLDBK(ICELL,IM,IY,PTBL,NPT,MXX,BKTBL,NBK) 
CONTINUE 
DISPLAY BKTBL PRIOR TO SORTING 
WRITE(6,101) (I,(BKTBL(1,J),J=1,3), I=1,NBK) 
CALL SORT3(BKTBL,NBK,MXX) 
WRITE(6,101) (I,(BKTBL(1,J),J=1,3), I=1,NBK) 
SUMMARIZE SINV,SY INTO INV,Y FOR MATCHING ENTRIES IN BKTBL 
CALL SUMBK(BKTBL,NBK,MXX,SINV,SY,INV,Y,MXY) 
WRITE(6,102) (I,(INV(1,J),J=1,MXY),(BKTBL(1,J),J=1,2), I=1,NBK) 
WRITE(6,102) (1,( Y(1,J),J=1,MXY),(BKTBL(1,J),J=1,2), I=1,NBK) 
FORMAT(14, 316) 
FORMAT(I4, 10F7.2,10X,214) 
FORMAT(/15,10F7.2) 
FORMAT( 5X,10F7. 2) 
END 


€ kki rri kri 


SUBROUTINE BLDBK(ICELL, IM,IY,PTBL,NPT,MXX,BKTBL,NBK) 
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MC8 107 


MC81112! 
MC81113! 
MC811141 


INTEGER*2 PTBL(MXX, 3), BKTBL(MXX,3) MC811150 
C --- RECORD ALL ENTRIES IN PTBL WITH MATCHING IM,IY IN BKTBL MC811160 
DO 10 I=1,NPT MC811170 
IF(PTBL(1,1). EO. IM . AND. PTBL(1,2). EO. IY) THEN MC811180 

“E INSTALL WITH CELL ID, IW & POINTER MC811190 
NBK=NBK+1 MC811200 

BKTBL(NBK, 1)=ICELL MC811210 

BKTBL(NBK, 2)=PTBL(I,3) MC811220 

BKTBL(NBK, 3)=1 MC811230 

ENDIF MC811240 

10 CONTINUE MC811250 
END MC811260 

ara ese qes MC811270 
SUBROUTINE SORT3(T,N,MXX) MC811280 

C --- INPLACE SORT USING SHELL ALGORITHM *******+* MG811290 
C --- SORTS ON 1ST 2 COLS. OF T & DOES SAME REORDERING ON 3RD°COLUMN MC811300 
INTEGER*2 T(MXX,3), ITEMP MC811310 
INTEGER GAP MC811320 
LOGICAL EXCH MC811330 

C MC811340 
GAP=(N/2) MC811350 

Ð IF (GAP. EQ. 0) GO TO 500 MC811360 
10 CONTINUE MG 85 20 
EXCH=. FALSE. MC811380 

K=N-GAP MC811390 

DO 200 I=1,K MC811400 

KK=I+GAP MC811410 

IE(T(1,1) GIMRGKK 1) OR: MC811420 

* (116 1). EO. DEKK 1) AND TEM? BEL DOEK, 2)) ) THEN MCS 1450 

au pO TO MC811440 

IT2=T(1,2) MC811450 

IT) MC811460 

T(I,1)=T(KK,1) MC811470 

T(1,2)=T(KK,2) MC811480 

TT Som RD MC811490 

TCKK,1)=IT1 MC811500 

T(KK,2)=IT2 MC811510 

T(KK,3)=IT3 MC811520 

EXCH=. TRUE. MES 11530 

ENDIF MC811540 

200 CONTINUE MC811550 
IF (EXCH) GO TO 10 MC811560 
GAP=(GAP/2) MC811570 

GO TO 5 MC811580 

500 CONTINUE MES 11590 
RETURN MC811600 

END MC811610 

C tedededede de dede dede dede deseo MC8 11620 
SUBROUTINE SUMBK( BKTBL,NBK,MXX,SINV,SY,INV,Y,MXY) MC811630 

C --- CREATE AGGREGATED ARRAYS INV,Y FROM CELL & 3RD DIM. INFO. IN BKTBLMC811640 
REAL INV(MXX,MXY) ,Y(MXX,MXY), SINV(MXX,MXY), SY(MXX ,MXY) MC811650 
INTEGER*2 BKTBL(MXX, 3) MC811660 
REAL*8 TINV,TY MC811670 

IP=0 MC811680 

Ils-1 MC811690 

I2s-1 MC811700 
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TINV=0 MC811. 


TY=0 MC811. 

DO 10 I71,NBK MC811, 
IF(BKTBLCI, 1). NE. 11 OR ERTRNÇI,2)- NES 12) EE MC811, 

Ee = CHANGE OF CELL,IW IDENTIFIERS MC811, 
IP=1P+1 MC811/ 

I1=BKTBL(1,1) MC811, 

I2=BKTBL(I, 2) MC811, 

DO 15 J=1,MXY MC811, 

INVCIP,J)=0 MC8118 

YCIP,J)=0 MC8118 

15 CONTINUE MC8118 
BKTBL(IP,1)-11 MC8118 

BKTBL(IP,2)=12 MC8118 

ENDIF MC8118 

C ==- ACCUMULATE MC8118: 
I3=BKTBL(1,3) MC8118 

DO 20 J=1,MXY MC8118 
INVCIP,J)=INVCIP,J)+SINVCI3 ,J) MC8118 





Y(IP,J)= Y(IP,J)+ SY(13,J) MC8119: 
TINV=TINV+SINV(13,J) MC8119: 

TY= TY+ SY(13,J) MC8119) 

20 CONTINUE MC8119) 
10 CONTINUE MC8119: 


WRITE(6,%) '==== TOTAL INV,Y AFTER BREAKDOWN=",TINV,TY MC8119. 


C MC81191 
NBK=IP e MC8119' 

C --- FIX INV. ENTRIES LOWER THAN CORRESP. LOSSES MC8119% 
DO 40 I=1,NBK MC81194 

DO 30 J=1,MXY MC8120( 
TEGINVET J). LI. YET NS TC) MC8120: 

30 CONTINUE MC8120: 
40 CONTINUE MC8 120: 
END MC81204 
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C. ESTIMATION SUBROUTINES 


SUBROUTINE MC87BZ(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) MC800010 
C --- CONDUCTS FIRST FIVE ESTIMATION METHODS MC800020 
C MC800030 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC800040 
REAL XTB(MXX), VXTB(MXX), A(MXX) MC800050 
C MC800060 
CALL EBTSI(INV,Y,NRC,NYR,XTB, VXTB, XEB,MXX,MXY) MC800070 
PRINT *, 'COMPLETED EBTS1' MC800080 
C MC800090 
CALL EBTS2(INV,X,NRC,NXR,XTB,VXTB,XEB,MXX MXX) MC800100 
PRINT *, ‘COMPLETED EBTS2' MC800110 
C MC800120 
CALL EBOS1(INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) MC800130 
PRINT *, 'COMPLETED EBOS1' MC800140 
C MC800150 
CALL EBOS2( INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) MC800160 
PRINT *, 'COMPLETED EBOS2' MC800170 
C MC800180 
CALL EMTS(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) MC800190 
PRINT $, 'COMPLETED EMTS' MC800200 
C MC800210 
END MC800220 
C MC800230 
AS AA MC GOOD 3 0 
C MC800250 
SUBROUTINE EBTS1( INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) MC800260 
C --- TRANSFORMED SCALE, TIME INDEPENDENT VARIANCE METHOD MC800270 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC800280 
REAL XTB(MXX), VXTB(MXX) MC800290 
REAL MAXL,MINL,L,MAXCHI , MINCHI MC800300 
INTEGER T, VYR MC800310 
DATA AA/1.6835/, B1/-.8934/, B2/.9881/ MC800320 
MAXL= -1000.0 MC800330 
MINL= 1000.0 MC800340 
SUML= 0.0 - MC800350 
KLSUM=0 | MC800360 
MAXCHI= -1000.0 MC800370 
MINCHI= 1000.0 MC800380 
SUMCHI= 0.0 MC800390 
KSUM=0 MC800400 
SUMMAD=0. 0 MC800410 
KPSUM=0 MC800420 
WERE (11.32) ' MC800430 
WRITE(11,21)'EMP BAYES TRANS SCALE - TIME DEP VAR: ' MC800440 
WRITE(11,21)'MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' MC800450 
WRITE(11,28)'FRACTION CELLS','FRACTION MAD' MC800460 
WRITE(11,29)'VALID YR','K','WITH UNDERAGE' ,'FROM UNDERAGE ' , MAD' MC800470 
WRITE(11,30) MC800480 
C MC800490 
C --- LOOP THROUGH VALIDATION YEARS MC800500 
DO 280 VYR=1, NYR MC800510 
C --- LOOP THROUGH CELLS MC800520 
DO 260 IN=1, NRC MC800530 
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o O AN 


C 


200 


260 


T=0 
SUMXT=0 
SUMVAR20 
LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 
DO 200 IT=1, NYR 
IF(IT .NE. VYR) THEN 
IFCINVCIN,IT) .NE. 0) THEN 
=FTT( INV(IN,IT), YCIN,IT)) 
C=SORT(0. 5+INV( IN, IT)) 
XX=X+C* (3. 141592654/2.0) 
XT-X/C 
T=T+1 
SUMXT=SUMXT+XT 
IF(XX .LT. 1.001) XX=1.001 
VARX=AA*(XX**B1)*(XX-1)**B2 
IF(VARX .GT. 1.0) VARX=1.0 
VARXT=VARX/ (0. 5+INV( IN, IT) ) 
SUMVAR=SUMVAR+VARXT 
ENDIF 
ENDIF 
CONTINUE 
XTB( IN)=SUMXT/T 
VXTB( IN)=SUMVAR/T**2 
CONTINUE 


CONDUCT ALGORITHM TO FIND XEB 
CALL EBITER(NRC,XTB,VXTB,XEB,MXX,VYR) 


COMPUTE MEAN SQUARED ERROR 
CALL MSECINV,Y,NRC,NYR,VYR,XEB,L,MXX,MXY,KL) 
IFCL . LT. MINL) THEN 
MINL=L 
MINLK=KL 
MINLYR=VYR 
ELSE IF(L .GT. MAXL) THEN 
MAXL=L 
MAXLK=KL 
MAXLYR=VYR 
ENDIF 
SUML=SUML+L*KL 
KLSUM=KLSUM+KL 


C --- INVERT XEB TO ORIGINAL SCALE 


C 


CALL INVERT(NRC,XEB,MXX) 


C --- COMPUTE MAD AND CHI SQUARE 


Se 


CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 
FCELLU, FMADU, PMAD,KP) 

IF(CHI . LT. MINCHI) THEN 
MINCHI=CHI 
MNCHI K=K 
MNCHYR=VYR 

ELSE IF(CHI .GT. MAXCHI) THEN 
MAXCHI=CHI 
MXCHIK=K 
MXCHYR=VYR 
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MC8005: 


MC80087 
MC80088 
MC80089 
MC80090. 
MC80091] 
MC80092: 
MC80093: 
MC80094: 
MC80095% 
MC80096t 
MC8009 71 
MC80098t 
MC80099( 
MC801004 
MC8010 1( 
MC801021 
MC80103( 
MC80104( 
MC80105( 
MC80106( 
MC80107( 
MC80108( 
MC80109( 


ENDIF MC801100 
SUMCHI=SUMCHI+CHI*K MC801110 
KSUM=KSUM+K MC801120 
KPSUM=KPSUM+KP MC801130 
SUMMAD=SUMMAD+PMAD*KP MC801140 

WRITE( 11,31) VYR,KP,FCELLU,FMADU,PMAD MC801150 

280 CONTINUE MC801160 
G MC801170 
C --- WRITE RESULTS TO OUTPUT FILE MC801180 
AVGL=SUML/KLSUM MC801190 
AVGCHI=SUMCHI /KSUM MC801200 
AVGMAD=SUMMAD / KPSUM MC801210 
WRITE(11,19)'AVG MAD = ',AVGMAD MC 801220 
WRITE( 11,21) ‘CHI SQUARE (ORIG SCALE): | MC801230 
WRITE(11,26) MIN CHI = ",MINCHI,'K = ' ,MNCHIK, 'VALID YR = ' ,MNCHYRMC801240 
WRITE(11,26)'MAX CHI = „MAXCHI, 'K = ' ,MXCHIK, ' VALID YR = ',MXCHYRMC801250 
WRITE(11,26)' AVG CHI = É MC801260 
WRITE(11,21)'MEAN SQUARED ERROR (TRANS SCALE): ' MC801270 
WRITE( 11,25) MIN MSE = ',MINL,'K = ',MINLK,'VALID YR = ',MINLYR MC801280 
WRITE(11,25) MAX MSE = ",MAXL, 'K = ' ,MAXLK,' VALID YR = ",MAXLYR  MC801290 
WRITE( 11,27) AVG MSE = ' ,AVGL MC801300 

19 FORMAT(38X,A,F5.3) MC801310 
21 FORMAT(/1X,A) MC801320 
DS. BORMAT(IX,A F6.3,5X,A,13,5X,A,12) MC801330 
26 FORMAT(1X,A,F9.3,5X,A,13,5X,A,12) MC801340 
27 FORMAT(1X,A,F6.3/) MC801350 
28 FORMAT(17X,A,2X,A) MC801360 
29 FORMAT(1X,A,4X,A,3X,A,2X,A,3X,A) MC801370 
BO TFORATCIX,8('-'), 2X, 4C -'),2X,14('-'"), 2X, 13("-"5),2X,5("-"J) MC801380 
31 FORMAT(1X,15,18,8X,F5.3,10X,F5.3,6X,F5.3) MC801390 
32 FORMAT(1X,A) MC801400 
RETURN MC801410 

END MC801420 
MC801430 

See EET TE EIE EE He dese TE EE IE KIE EK IE IE KIE IE kkk Fe NE EE He WEE He Ie TE He ENE eI rede MOBO 1440 
C MC801450 
SUBROUTINE EBTS2(INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) MC801460 

C --- TRANSFORMED SCALE, TIME INDEPENDENT VARIANCE METHOD MC801470 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC801480 
REAL XTB(MXX), VXTB(MXX) MC801490 
REAL MAXL,MINL,L,MAXCHI ,MINCHI MC801500 
INTEGER T, VYR MC801510 
MAXL= -1000. 0 MC801520 
MINL= 1000.0 MC801530 
SUML= 0.0 MC801540 
KLSUM=0 MC801550 
MAXCHI= -1000. 0 MC801560 
MINCHI= 1000.0 MC801570 
SUMCHI= 0.0 MC801580 
KSUM=0 MC801590 
SUMMAD=0. 0 MC801600 
KPSUM=0 MC801610 
WRITE(11,21)'EMP BAYES TRANS SCALE - TIME INDEP VAR: ' MC801620 
WRITE(11,21)'MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' MC801630 
WRITE(11,28)'FRACTION CELLS','FRACTION MAD' MC801640 
WRITE(11,29)'VALID YR','K','WITH UNDERAGE','FROM UNDERAGE' ,'MAD' MC801650 


7 


C 
od 


C oos 


300 


360 


WRITE( 11, 30) 


LOOP THROUGH VALIDATION YEARS 
DO 380 VYR=1, NYR 
LOOP THROUGH CELLS 
DO 360 IN=1, NRC 
T=0 
SUMXT=0 
SUMXT2=0 
LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 
DO 300 IT=1, NYR 
IF(IT .NE. VYR) THEN 
IFCINVCIN,IT) .NE. 0) THEN 
X=FTTCINVCIN, IT), YCIN,IT)) 
XT=X/SQRT(O. 5+INVCIN,IT)) 
=T+1 
SUMXT=SUMXT+XT 
SUMXT2=SUMXT2+XT**2 
ENDIF 
ENDIF 
CONTINUE 
XTB( IN)=SUMXT/T 
VXTB( IN)=( (T*SUMXT2) -( SUMXT**2) ) /((T-1)*T#*2) 
CONTINUE 


CONDUCT ALGORITHM TO FIND XEB 
CALL EBITER(NRC,XTB,VXTB,XEB,MXX,VYR) 


COMPUTE MEAN SQUARED ERROR 
CALL MSECINV,Y,NRC,NYR,VYR,XEB,L,MXX,MXY,KL) 
IF(L .LT. MINL) THEN 
MINL=L 
MINLK=KL 
MINLYR=VYR 
ELSE IF(L .GT. MAXL) THEN 
MAXL=L 
MAXLK=KL 
MAXLYR=VYR 
ENDIF 
SUML=SUML+L*KL 
KLSUM=KLSUM+KL 


INVERT XEB TO ORIGINAL SCALE 
CALL INVERT(NRC,XEB,MXX) 


COMPUTE MAD AND CHI SQUARE 
CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 
FCELLU, FMADU, PMAD,KP) 
IF(CHI .LT. MINCHI) THEN 
MINCHI=CHI 
MNCHIK=K 
MNCHYR=VYR 
ELSE IF(CHI .GT. MAXCHI) THEN 
MAXCHI=CHI 
MXCHIK=K 
MXCHYR=VYR 
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MC8019 
MC8019: 
MC8019, 
MC8019: 


ENDIF MC802220 
SUMCHI=SUMCHI+CHI*K MC802230 
KSUM=KSUM+K MC802240 
KPSUM=KPSUM+KP d MC802250 
SUMMAD=SUMMAD+PMAD*KP MC802260 

WRITE(11,31) VYR,KP,FCELLU, FMADU, PMAD MC802270 

380 CONTINUE MC802280 

€ MC802290 
C --- WRITE OUTPUT TO FILE MC802300 
AVGL=SUML/KLSUM MC802310 
AVGCHI=SUMCHI/KSUM MC802320 
AVGMAD=SUMMAD/KPSUM MC802330 
WRITE(11,19)'AVG MAD = ' , AVGMAD MC802340 
WRITE(11,21)'CHI SQUARE (ORIG SCALE): ' MC802350 
WRITE(11,26)'MIN CHI = ',MINCHI,'K = ',MNCHIK,'VALID YR = ' ,MNCHYRMC802360 
WRITE(11,26)'MAX CHI = ',MAXCHI,'K = ',MXCHIK,'VALID YR = ' ,MXCHYRMC802370 
WRITE( 11,26) ' AVG CHI = ' ,AVGCHI MC802380 
WRITE(11,21)'MEAN SQUARED ERROR (TRANS SCALE): ' MC802390 
WRITE(11,25)'MIN MSE = ',MINL,'K = ',MINLK,'VALID YR = ',MINLYR MC802400 
Via 25) MAX MSE = ',MAXL,'K = ‘ MAXLK, VALID YR = ',MAXLYR MG802410 
WRITE(11,27)'AVG MSE = ',AVGL MC802420 

19 FORMAT(38X,A,F5. 3) MC802430 

21 FORMAT(/1X,A) MC802440 

5 MAORMAR(1X,A,E6. 3,5X,A,13,5X,4,12) MC802450 

26 PORMATCIX,A,F9.3,5X,A,13,5X,A,12) MC802460 

27 FORMAT(1X,A,F6. 3/) MC802470 

28 FORMAT(17X,A,2X,A) MC802480 

29 FORMAT(1X,A,4X,A,3X,A,2X,A,3X,A) MC802490 

30 FORMAT(1X,8('-'),2X,4(' -*'),2X,14(" -'),2X,13('-'),2X,5('-')) MC802500 

CONAN 15 TE, 8X F5.3  10X F5.3, 6X,F5. 3) MC802510 

RETURN MC802520 
END MC802530 
E MC802540 
Vee ede dee Ve Vee Vee Vee Ve dee Te Ve ese Vee ede ee Vee Ve Ve ede ede dee dee Vee Vee Te Vee Vee ee Ve Ve Tee ET Vee ee ee EE RTEMCROI 550 
E MC802560 
SUBROUTINE EBOSI(INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) MC802570 
C --- ORIGINAL SCALE, TIME DEPENDENT VARIANCE METHOD MC802580 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC802590 
REAL XTB(MXX), VXTB(MXX) MC802600 
REAL MAXCHI, MINCHI MC802610 
INTEGER T, VYR MC802620 
MAXCHI= -1000. 0 MC802630 
MINCHI= 1000. 0 MC802640 
SUMCHI= 0.0 MC802650 
KSUM=0 MC 802660 
SUMMAD=0. 0 MC802670 
KPSUM=0 MC802680 
WRITE(11,21)'EMP BAYES ORIG SCALE - TIME DEP VAR: ' MC802690 
WRITE(11,21)'MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' MC802700 
WRITE(11,28)'FRACTION CELLS','FRACTION MAD' MC802710 
WRITE(11,29)'VALID YR','K','WITH UNDERAGE' ,'FROM UNDERAGE' ,'MAD' MC802720 
WRTTE(1 30) MC802730 
E MC802740 
C --- LOOP THROUGH VALIDATION YEARS MC802750 
DO 480 VYR=1, NYR MC802760 
C --- LOOP THROUGH CELLS MC802770 


77 


AN io 


DO 460 IN71, NRC MC802 





T=0 MC802' 
SUMXT=0 MC802i 
SUMVAR=0 MC8021 

C --- LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) MC8021 
DO 400 IT=1, NYR MC802i 

IF(IT .NE. VYR) THEN MC802i 
IFCINVCIN,IT) .NE. 0) THEN MC8021 
PHATSV(IN,IT)/INV(IN,IT) MC802% 

T=T+1 MC802E 
SUMXT=SUMXT+PHAT MC802f 

IF(PHAT .GT. 0.0) THEN MC802£ 
SUMVAR=SUMVAR+PHAT*( 1-PHAT) /INV(IN, IT) MC802° 

ELSE MC802€ 
SUMVAR=SUMVAR+1/(INV(IN,IT)+1)**2 MC8025 

ENDIF > MC8025 

ENDIF MC802S 

ENDIF MC8029 

400 CONTINUE MC8029 
XTB(IN)=SUMXT/T MC8029 

VXTB( IN)=SUMVAR/T**2 MC8029 

460 CONTINUE MC8029 
MC8030 

--- CONDUCT ALGORITHM TO FIND XEB MC8030 
CALL EBITER(NRC,XTB, VXTB, XEB,MXX, VYR) MC8030 
MC8030 

--- COMPUTE MAD AND CHI SQUARE MC8030 
CALL OSMOE(INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, MC8030 

* FCELLU, FMADU, PMAD, KP) MC8030 
IF(CHI . LT. MINCHI) THEN MC8030 
MINCHI=CHI MC8030 

MNCHIK=K MC8030 
MNCHYR=VYR MC8031; 

ELSE IF(CHI .GT. MAXCHI) THEN MC8031: 
MAXCHI=CHI MC8031. 

MXCHIK=K MC8031: 
MXCHYR=VYR MC8031! 

ENDIF MC8031: 
SUMCHI=SUMCHI+CHI*K MC80311 
KSUM=KSUM+K MC8031. 
KPSUM=KPSUM+KP MC8031. 
SUMMAD=SUMMAD+PMAD*KP MC8031! 

WRITE( 11,31) VYR,KP,FCELLU, FMADU, PMAD MC8032( 

480 CONTINUE MC8032: 
MC8032 

C --- WRITE OUTPUT TO FILE MC8032: 
AVGCHI=SUMCHI /KSUM MC80324 
AVGMAD=SUMMAD/KPSUM MC8032:! 
WRITE(11,19)'AVG MAD = ' ,AVGMAD MC8032€ 
WRITE(11,21)'CHI SQUARE (ORIG SCALE): ' MC8032; 
WRITE(11,26)'MIN CHI = ',MINCHI,'K = ',MNCHIK,'VALID YR = ' ,MNCHYRMC8032¢ 
WRITE(11,26)'MAX CHI = ",MAXCHI,'K = ',MXCHIK, VALID YR = ' MXCHVRMC80320 
WRITE(11,27)'AVG CHI = ",AVGCHI MC8033( 
19 FORMAT(38X,A,F5. 3) MC80331 
21 FORMAT(/1X,A) MC80332 
26 FORMATO]JN A FO. 3, 5X,A 13,5X A 12) MC8033:) 
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27 FORMAT(1X,A,F9.3/) MC803340 
28 FORMAT(17X,A,2X,A) MC803350 
29 FORMAT(1X,A,4X,A,3X,A,2X,A,3X,A) MC803360 
op BOELA IS EN - ),2X,40'-"), 2x4 -'), e )a2X.50°=')) MC803370 
31 FORMAT(1X,15,18,8X,F5.3,10X,F5.3,6X,F5. 3) MC803380 
RETURN MC803390 

END MC803400 

C MC803410 
ee Ee IE EE EE IE EE EE ee EE ee EE EE kkkikkkikkkkkkkikkkrrMC8BO3 420 
6 MC803430 
SUBROUTINE EBOS2( INV,Y,NRC,NYR,XTB,VXTB,XEB,MXX,MXY) MC803440 

C --- ORIGINAL SCALE, TIME INDEPENDENT VARIANCE METHOD MC803450 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC803460 

REAL XTB(MXX), VXTB(MXX) MC803470 

REAL MAXCHI, MINCHI MC803480 
INTEGER T, VYR MC803490 
MAXCHI= -1000.0 MC803500 
MINCHI= 1000.0 MC803510 
SUMCHI= 0.0 MC803520 
KSUM=0 MC803530 
SUMMAD20. 0 MC803540 
KPSUM20 MC803550 
WRITE(11,21)'EMP BAYES ORIG SCALE - TIME INDEP VAR: ' MC803560 
WRITE( 11,21) ‘MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' MC803570 
WRITE(11,28)'FRACTION CELLS','FRACTION MAD' MC803580 
WRITE(11,29)'VALID YR','K','WITH UNDERAGE','FROM UNDERAGE','MAD' MC803590 
WRITE(11,30) 2 MC803600 

e MC803610 
C --- LOOP THROUGH VALIDATION YEARS MC803620 
DO 580 VYR=1, NYR MC803630 

C --- LOOP THROUGH CELLS MC803640 
DO 560 IN=1, NRC MC803650 

T=0 MC803660 

SUMXT=0 MC803670 
SUMVAR=0 MC803680 

SUMY=0 MC803690 
SUMINV=0 MC803700 

C --- LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) MC803710 
DO 500 IT=1, NYR MC803720 

IF(IT .NE. VYR) THEN MC803730 
IF(INV(IN,IT) .NE. 0) THEN MC803740 

PHAT=Y(IN, IT) /INV(IN, IT) MC803750 
SUMXT=SUMXT+PHAT MC803760 

SUMY=SUMY+Y( IN, IT) MC803770 
SUMINV=SUMINV+INV (IN, IT) MC803780 

T=T+1 MC803790 
SUMVAR=SUMVAR+1. 0/INV( IN, IT) MC803800 

ENDIF MC803810 

ENDIF MC803820 

500 CONTINUE MC803830 
XTB( IN)=SUMXT/T MC803840 

IF(SUMY .GT. 0.0) THEN MC803850 
PTILDE=SUMY / SUMINV MC803860 

VXTB( IN)=( PTILDE*( 1-PTILDE )*SUMVAR) /T**2 MC803870 

ELSE MC803880 

VXTB( IN)=SUMINV**SUMVAR/(( (1+SUMINV )***2 ) *To*2 ) MC803890 
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ENDIF MC 803: 





560 CONTINUE MC803: 

E MC803: 
C --- CONDUCT ALGORITHM TO FIND XEB MC803! 
CALL EBITER(NRC, XTB,VXTB,XEB ,MXX, VYR) MC803! 

C MC803: 
C --- COMPUTE MAD AND CHI SQUARE MC 803“ 
„CALL OSMOE( INV,Y,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, MC803¢ 
FCELLU, FMADU , PMAD ‚KP) MC803* 

IF(CHI .LT. MINCHI) THEN MC8034 
MINCHI=CHI MC804( 

MNCHIK=K ' MC804( 

MNCHYR=VYR MC804( 

ELSE IF(CHI .GT. MAXCHI) THEN MC804( 
MAXCHI=CHI MC 804( 

MXCHIK=K MC804( 

MXCHYR=VYR MC804( 

ENDIF MC804( 
SUMCHI=SUMCHI+CHI*K MC804( 
KSUM=KSUM+K MC804( 
KPSUM=KPSUM+KP MC8041 
SUMMAD=SUMMAD+PMAD*KP MC8041 
WRITE(11,31) VYR,KP,FCELLU, FMADU , PMAD MC804] 

580 CONTINUE MC8041 

ë MC8041 
C --- WRITE OUTPUT TO FILE MC8041 
AVGCHI=SUMCHI/KSUM MC8041 
AVGMAD=SUMMAD/KPSUM MC8041 
WRITE(11,19)' AVG MAD = ' , AVGMAD MC8041 
WRITE(11,21)'CHI SQUARE (ORIG SCALE): ' MC8041 
WRITE(11,26)'MIN CHI = ',MINCHI,'K = ' ,MNCHIK, 'VALID YR = ' ,MNCHYRMC8042 
WRITE(11,26)'MAX CHI = ' ,MAXCHI,'K = ',MXCHIK,' VALID YR = ' ,MXCHYRMC8042 
WRITE(11,27)'AVG CHI = ' ,AVGCHI MC8042 

19 FORMAT(38X,A,F5.3) MC8042 
21 FORMAT(/1X,A) MC8042 
26 FORMAT(IX,A,F9.3,5X,A,13,5X,A4,12) MC804.2; 
27 FORMAT(1X,A,F9.3/) MC804 2: 
28 FORMAT(17X,A,2X,A) 2 MC8042: 
29 FORMAT(1X,A,6X,A,3X,A,2X,A,3X,A) | MC8042: 
30 FORMAT( 1X,8( ' - ' ) „2X,&4( ' =") 22, 140 EE MC8042: 
31 FORMAT( 1X,15,18,8X,F5.3,10X,F5.3,6X,F5.3) MC8043' 
RETURN MC8043? 

END MC8043: 

$ MC8043; 
ede IE dek e de dek dee zk ese de ee de dere Ie Fe Ie He Te Te He se Te e de e de de de e te e de de de e de de de kee EI LA NN NANA NANA ese vk Ter He MC ROG 33 
C MC8043; 
SUBROUTINE EMTS(INV,Y,NRC,NYR,XTB,VXTB,XEB,A,MXX,MXY) MC8043! 

C --- EFRON-MORRIS METHOD MC8043' 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC8043i 

REAL XTB(MXX), VXTB(MXX), A(MXX) MC80434 

REAL MAXL,MINL,L,MAXCHI ,MINCHI MC80441 
INTEGER T, VYR MC8044, 

DATA AA/1.6835/, B1/-.8934/, B2/.9881/ MC8044: 
MAXL= -1000.0 MC8044: 
MINL= 1000.0 MC80444 
SUML= 0.0 MC8044! 


so 


C 
C 


aa 


aa 


200 


260 


KLSUM=0 

MAXCHI= -1000.0 

MINCHI= 1000.0 

SUMCHI= 0.0 

KSUM=0 

SUMMAD=0. 0 

KPSUM=0 

WRITE(11,21)'EFRON-MORRIS TRANS SCALE - TIME DEP VAR: ' 
WRITE(11,21)'MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' 
WRITE(11,28)'FRACTION CELLS','FRACTION MAD' 
WRITE(11,29)'VALID XR','K','WITH UNDERAGE','FROM UNDERAGE’ ,'MAD' 
WRITE( 11,30) 


LOOP THROUGH VALIDATION YEARS 
DO 280 VYR=1, NYR 
LOOP THROUGH CELLS 
DO 260 IN=1, NRC 
T=0 
SUMXT=0 
SUMVAR=0 
LOOP THROUGH YEARS OF DATA TO COMPUTE XTB AND VAR(XTB) 
DO 200 IT=1, NYR 
IF(IT .NE. VYR) THEN 
IF(INV(IN,IT) .NE. 0) THEN 
X-FTT(INV(IN,IT), X(IN,IT)) 
C=SQRT(0. 5+INV(IN,IT)) 
XX=X+C* (3. 141592654/2. 0) 
XT=X/C 
T=T+1 
SUMXT=SUMXT+XT 
IF(XX .LT. 1.001) XX=1.001 
VARX=AA*(XX**B1)*(XX-1)**B2 
IF(VARX .GT. 1.0) VARX=1. 0 
VARXT=VARX/(0. 5+INV(IN,IT)) 
SUMVAR=SUMVAR+VARXT 
ENDIF 
ENDIF 
CONTINUE 
XTB( IN)=SUMXT/T 
VXTB( IN)=SUMVAR/T**2 
CONTINUE 


CONDUCT ALGORITHM TO FIND XEB 
CALL EMITER(NRC,XTB,VXTB,XEB,A,MXX,VYR) 


COMPUTE MEAN SQUARED ERROR 
CALL MSE( INV,Y,NRC,NYR,VYR, XEB, L,MXX,MXY,KL) 
IF(L .LT. MINL) THEN 
MINL=L 
MINLK=KL 
MINLYR=VYR 
ELSE IF(L .GT. MAXL) THEN 
MAXL=L 
MAXLK=KL 
MAXLYR=VYR 
ENDIF 
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MC804460 
MC804470 
MC804480 
MC804490 
MC804500 
MC804510 
MC804520 
MC804530 
MC804540 
MC804550 
MC804560 
MC804570 
MC804580 
MC804590 
MC804600 
MC804610 
MC804620 
MC804630 
MC804640 
MC804650 
MC804660 
MC804670 
MC804680 
MC804690 
MC804700 
MC804710 
MC804720 
MC804730 
MC804 740 
MC804750 
MC804 760 
MC804770 
MC804780 
MC804790 
MC804800 
MC804810 
MC804820 
MC804830 
MC804840 
MC804850 
MC804860 
MC804870 
MC804880 
MC804890 
MC804900 
MC804910 
MC804920 
MC804930 
MC804940 
MC804950 
MC804960 
MC804970 
MC804980 
MC804990 
MC805000 
MC805010 


E 
C --- INVERT XEB TO ORIGINAL SCALE 
CALL INVERT(NRC,XEB,MXX) 
C 
C --- COMPUTE MAD AND CHI SQUARE 
„CALL OSMOE( INV, Y ,NRC,NYR,VYR,XEB,CHI,K,MXX,MXY, 
FCELLU, FMADU, PMAD ‚KP) 
IECCHT SET MINCHI) THEN 
MINCHI=CHI 
MNCHIK=K 
MNCHYR=VYR 
ELSE IF(CHI .GT. MAXCHI) THEN 
MAXCHI=CHI 
MXCHIK=K 
MXCHYR=VYR 
ENDIF 
SUMCHI=SUMCHI+CHI*K 
KSUM=KSUM+K 
KPSUM=KPSUM+KP 
SUMMAD=SUMMAD+PMAD*KP 
WRITE(11,31) VYR,KP,FCELLU,FMADU, PMAD 
280 CONTINUE 
C 
C --- WRITE OUTPUT TO FILE 
AVGL=SUML/KLSUM 
AVGCHI=SUMCHI / KSUM 
AVGMAD=SUMMAD /KPSUM 
WRITE(11,19)'AVG MAD = ' ,AVGMAD 
WRITE(11,21)'CHI SOUARE (ORIG SCALE): ' 
WRITE(11,26) MIN CHI = ",MINCHI, K = ' 
WRITE(11,26)'MAX CHI = ',MAXCHI, 'K 
WRITE(11,26)'AVG CHI = ' ,AVGCHI 
WRITE(11,21)'MEAN SQUARED ERROR (TRANS SCALE): ' 
WRITE( 11,25) MIN MSE = " MINL, 'k 2 ' MINLK, VALID YR 
WRITE(11,25)'MAX MSE = ',MAXL, 'K = ' ,MAXLK, 'VALID YR 
WRITE(11,27)'AVG MSE = ',AVGL 
19 FORMAT(38X,A,F5. 3) 
21 FORMAT(/1X,A) 
25. FORMATCIX A. F6. 3,5% AL TE EER 
26 EORMAT( IX A,F9.3 EX A 15 SA 12) 
27 FORMAT( 1X,A,F6. 3) 
28 FORMAT(17X,A,2X,A) 
29 "FORMAT(1X A £X,A,3X ENE dT 
30 FORMATC(1X, BC '-'),2X AC('sD. 2% dd ad Pa 
31 FORMAT(1X,15,18,8X,F5. 3,10X,F5. 3,6X,F5. 3) 
RETURN 
END 
© 
EE EE IE IE EIE EE EIE IE de À Ve e ele R IE IE II KIE EIE EIE EIE IE HEE TENE HIE TE HEE WEE MEME HEME HB WKAM C BO 552 
C 
SUBROUTINE EBITER(NRC,XTB,VXTB,XEB,MXX,VYR) 
C --- ITERATIVE ALGORITHM TO SOLVE FOR XEB 


SUML=SUML+L*KL 
KLSUM=KLSUM+KL 


REAL XTB(MXX), VXTB(MXX), XEB(MXX) 


INTEGER VYR 
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,MNCHIK, ' VALID 
' ,MXCHIK, ' VALID 


YR 
MK 


' ,MINLYR 
' MAXLVR 


zn) 


' MNCHVRMC8053: 
' MXCHVRMC80534 





MC8051 
MC8052. 
MC8052 
MC805 2: 


MC8053: 
MC8053€ 
MC8053/ 
MC80536 
MC8053S 


MC80553 
MC80554 
MC80555 
MC80556 
MC80557 


A=0 MC805580 
ITER-0 MC805590 

100 CONTINUE MC805600 
ITER=ITER+1 MC805610 
IF(1TER .GT. 100) PRINT $, 'EBITER GT 100' MC805620 

AO=A MC805630 
SUMALK=0 MC805640 

C MC805650 
C --- SUM THE ALPHAS MC805660 
DO 200 I=1,NRC MC805670 
SUMALK=SUMALK+1/(A+VXTB(1)) MC805680 

200 CONTINUE MC805690 

C MC805700 
=== COMPUTE XBB MC805710 
XBB=0 MC805720 

DO 300 I=1,NRC MC805730 
ALPHA=1/(A+VXTB(1)) MC805740 
GAMMA=ALPHA / SUMALK MC805750 
XBB=XBB+GAMMA*XTB(1) MC805760 

300 CONTINUE MC805770 

C MC805780 
C --- UPDATE VALUE OF A MC805790 
SUMNUM=0 MC805800 
SUMDEN=0 MC805810 

DO 400 I=1,NRC MC805820 
ALPHA=1/(A+VXTB(I)) MC805830 
SUMNUM=SUMNUM+ALPHA*(XTB(1)-XBB)**2 MC805840 
SUMDEN=SUMDEN+( (XTB( 1) -XBB)*"*2)*ALPHA**2 MC805850 

400 CONTINUE MC805860 
A=A-(NRC-1-SUMNUM) /SUMDEN MC805870 

IF(A . LE. 0) THEN MC805880 

A=0 MC805890 

GO TO 500 MC805900 

ENDIF MC805910 
IF(ABS(A-A0) .GT. 0.0001) GO TO 100 MC805920 

500 CONTINUE MC805930 

C MC805940 
C --- ITERATIONS CONVERGED, COMPUTE XEB MC805950 
DO 600 I-1,NRC MC805960 
=XTB(I) MC805970 

V=VXTB( I) MC805980 

XEB( I)=(A*X) /(A+V)+(V*XBB) / (A+V) MC805990 

600 CONTINUE MC806000 
RETURN MC806010 

END MC806020 

C MC806030 
Vee de eed dee dee Vee ee de e de de dee e de de dee e ele de ee de de dee de de de de de dee de de de e de de de de elek de de de dee de de de dee ee dee de de dee M CG OG OD 0 
C MC806050 
SUBROUTINE EMITER(NRC,XTB,VXTB,XEM,A,MXX,VYR) MC806060 

C --- ITERATIVE ALGORITHM TO SOLVE FOR XEB FOR EFRON-MORRIS METHOD MC806070 
REAL XTB(MXX), VXTB(MXX), XEM(MXX), A(MXX) MC806080 
INTEGER VYR NC806090 

DO 100 I=1,NRC MC806100 
ACI)=O MC806110 
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100 


200 


300 
311 


399 


400 


444 


CONTINUE 


SUM THE ALPHAS 

SUMALK=0 

DO 200 I=1,NRC 
SUMALK=SUMALK+1/(ACI)+VXTB(I)) 

CONTINUE 


COMPUTE XHAT 

XHAT-0 

DO 300 I=1,NRC 
ALPHA=1/(A(I)+VXTB(1)) 
GAMMA=ALPHA / SUMALK 
XHAT=XHAT+GAMMA*XTB( 1) 

CONTINUE 

XHATP=XHAT 

I=1 

AP=A(I) 

S=(XTB( 1) -XHAT)**2 


COMPUTE SN AND SD 
SN=0 
SD=0 
DO 400 J=1,NRC 
TES NE. T) THEN 
DEN=( A( J )+VXTB( J) )**2 
SN=SN+((XTB(J)-XHAT)**2-VXTB(J) ) /DEN 
SD=SD+1/DEN 
ENDIF 
CONTINUE 


NEWTON-RAPHSON ITERATIONS TO SOLVE FOR A 
AD=A( I )+VXTB( I ) 
GNUM=S - 3* VXTB ( I ) +SN*AD**2 
GDEN=3+SD-AD***2 
GPRM=2*AD** SN /GDEN- 2**GNUM*AD:*SD/GDEN**2 
G=GNUM/GDEN 
ACI)=ACI)-(CACI)-G)/(1-GPRM)) 
IF(A(I) .LE. 0.) THEN 
A(1)=0. 0 
I=I+1 
IF(I . LE. NRC) THEN 
GO TO 333 
ELSE 
GO TO 555 
ENDIF 
ENDIF 
IF(ABS(A(I)-AP) .LE. 0.0001) THEN 
I=I+1 
IF(I .LE. NRC) THEN 
GO TO 333 
ELSE 
GO TO 555 
ENDIF 
ELSE 
AP=A(I) 
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MC8065: 
MC8065: 
MC8065: 
MC80652 
MC80655 


GO TO 444 MC806680 

ENDIF MC806690 

C MC806700 
G --- TEST FOR CONVERGENCE: ABS(S-SP) LT EPSILON MC806710 
555 SUMALK=0 MC806720 
DO 600 J=1,NRC MC806730 
SUMALK=SUMALK+1/(A(J)+VXTB( J)) MC806740 

600 CONTINUE > MC806750 
XHAT=0 MC806760 

DO 700 J=1,NRC MC806770 
ALPHA=1/( AC J)+VXTB(J)) MC806780 
GAMMA=ALPHA /SUMALK MC806790 
XHAT=XHAT+GAMMA*XTB(J) MC806800 

700 CONTINUE MC806810 
DO 800 J=1,NRC MC806820 
S=(XTB(J)-XHAT)**2 Va MC806830 
SP=(XTB(J)-XHATP)**2 MC806840 
IPQABSOSE5SP) .GT. 0.0001) GO TO 311 MC806850 

800 CONTINUE MC806860 

C MC806870 
C --- ITERATIONS CONVERGED, COMPUTE XEM MC806880 
DO 950 K=1,NRC MC806890 

SD=0 MC806900 

DO 900 J=1,NRC MC806910 

IF(J .NE. K) SD=SD+1/(A(J)+VXTB( J))**2 MC806920 

900 CONTINUE MC806930 
AD=A(K)+VXTB(K) MC806940 
DSTAR=3+( AD***2 )*SD MC806950 

B=( 1. -4. /DSTAR)*VXTB(K) /AD MC806960 

IF(B .GT. 1.0) B=1.0 MC806970 

IF(B .LT. 0.0) B=0.0 MC806980 
XEM(K)=XHAT+( 1-B)*(XTB(K) -XHAT) MC806990 

950 CONTINUE MC807000 
END MC807010 

C MC807020 
A A A AA RA A (01-10 F Á 0 kc 1 9) 
C MC807040 
SUBROUTINE MSE(INV,Y,NRC,NYR,VYR,XEB,L,MXX,MXY,KL) MC807050 

C --- COMPUTES MEAN SQUARED ERROR MOE MC807060 
REAL INV(MXX,MXY), Y(MXX,MXY), XEB(MXX) MC807070 

REAL L,MU MC807080 
INTEGER VYR MC807090 
SUMSE=0 MC807100 

KL=0 MC807110 

DO 100 I=1,NRC MC807120 

EIN LIN YR2 O @ GT.¿0. 0) THUEN MC807130 
=FTTCINVCI,VYR),YCI,VYR)) MC807140 
MU=X/SORT(0. 5+INV(1, VYR)) MC807150 
SUMSE=SUMSE+(XEB( 1) -MU)**2 MC807160 

KL=KL+1 MC807170 

ENDIF MC807180 

100 CONTINUE MC807190 
L=SUMSE /KL MC807200 
RETURN MC807210 

END MC807220 

C MC807230 
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Heme Here deve e ke kee ee ke e ke ke NA ANNA KA AA ANN NN AAN AN AN AK ANNAN AA AAN AAN A AA NA AZN N Fee He ever HMC BO 7 Z 


MC8072 

SUBROUTINE OSMOE( INV,Y,NRC ,NYR,VYR,XEB,CHI,K,MXX,MXY, MC807 2 

* FCELLU,FMADU, PMAD,KP) MC8072 

C =-- COMPUTES MAD AND CHI SOUARE MOES MC8072 
REAL INV(MXX ,MXY) , Y(MXX,MXY), XEB(MXX) MC8072 
INTEGER VYR MC80 78 
CHI=0. 0 MC8073 

K=0 MC8073 
SUMPM0=0. 0 MC8073 
SUMPMU=0. 0 MC8073 
KPMO=0 MC8073 
KPMU=0 MC8073 

DO 100 I71,NRC MC8073 
P=XEB(I) MC8073 
E=P*INV(I,VYR) MC8073 
A=Y(I,VYR) MC8074 

C MC8074 
C see COMPUTE MAD FOR THIS CELL MC8074 
IFCINVCI, VIR) .GT. 0.0) THEN MC8074 
PA7A/INV(I,VXR) MC8074 

IF(P .GT. PA) THEN MC8074. 
SUMPMO=SUMPMO+(P-PA) MC8074 

KPMO=KPMO+1 MC8074. 

ELSE MC8074: 
SUMPMU=SUMPMU+( PA-P ) MC8074. 

KPMU=KPMU+1 3 MC8075: 

ENDIF ` MC8075: 

ENDIF MC8075' 

C MC8075 
C ses COMPUTE CHI SQUARE FOR THIS CELL MC8075: 
IF(E «NE. 0.0 AND ) Paani 7.0) THEN MC8075 

K=K+1 MC8075: 
CHI=CHI+((A-E)*%*2)/(E*(1-P)) MC8075 

ENDIF MC8075:/ 

100 CONTINUE MC8075' 

C MC8076t 
C --- COMPUTE WEIGHTED AVERAGES MC8076; 
KP=KPMO+KPMU MC8076 
FCELLU=REAL(KPMU)/REAL(KP) MC8076: 
FMADU=SUMPMU/( SUMPMU+SUMPMO) MC80764 
PMAD=( SUMPMU+SUMPMO) /KP MC8076! 
RETURN MC8076t 

END MC8076, 

C MC8076t 
dedede dede dede de dede dede dede dede dede de dede dede dede dede de dede dede dede dede dede dede de dede dede de de de dede dede dede de dede dede dede dede de dede dede deite eMOBO7 6€ 
C MC8077( 
SUBROUTINE INVERT(NRC ,XEB,MXX) MC80771 

C --- INVERT XEB TO ORIGINAL SCALE MC80772 
REAL XEB(MXX) MC8 0776 

DO 100 I=1, NRC MC80774 

P=0, 5*( 1+SIN( XEB( I ) ) ) MC8077: 

TEM E ET O: O) SHEN MC8077€ 

P=0.0 MC8077 1 

ELSE IF (P .GT. 1.0) THEN MC8077€ 

P=1.0 MC8077$ 
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ENDIF MC807800 
XEB(I)=P MC807810 

100 CONTINUE MC807820 
RETURN MC807830 

END MC807840 

C MC807850 
EE EIE EIE II EK Je Je ae ge ge ge ge ze ge Je ze ge Je fe CK IKK KET CK $e Je Je ge Je Je ge Je içe ge $e e 3e 39 e e 3e zM Ç 8 O 7 8 60 
G MC807870 
FUNCTION FTT(INV,Y) MC807880 

C --- CONDUCTS FREMAN-TUKEY TRANSFORM MC807890 
REAL INV,Y MC807900 

TEMP =-1. + 2.*Y/(1.+INV) MC807910 
TEMP1=-1. + 2.*(1.+Y)/(1.+INV) MC807920 
IF(ABS(TEMP). GT. 1 .OR. ABS(TEMP1).GT. 1) THEN MC807930 
WRITE(6,*) 'FTT ERROR INV,Y=' INV, Y,TEMP,TEMP1 MC807940 

FTT=1 MC807950 

RETURN MC807960 

ENDIF MC807970 
FTT=SQRT(. 5+INV)*. 5*(ASINCTEMP) + ASIN(TEMP1) ) MC807980 

END MC807990 
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D. VECTOR METHOD SUBROUTINE 


* 


SUBROUTINE MC87V(INV,Y,MXX,NYR,NRC,XTBJI,DELTA,X,XVYR,VYRINV, 
VYRY,BSTAR,S,GAMMA,XBBJ,EVAL,MXP,MXK,BKTBL,NBK,NSC,NCSR, ISFLAG) MC8000 


C --- VECTOR METHOD 


80 
90 


REAL INV(MXX,NYR), Y(MXX,NYR) 

REAL XTBJI(MXP,MXK), DELTA(MXP,MXK), X(MXP,MXK) 
REAL XVYR(MXP,MXK), VYRINV(MXP,MXK), VYRY(MXP,MXK) 
REAL BSTAR(MXP,MXP), S(MXP,MXP), GAMMA(MXP,MXP) 
REAL XBBJ(MXP), EVAL(MXP) 

INTEGER*2 BKTBL(MXX, 3) 


REAL MAXL,MINL,L,MAXCHI ,MINCHI ,MO,MU, MAD 
INIEGER T, YYR TE 


MAXL= -1000.0 

MINL= 1000. 0 

SUML= 0.0 

KPSUM=0 

MAXCHI= -1000. 0 

MINCHI= 1000.0 

SUMCHI= 0.0 

KCSUM=0 
WVRITEF 32) 3 
WRITE(11,21)' EMP BAYES TRANS SCALE - VECTOR CASE: ' 
IF (ISFLAG .EQ. 1) THEN 

P=NSC 
WRITE(11,21)'VECTOR IS BY SERVICE COMPONENT ' 
ELSE 
P=NCSR 

WRITE(11,21)'VECTOR IS BY COMMISSIONING SOURCE’ 
ENDIF 
WRITE(11,22) 'K=' ,NRC, 'P=',P, 'KP= ',(NRC*P) 
WRITE(11,21)'MEAN ABSOLUTE DEVIATION (ORIG SCALE): ' 
WRITE(11,28)'FRACTION CELLS','FRACTION MAD' 


WRITE(11,29)'VALID YR','KP','WITH UNDERAGE' ,' FROM UNDERAGE ' , 'MAD' MC80034 


WRITE(11,30) 

K=NRC 

IF(K .LE. (P+2)) THEN 
WRITE(6,*)'*** ERROR IN VECTOR CASE: P+2 GT K ***' 
STOP 

ENDIF 


CONDUCT VALIDATION E 
DO 999 VYR=1, NYR 
DO 90 J=1,MXP 
DO 80 I=1,MXK 
MOBIC) 1) —020 
DELTA( J,1)=0. 0 
XVYR(J,1)=0. 0 
CONTINUE 
CONTINUE 
KMKG=BKTBL( 1,1) 
NRC=1 


--- LOOP THROUGH CELLS IN VECTOR FORM 
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MC800C 
MC800C 


MC8002. 
MC8002: 


MC80035] 
MC8003t 
MC8 003: 


100 
110 


C --- 


130 


200 


210 


DO 130 


IF(BKTBL(1,1) .NE. KMKG) NRC=NRC+1 


I=1,NBK 


DO 100 J=1,P 

IF(BKTBL(1,2) .EQ. J) THEN 
JP=J 

GO TO 110 


ENDI 


P 


CONTINUE i 


WRITE(6,*) '*** ERROR IN P VECTOR ASSIGNMENT **x' 


T=0 


SUMXT=0 
SUMVAR=0 


LOOP THROUGH YEARS OF DATA TO SOLVE FOR XTB AND VAR(XTB) 


DO 120 
IF(I 


TE TINV( T, IT) ONES 0) THEN 
FIJEFTIVGINEMT; IT), Y(I,IT)) 


EN 
ENDI 


IT-1, NVR 
T .NE. VYR) THEN 


C=0. 5+INV(I, IT) 
XT=XIJ/SQRT(C) 
SUMXT=SUMXT+XT 
SUMVAR=SUMVAR+1/C 
T=T+1 

DIF 

E 


CONTINUE 


XTBJI( 


STORE VARIANCE MATRIX IN DELTA MATRIX (TEMPORARY) 
DELTA( JP , NRC )=SUMVAR/T**2 
GET VALIDATION YEAR ESTIMATE, INVENTORY AND ATTRITION INFO 


JP,NRC)2SUMXT/T 


IFC INV Gi (KT GT 050) THEN 


XIJ=FTTV(INV(1,VYR), Y(I,VYR)) 


XT=XIJ/SQRTCO. 5+INV(I,VYR)) 
XVYR( JP ,NRC )=XT 


ENDIF 


VYRINV(JP,NRC)=INV(I,VYR) 


VYRY( JP ,NRC)=Y(1,VYR) 


KMKG=B 
CONTINUE 
TEC NE 


WRITE(6,*) '*%* ERROR IN VECTOR CASE: K ONE NRC ***' 


ENDIF 


COMPUTE 
DO 210 
SUMXTB 


IS TEA 


. NRC) THEN 


XBB SUB J 
J=1,P 
=0. 0 


DO 200 I=1,K 


SUMXTB=SUMXTB+XTBJI(J,I) 


CONTINUE 
XBBJ(J)-SUMXTB/K 


CONTINUE 


COMPUTE X SUB JI MATRIX, MAKE A COPY IN DELTA MATRIX (TEMPORARY) 


DO 230 


J=1,P 


DO 220 I=1,K 


de )—CATBII( J ,1)-XBBJCJ))*SQRTCDELTA(J,1)) 


DELTA(J,1)=X(J,1) 
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MC800540 
MC800550 
MC800560 
MC800570 
MC800580 
MC800590 
MC800600 
MC800610 
MC800620 
MC800630 
MC800640 
MC800650 
MC800660 
MC800670 
MC800680 
MC800690 
MC800700 
MC800710 
MC800720 
MC800730 
MC800740 
MC800750 
MC800760 
MC800770 
MC800780 
MC800790 
MC800800 
MC800810 
MC800820 
MC800830 
MC800840 
MC800850 
MC800860 
MC800870 
MC800880 
MC800890 
MC800900 
MC800910 
MC800920 
MC800930 
MC800940 
MC800950 
MC800960 
MC800970 
MC800980 
MC800990 
MC801000 
MC801010 
MC801020 
MC801030 
MC801040 
MC801050 
MC801060 
MC801070 
MC801080 
MC801090 


220 
230 


aa 


OO 262.0) 


240 


250 
260 


270 
280 


290 
300 


310 
320 


C --- 


CONTINUE 
CONTINUE 


COMPUTE S MATRIX 
CALL MXYTF(P,K,X MXP jE, ÉIER 


DO EIGENANALYSIS OF S 
PUT EIGENVALUES INTO EVAL, EIGENVECTORS INTO GAMMA 
CALL EVCSF(P,S,MXP,EVAL, GAMMA, MXP) 


CREATE ESTAR INVERSE 
KP2=K-P-2 
DO 240 J=1,P 
IF(EVAL(J) .LT. KP2) THEN 
EVAL( J)=KP2 
ENDIF 
EVAL(J)=1. 0/EVAL(J) 
CONTINUE 
DO 260 I=1,P 
DO 250 J=1,P 
IF(I .EQ. J) THEN 
BSTAR(1,J)=EVAL(J) 
ELSE 
BSTAR(I,J)=0. 0 
ENDIF 
CONTINUE 
CONTINUE 


CREATE BSTAR = I - (K-P-2) S TILDE INVERSE 
CALL MRRRR(P,P,GAMMA,MXP,P,P,BSTAR,MXP,P,P,S,MXP) 
CALL MXYTF(P,P,S,MXP,P,P,GAMMA,MXP,P,P,BSTAR,MXP) 
DO 280 I=1,P 
DO 270 J=1,P 
BSTAR(1,J)=KP2*BSTAR(I,J) 
TE(T Do J) THEN 
BSTAR(1,J)21. O-BSTAR(1,J) 
ELSE 
BSTAR(1,J)20. O-BSTAR(1,J) 
ENDIF 
CONTINUE 
CONTINUE 


COMPUTE DELTA SUB JI 
DO 300 J=1,P 
DO 290 1=1,K 
X(J,1)=XTBJI(J,1)-XBBJ(J) 
CONTINUE 
CONTINUE 
CALL MRRRR(P,P,BSTAR,MXP,P,K,X,MXP,P,K,XTBJI,MXP) 
DO 320 J=1,P 
DO 310 I=1,K 
DELTA(J, 1)=XBBJ(J)+XTBJI(J,1) 
CONTINUE 
CONTINUE 


COMPUTE MSE 
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MC8012. 
MC8012. 
MC8012' 
MC8013! 
MC8013: 
MC8013: 
MC8013. 
MC8013: 
MC8013: 


330 
340 


C --- 


350 
360 


C --- 


KP=0 
SUMSE=0. 0 
DO 340 J=1,P 
DO 330 I=1,K 
IF(INV(I,VYR) .GT. 0.0) THEN 
SUMSE=SUMSE+(DELTA(J,1)-XVYR(J,1))**2 
KP=KP+1 
ENDIF 
CONTINUE 
CONTINUE 
L=SUMSE /KP 
IF(L .LT. MINL) THEN 
MINL=L 
MINLKP=KP 
MINLYR=VYR 
ELSE IF(L .GT. MAXL) THEN 
MAXL=L 
MAXLKP=KP 
MAXLYR=VYR 
ENDIF 
SUML=SUML+L*KP 
KPSUM=KPSUM+KP 


INVERT DELTA SUB JI BACK TO ORIGINAL SCALE 
MOMJOONTJET,P 
DO 350 I=1,K 
PHAT=0. 5*(1+SIN(DELTA(J,I))) 
TE (PHAT .LT. 0.0) RHEN 


PHAT=0. 0 
ELSESIE (PHAT SGT. 180), MEN 
PHAT=1.0 
ENDIF 
DELTA(J,I)-PHAT 
CONTINUE 
CONTINUE 
COMPUTE CHI SOUARE AND MAD 
CHI=0. 0 
KCHI=0 
SUMPMO=0. 0 
SUMPMU=0. 0 
KPMO=0 
KPMU=0 


DO 410 J=1,P 
DO 400 I=1,K 
PHAT=DELTA( J, I) 
E=PHAT*VYRINV( J, 1) 
A=VYRY(J,1) 
IF(VYRINV(J,1) .GT. 0.0) THEN 

PACT=A/VYRINV(J,1) 

IF(PHAT .GT. PACT) THEN 
SUMPMO=SUMPMO+( PHAT-PACT) 
KPMO=KPMO+1 

ELSE 
. SUMPMU=SUMPMU+( PACT-PHAT) 
KPMU=KPMU+1 
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MC801660 
MC801670 
MC801680 
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MC801700 
MC801710 
MC801720 
MC801730 
MC801740 
MC801750 
MC801760 
MC801770 
MC801780 
MC801790 
MC801800 
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MC801860 
MC801870 
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MC801890 
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MC801920 
MC801930 
MC801940 
MC801950 
MC801960 
MC801970 
MC801980 
MC801990 
MC802000 
MC802010 
MC802020 
MC802030 
MC802040 
MC802050 
MC802060 
MC802070 
MC802080 
MC802090 
MC802100 
MC802110 
MC802120 
MC802130 
MC802140 
MC802150 
MC802160 
MC802170 
MC802180 
MC802190 
MC802200 
MC802210 


ENDIF MC8022 


ENDIF MC8022 
IF(E.NE. 0.0 . AND. PHAT.NE. 1.0) THEN MC8022 
KCHI=KCHI+1 MC8022 
CHI=CHI+((A-E)**2)/(E*(1-PHAT)) MC8022 

ENDIF MC8022 

400 CONTINUE MC8022 
410 CONTINUE MC8022 
KMAD=KPMO+KPMU MC8023 
FCELLU=REAL( KPMU) /REAL( KMAD) MC8023 
FMADU=SUMPMU/ ( SUMPMU+SUMPMO ) MC8023. 
PMAD=( SUMPMU+SUMPMO) /KMAD MC8023 
IF(CHI . LT. MINCHI) THEN MC8023 
MINCHI=CHI MC8023. 
MNCHIK=KCHI MC8023 
MNCHYR=VYR > MC8023 
ELSE IF(CHI .GT. MAXCHI) THEN MC8023 
MAXCHI=CHI MC8023 
MXCHIK=KCHI MC8024. 
MXCHYR=VYR MC8024. 
ENDIF MC8024. 
SUMCHI=SUMCHI+CHI*KCHI MC8024: 
KCSUM=KCSUM+KCHI MC8024. 
WRITE(11,31) VYR,KMAD,FCELLU, FMADU, PMAD MC8024, 
999 CONTINUE MC8024: 
AVGL=SUML/KPSUM MC8024. 
AVGCHI=SUMCHI /KCSUM MC8024: 
WRITE(11,21)'CHI SQUARE (ORIG SCALE): ' MC8024 
WRITE(11,26)'MIN CHI = ',MINCHI,'KP = ' ,MNCHIK, MC80251 
ve ' VALID YR = ' ,MNCHYR MC8025 
WRITE(11,26)'MAX CHI = ',MAXCHI,'KP = ',MXCHIK, MC8025: 
* ' VALID YR = ' ,MXCHYR MC8025: 
WRITE(11,26)'AVG CHI = ' ,AVGCHI MC8025: 


WRITE(11,21)'MEAN SQUARED ERROR (TRANS SCALE): ' MC8025! 
WRITE(11,25)'MIN MSE ' MINL,'KP = ",MINLKP,' VALID YR = ',MINLYR MC8025t 
WRITE(11,25)'MAX MSE = ',MAXL,'KP = ",MAXLKP,' VALID YR = ",MAXLYR MC8025' 







WRITE( 11,27) AVG MSE = ',AVGL MC80251 
21 FORMAT(/1X,A) MC8025$ 
22 -FORMAT SA 13.59) MC8026t] 
25 JKORMAMQGIX A F6.3. 5% AIS SR 2) MC8026; 
26 FORMAICIX,A,F9.3,5X,A,I3,S5K A 172) MC8026; 
27 FORMAT(1X,A,F6.3/) MC8026: 
28 SEORNAT( 17X A 2% A) MC80264 
29 FORMAT(1X,A,3X,A,3X,A,2X,A,3X,A) MC8026: 


30 FORMAT(1X,8("-"),2X,4('"-"),2X,14("-"),2X 180 3824511) MC8026€ 
31 FORMAT(1X,I5,18,8X,F5.3,10X,F5.3,6X,F5.3) MC80267 
32 FORMAT(1X,A) MC8026€ 


WRITE(6,*)'COMPLETED VECTOR CASE' MC8026¢ 
END MC8027( 
C MC80271 
Heese vere ve de dere de de Fe ve ve Me Fe se Me e de ele e dee de ee dee dee ee de de ee dee He de dee de MEME EK EE KERKE WEE TE VEN ee We ME Me Hee He M CBO? 7 Z 


C MC80278 
FUNCTION FTTV(INV,Y) MC80274 
C --- CONDUCTS FREMAN-TUKEY TRANSFORM MC80273 
REAL INV,Y MC8027€ 


TEMPS E Ed MC80277 
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Il 2 + (1. +Y)/( 1. +INV) 

IF(ABS(TEMP). GT. 1 .OR. ABS(TEMP1). GT. 1) THEN 
e e 'FTT ERROR INV,Y=' ,INV,Y,TEMP,TEMP1 
RETURN 

ENDIF 

Geer `` ITEE + ASIN(TEMP1)) 
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MC802780 
MC802790 
MC802800 
MC802810 
MC802820 
MC802830 
MC892840 
MC802850 


E. EXEC PROGRAM 


CP LINK MVS 103TmIUSTRR 

ACC 103 K 

detededede deseo dede dede dede dede de dede de de dede dede dede dede dede dede 

FIL * CLEAR 

FIL 01 K DSN F0968 MCOR87 DATA (RECFM FB LRECL 69 BLOCK 17940 
FIL 02 DISK MC87 TEMP 

FIL 06 81 (RECFM FBA LRECL 150 

FIL 25 DISK MCLASS PG15 (RECFM F LRECL 25 
FIL 27 DISK MCLASS PG17 (RECFM F LRECL 25 
FIL 29 DISK MCLASS PG19 (RECFM F LRECL 25 
FIL 30 DISK MCLASS PG20 (RECFM F LRECL 25 
FIL 31 DISK MCLASS PG21 (RECFM F LRECL 25 
FIL 32 DISK MCLASS PG22 (RECFM F LRECL 25 


&BEGSTACK 

30.0 /* AVG INV THRESHOLD T */ 

30 /* NO. CELLS THRESHOLD K */ 

13 /* MOS (ONLY 1) */ 

4 /* YCS (ONLY 1) */ 

15 /' GRADE (ONLY 1) */ 

3 128 /* NO. SVC COMPS AND ARRAY(1-REG,2-AUGREG,3-RES,4=1+2,5=ALL */ 
1 16 /* NO. COMM SRCS AND ARRAY(1-15, 16=ALL) 

1 /** 3RD DIMENSION (O=NONE, 1=SVC, 2=CS) 

SEND 


LOAD MC87 (START CLEAR 
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F. SAMPLE DATA FILE 


1 27 3 As 6 7 
15 0 12128 1 0102 
15 0, 3 TZS 5 0.42 
15 US R 6 0. 75 
15 £ 10 4 0715 
ÉS 1 11 1 0210 
15 0513. 10 2 0. 15 
15 0 272,53 1 0702 
15 01 23987 4 il 0. 02 
15 OA 6 0557 
15 0 4 2 10 2 9320 
15 ow 3 9 2 0. 07 
15 0 4 1 11 1 0. 10 
15 0 43 4 2 0.05 
15 0 43 10 1 ES 
15 Vz 3 0827 
15 OU. 2 io 2 0. 20 
15 OF kl ili 007 
15 OO 1 0202 
15 Lð 2 35 1 0.05 
15 5 1721 2 0. 20 
15 1 3 1 0. 05 
15 ES SW 1 OSLO 
15 Ns os 3 VIE 
15 18857 3 5 0230 
15 MD 3 2 0.65 
1:5 ae 35538 2 5 0.45 
15 IS Ð 0012 
15 VS 2 15 il 005 
15 j 1 T1 3 0.12 
15 Po bo 2 0. YY 
15 1 1 0502 
15 e Seed E 1 OO 
15 J 2:5 es) 2 012 
15 l 43 2 5 0. 40 
15 1 42 3 2 B22 
15 1 & 1 1 d OSOZ 
15 1 42 9 2 0.07 
15 I em 1 002 
15 Lasso 3° 7 5 025 
15 1 4310 l 0. 02 
15 | 243,73 z 0. 05 
15 1 4110 1 002 
15 jJ 3.12 1 0202 
15 1 42 12 1 0.07 


(remaining entries omitted) 


Column descriptions: 


1 - grade 5 - commissioning source 
27 = HOS 6 - number of records 
SAS 7 - total average inventory 


4 - service component 
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C ses 


NET 


PROGRAM TO CREATE INVENTORY DATA FILE BY GRADE 
PARAMETER (MXX=20000, MXY=10) 

CLASSIF. TABLE: GRADE, MOS, YCS, SVC, CS 

INTEGER*2 PTRTBL(MXX, 5), NRECS(MXX) 

REAL AINV(MXX) 

INTEGER TVPE,VCS,PG,MOS,SEX,CS,EDLV,SVC,MOS1,MOS2 ,RACE 
INTEGER DATA(MXV), SPG 

CHARACTER*7 CITLS 

DATA AINV/MXX*0. /, NRECS/MXX*O/ 


WRITE(5,*) 'ENTER PG' 
READ(5,*) SPG 
WRITE(5,*) 'PG TO USE=' ,SPG 


ICR=0 
NRC=0 
NG=0 
DO 10 1=1,999999 
READ( 1,100,END=999) TYPE,YCS,PG,MOS,SEX,CS, EDLV, SVC,MOS1,MOS2, 
* RACE , CITLS ,DATA 
ICR=ICR+1 


e ese 


NET 


10 


999 


20 
100 
101 


s 


C --- 


C --- 


CLASSIFY ALL RECORDS TYPE 0 

IF(TYPE.GT. 0) GO TO 999 

ADD NEW RECORD TO TABLE 

IF(PG. EQ. SPG) CALL ADDTBL(PG,MOS,YCS,SVC,CS, DATA,MXY, PTRTBL, 


* MXX, NRC,AINV,NRECS) 


IF( PG. EQ. SPG) NG=NG+1 
IF(MOD(ICR,5000).EQ.0) WRITE(6,*) ' ICR,NRC=",ICR,NRC 


CONTINUE 

CONTINUE 

WRITE(6,*) ' ' 

WRITE(6,*) 'TOTAL RECORDS READ =' ICR 
WRITE(6,*) 'TOTAL RECORDS ACCEPTED =' ,NG 

WRITE(6,*) 'TOTAL INVENTORY COMBINATIONS =',NRC 


DO 20 I=1,NRC 

WRITE(2,101) (PTRTBL(1,J),J=1,5), NRECS(I) ,AINV(I) 
CONTINUE 

FORMAT(312,13,11,12,211,213,11,47, 1X, 1014) 
FORMAT(12,14,13,12,13, 14, F7.2) 

END 


SUBROUTINE ADDTBL(PG,MOS,YCS,SVC,CS, DATA,MXY, PTRTBL,MXX, NRC, 
AINV,NRECS) 
SET INVENTORY POINTER FOR THIS ENTRY AND ACCUMULATE 
INTEGER*2 PTRTBL(MXX, 5), NRECS(MXX) 
REAL AINV(MXX) 
INTEGER YCS,PG,MOS,CS,SVC 
INTEGER DATA(MXX) 
MINV=GETINV(PTRTBL, MXX,NRC, PG,MOS,YCS,SVC,CS) 
IF(MINV .EQ. 0) THEN 
NEW COMBINATION 
NRC=NRC+1 
IF(NRC .GT. MXX) THEN 
WRITE(6,%) '*%% ERROR - TOO MANY INV. COMBINATIONS ' NRC 
STOP 
ENDIF 
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MC8002: 

MC8002: 
MC8002: 
MC8002: 
MC8002! 
MC8002¢ 
MC8002:1 


Moe a 
| 


| 


110 


MINV=NRC 
PTRTBL(MINV, 1)=PG 
PTRTBL(MINV, 2)=MOS 
PTRTBL(MINV, 3)=YCS 
PTRTBL(MINV, 4)=SVC 
PTRTBL(MINV, 5)=CS 
NRECS(MINV)20 

ENDIF 

AI=0 

DO 110 IT=1,MXY 

AI=AI + FLOAT(DATA(IT)) 

CONTINUE 

AINV(MINV)=AINV(MINV) + .25*AI/MXY 

NRECS(MINV)=NRECS(MINV) + 1 

END 


FUNCTION GETINV(PTRTBL, MXX,NRC, PG,MOS,YCS,SVC,CS) 
FIND LOCATION OF MATCHING INVENTORY ENTRY FOR A LOSS 
INTEGER*2 PTRTBL(MXX, 5) 
INTEGER YCS,PG,MOS,CS,SVC 
DO 10 I21,NRC 
IF(PTRTBL(I, 1) .EQ. PG .AND. 
PTRTBL(I, 2) .EQ. MOS .AND. 
PRRRBRCT, 3) BO. YCS AND: 
PTRTBL(I, 4) .EQ. SVC .AND. 
PTRTBL(I, 5) .EQ. CS ) THEN 


GETINV=I 
RETURN 
ENDIF 
10 CONTINUE 
GETINV=0 
END 


97 


MC800570 
MC800580 
MC800590 
MC800600 
MC800610 
MC800620 
MC800630 
MC800640 
MC800650 
MC800660 
MC800670 
MC800680 
MC800690 
MC800700 
MC800710 
MC800720 
MC800730 
MC800740 
MC800750 
MC800760 
MC800770 
MC800780 
MC800790 
MC800800 
MC800810 
MC800820 
MC800830 
MC800840 
MC800850 
MC800860 
MC800870 
MC800880 


APPENDIX C. SAMPLE OUTPUT “ 


A. GENERAL 

This appendix contains sample output from the computer program. A sample out- 
put for test cases one through 30 which use the first five estimation methods is shown 
in paragraph B. A sample output for the vector test cases is shown in paragraph C. 
“These examples show the output that is produced by the WRITE statements for file 
definition 11, e.g., WRITE(11,101). The program also contains several WRITE and 
PRINT statements that provide interactive information to the user via the terminal 
screen, e.g., WRITE(6,*), WRITE(5,*) and PRINT *. This interactive output is 


omitted. 


B. SAMPLE OUTPUT (TEST CASES 1-30) 


TEST CASE INPUT PARAMETERS: 

INVENTORY THRESHOLD= 30.0 THRESHOLD NO. OF CELLS= 30 
MOS= 13 YCS= 4 GRADE= 15 

SERVICE COMPONENTS= ime 3 

COMM SOURCES= l 2 34 53 6 7 8 9 10°21 12 13 LIE 


EXPANSION INFORMATION: 
ACTUAL NO: OF CELLS USED: 24 
MOS GROUP # 1 YCS'S USED= 
A 5 
LARGE MOS GROUP #1 YCS'S USED= 
4 5 


MAJOR MOS GROUP #1 YCS'S USED= 
4 5 


EMP BAYES TRANS SCALE - TIME DEP VAR: 


MEAN ABSOLUTE DEVIATION (ORIG SCALE): 
FRACTION CELLS FRACTION MAD 
VALID YR K WITH UNDERAGE FROM UNDERAGE MAD 


l 24 0.458 0.478 OZ 
2 24 02250 0.167 0: 099 
3 24 0.542 0. 441 0. 098 
4 24 0.355 (5 02089 
5 24 0.417 Ü 35 2 0. 072 
6 24 O1 25 02.053 0.082 
7 24 0. 208 0. 138 07099 
8 24 0.417 0.472 0 077 
9 24 0.655 0. 943 0. 181 
10 24 0.1333 07352 OI 

AVG MAD = 0.102 
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CHI SQUARE (ORIG SCALE): 


MIN CHI = 48.590 K = 24 VALID YR = 8 
MAX CHI = 329.334 K = 24 VALID YR = 9 
AVG CHI = 99. 791 

MEAN SQUARED ERROR (TRANS SCALE): 

MIN MSE = 0.033 K = 24 VALID YR = 4 
MAX MSE = 0.205 K= 24 VALID YR = 9 
AVG MSE = 0.079 


EMP BAYES TRANS SCALE - TIME INDEP VAR: 


MEAN ABSOLUTE DEVIATION (ORIG SCALE): 
FRACTION CELLS FRACTION MAD 
VALID YR K WITH UNDERAGE FROM UNDERAGE MAD 


1 24 0.458 0. 504 0.127 
2 24 0. 250 0. 204 0. 094 
> 24 0563 0.477 0.101 
4 24 0.575 0. 432 0. 071 
5 24 0.417 0.396 0.073 
6 24 0. 083 0. 061 0. 076 
Z 24 0.250 9.152 0.094 
8 24 0.458 03517 0.075 
9 24 0. 792 0.952 0. 183 
10 24 02875 0.961 0.115 
AVG MAD = 0.101 
CHI SQUARE (ORIG SCALE): 
MIN CHI = 45.452 K = 24 VALID YR = 6 
MAX CHI = 344, 445 K = 24 VALID YR = 9 
AVG CHI = 99. 284 
MEAN SQUARED ERROR (TRANS SCALE): 
MIN MSE = 0.036 K = 24 VALID YR = 4 
MAX MSE = 0.213 K = 24 VALID YR = 9 
AVG MSE = 0.078 


EMP BAYES ORIG SCALE - TIME DEP VAR: 


MEAN ABSOLUTE DEVIATION (ORIG SCALE): 
FRACTION CELLS FRACTION MAD 
VALID YR K WITH UNDERAGE FROM UNDERAGE MAD 


iL 24 0.458 0.478 02127 
2 24 0. 250 02195 0.097 
3 24 0. 542 0. 461 01095 
4 24 09998 07297 0. 070 
5 24 0.417 03375 0. 074 
6 24 0, 125 0.057 0.079 
7 24 0. 208 0. 161 0.102 
8 24 0. 417 0. 488 02079 
9 24 0.833 0.944 02185 
10 24 02533 0.952 05117 
AVG MAD = 0.102 
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CHI SQUARE (ORIG SCALE): 


MIN CHI = 49. 207 K= 24 VALID YR = 
MAX CHI = 340.035 K = 24 VALID YR = 
AVG CHI = 100. 985 


EMP BAYES ORIG SCALE - TIME INDEP VAR: 


MEAN ABSOLUTE DEVIATION (ORIG SCALE): 
FRACTION CELLS FRACTION MAD 
VALID YR K WITH UNDERAGE FROM UNDERAGE 


1 24 0.458 0.479 
2 24 0. 250. 0.196 
3 24 0.542 0.462 
4 24 0.335 01599 
Ð 24 0.417 0577 
6 24 0125 08056 
7 24 0. 208 0. 161 
8 24 0. 417 0. 490 
9 24 0:83. 0.945 
10 24 02833 0. 955 
AVG MAD = 
CHI SQUARE (ORIG SCALE): 
MIN CHI = 48. 836 K= 24 VALID YR = 
MAX CHI = 339939 5 K= 24 VALID YR = 
AVG CHI = 100. 960 


EFRON-MORRIS TRANS SCALE - TIME DEP VAR: 


MEAN ABSOLUTE DEVIATION (ORIG SCALE): 
FRACTION CELLS FRACTION MAD 
VALID YR K WITH UNDERAGE FROM UNDERAGE 


ll 24 0.458 0,475 
2 24 0. 250 0.172 
2 24 0. 583 0.451 
4 24 0.375 0 365 
5 24 0. 333 0. 348 
6 24 0. 083 0. 046 
7 24 0. 208 05155 
8 24 0. 458 0. 469 
2 24 0. 708 0.958 
10 24 0.875 0.960 
AVG MAD = 
CHI SQUARE (ORIG SCALE): 
MIN CHI = &6. 301 K= 24 VALID YR = 
MAX CHI = 340.712 K = 24 VALID YR = 
AVG CHI = 101251 
MEAN SQUARED ERROR (TRANS SCALE): 
MIN MSE = 0.042 K= 24 VALID YR = 4 
MAX MSE = 0.211 K = 24 VALID YR = 9 
AVG MSE = 0.080 
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C. SAMPLE OUTPUT (VECTOR TEST CASES) 


TEST CASE INPUT PARAMETERS: 

INVENTORY THRESHOLD= 30.0 THRESHOLD NO. OF CELLS= 30 
MOS= 151 YCS= 7 GRADE= 17 

SERVICE COMPONENTS= LIZ 

COMM SOURCES= ee 8 9 10-11-12 13 Gout? 


EXPANSION INFORMATION: 
ACTUAL NO. OF CELLS USED= 8 
MOS GROUP # 8 YCS'S USED= 
7 
LARGE MOS GROUP #3 YCS'S USED= 
7 
MAJOR MOS GROUP #2 YCS'S USED= 
7 
EMP BAYES TRANS SCALE - VECTOR CASE: 


VECTOR IS BY SERVICE COMPONENT 
K= 8 P= 3 KP= 24 


MEAN ABSOLUTE DEVIATION (ORIG SCALE): 
FRACTION CELLS FRACTION MAD 
VALID YR KP WITH UNDERAGE FROM UNDERAGE MAD 


1 24 0.375 03019 0.186 
2 24 0.458 OS Os d 79 
3 24 0. 458 0.417 0150 
4 24 0. 458 0. 422 AZ 
5 24 0. 542 0. 708 0. 146 
6 24 0292 Fl 0. 142 
7 24 0. 230 04193 02726 
8 24 02975 0.434 02092 
9 24 0.458 0: 705 05161 
10 24 O. On 912 07202 

CHI SQUARE (ORIG SCALE): 

MIN CHI = 278827 KP = 24 VALID YR = 8 

MAX CHI = 165.694 KP = 24 VALID YR = 10 

AVG CHI = 61. 025 

MEAN SQUARED ERROR (TRANS SCALE): 

MIN MSE = 0.089 KP = 24 VALID YR = 8 

MAX MSE = 0.483 KP = 24 VALID YR = 10 

AVG MSE = 0.229 
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